Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,231,849 T→A 5.6% H2Q (CAT→CAA)  cdd → cytidine/deoxycytidine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,231,8490TA 5.6% 230.0 / 3.5 72H2Q (CAT→CAA) cddcytidine/deoxycytidine deaminase
Reads supporting (aligned to +/- strand):  ref base T (41/27);  new base A (2/2);  total (43/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

AGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGC  >  NC_000913/2231702‑2231997
                                                                                                                                                   |                                                                                                                                                    
aGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATc                                                                                                                                                     >  2:374825/1‑149 (MQ=255)
         cTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTc                                                                                                                                            <  4:84059/149‑1 (MQ=255)
            aaaTTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGTCAACGCCATGCATCCACGTTTTCaaa                                                                                                                                         >  2:141425/1‑149 (MQ=255)
                 tGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGc                                                                                                                                     <  2:432097/148‑1 (MQ=255)
                     aTGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTg                                                                                                                                <  1:374825/149‑1 (MQ=255)
                               cATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGc                                                                                                                      <  1:371933/149‑1 (MQ=255)
                                 tttttGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTg                                                                                                                                >  7:264900/1‑137 (MQ=255)
                                 tttttGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTg                                                                                                                                <  8:264900/137‑1 (MQ=255)
                                 tttttGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCgg                                                                                                                    <  4:238908/149‑1 (MQ=255)
                                      gATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCGACAACGGGTTCGTAAACTGATATCCCATTACATGATTATGAGGCACCGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAAc                                                                                                               >  6:391933/1‑149 (MQ=255)
                                         gtatgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTg                                                                                                            >  1:25281/1‑149 (MQ=255)
                                          tatgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTg                                                                                                            >  5:70912/1‑148 (MQ=255)
                                            tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCGCATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCaa                                                                                                         >  7:306025/1‑149 (MQ=255)
                                            tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCaa                                                                                                         <  3:153964/149‑1 (MQ=255)
                                              tttCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAAt                                                                                                        <  6:221326/148‑1 (MQ=255)
                                                  aCGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGc                                                                                                    <  6:118434/148‑1 (MQ=255)
                                                  aCGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCa                                                                                                   <  2:25281/149‑1 (MQ=255)
                                                    gcgTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCACCTTGCGGATAACTTGCAATCTGCACt                                                                                                 >  3:9542/1‑149 (MQ=255)
                                                                taatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCACATTACATGATTATGAGGCAACGCCAGGCATCCCCGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCt                                                                                     >  6:247992/1‑149 (MQ=255)
                                                                  atGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTgg                                                                                   >  4:345468/1‑149 (MQ=255)
                                                                     agaTACAGATCACATATAAAGCCACAACCGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCCCGTTTACAAACCGCTATTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTgg                                                                                   >  1:381839/1‑146 (MQ=255)
                                                                          cAGATCGCATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACaa                                                                            <  8:306025/148‑1 (MQ=255)
                                                                          cAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCCTCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGATaa                                                                            <  7:84302/148‑1 (MQ=255)
                                                                           aGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCTGAGAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAg                                                                           >  4:326018/1‑148 (MQ=255)
                                                                           aGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGt                                                                          <  5:219988/149‑1 (MQ=255)
                                                                            gATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCtttt                                                                                                                                 >  7:399832/1‑93 (MQ=255)
                                                                            gATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCtttt                                                                                                                                 <  8:399832/93‑1 (MQ=255)
                                                                            gATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTa                                                                         <  3:345468/149‑1 (MQ=255)
                                                                                 caTATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCAAGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTcc                                                                    <  1:141425/149‑1 (MQ=255)
                                                                                    ataAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCa                                                                                 >  1:263896/1‑133 (MQ=255)
                                                                                    ataAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCa                                                                                 <  2:263896/133‑1 (MQ=255)
                                                                                    ataAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCg                                                                 >  4:249844/1‑149 (MQ=255)
                                                                                      aaaGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCAAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCt                                                               >  1:8298/1‑149 (MQ=255)
                                                                                      aaaGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCAACCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCt                                                               >  6:369262/1‑149 (MQ=255)
                                                                                      aaaGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCAACCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCt                                                               >  6:375734/1‑149 (MQ=255)
                                                                                          acacaACGGTTTCGTAAAATGTTCTCACATTACATGATTATGAGGCAACGCCATGAATCCACTTTTTCAAACCGCTTTTGCACAACTTGCGGATAACTTTCAATCTGCACTTGAACCTATTCTTGCAGACAAGTACTTCCACGCttttt                                                           <  1:79862/148‑3 (MQ=255)
                                                                                             caACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGa                                                        <  6:70912/149‑1 (MQ=255)
                                                                                              aaCGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCAACCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGAc                                                       <  5:375734/149‑1 (MQ=255)
                                                                                              aaCGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCAACCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGAc                                                       <  5:369262/149‑1 (MQ=255)
                                                                                                cGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCg                                                     <  2:8298/149‑1 (MQ=255)
                                                                                                 gggTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCgg                                                    <  3:80158/149‑1 (MQ=255)
                                                                                                   gTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCggg                                                   >  1:209677/1‑148 (MQ=255)
                                                                                                        tAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGtt                                                                                                                                               <  2:191123/51‑1 (MQ=255)
                                                                                                        tAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGtt                                                                                                                                               >  1:191123/1‑51 (MQ=255)
                                                                                                        tAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGcc                                                                                                                              <  2:400175/68‑1 (MQ=255)
                                                                                                        tAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGcc                                                                                                                              >  1:400175/1‑68 (MQ=255)
                                                                                                        tAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTg                                                                                                            >  8:302199/1‑86 (MQ=255)
                                                                                                        tAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTg                                                                                                            <  7:302199/86‑1 (MQ=255)
                                                                                                        tAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGc                                               >  7:38416/1‑147 (MQ=255)
                                                                                                         aaaCTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACtt                                                                                                             >  2:122692/1‑84 (MQ=255)
                                                                                                         aaaCTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACtt                                                                                                             <  1:122692/84‑1 (MQ=255)
                                                                                                         aaaCTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGGTGACCGGGGAGCAAg                                            >  1:213122/1‑149 (MQ=255)
                                                                                                                 aTCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCg                                    >  6:93781/1‑149 (MQ=255)
                                                                                                                    ccATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCt                                  <  5:391933/148‑1 (MQ=255)
                                                                                                                     cATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCAAACTTGCGGATAACTTGCAATCTGCAATGGAACCTATTCTGGCAGACAAGTACTTCAACGCTTTGTTGACCGTGGAGCAAGTCTCATCGCTg                                 >  3:175820/1‑148 (MQ=255)
                                                                                                                       ttACAAGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGCAACCTATTCTGGCAGACAACTACATCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTg                                 >  7:121559/1‑146 (MQ=255)
                                                                                                                          cATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAgc                           >  2:429446/1‑149 (MQ=255)
                                                                                                                           aTGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAgcg                          >  7:260054/1‑149 (MQ=255)
                                                                                                                            tGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAgcg                          >  2:212644/1‑148 (MQ=255)
                                                                                                                            tGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAgcg                          >  3:467/1‑148 (MQ=255)
                                                                                                                             gATTATGAGGCAACGCCATGCATCCACGTTTTCAACCCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAgcgc                         >  4:248620/1‑148 (MQ=255)
                                                                                                                             gATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAgcgc                         >  5:412230/1‑148 (MQ=255)
                                                                                                                              aTTATGAGGCAACGCCATGCATCCACGTTTTCACACCGCTTTTGCCCAACTTGCCGATAACTTGCAATCTGCACTGGAACCTCTTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCACATCGCTGAAGAGCGCaa                       >  2:316778/1‑149 (MQ=255)
                                                                                                                                tATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACg                     >  1:125308/1‑149 (MQ=255)
                                                                                                                                   gAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACgggg                  >  6:123621/1‑149 (MQ=255)
                                                                                                                                            ccATGCATCCACGTGTTCACACGGCGTTTGCGCAACTTGCGGATACCTTGCAATCTGCACTGGAACCTATTCGGGCAGACAAGTACTTCCGGGCTTTGTTGCCCGGGGAGCACGTCTCGTCGCTGAAGAGCGCACCGGGGCGGGACg           >  8:109784/1‑147 (MQ=255)
                                                                                                                                             cATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACg           >  5:46743/1‑146 (MQ=255)
                                                                                                                                              aTGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAATTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGa       >  5:108561/1‑149 (MQ=255)
                                                                                                                                              aTGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGa       >  4:203360/1‑149 (MQ=255)
                                                                                                                                              aTGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGa       >  2:230011/1‑149 (MQ=255)
                                                                                                                                                 cATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCt                <  6:403242/137‑1 (MQ=255)
                                                                                                                                                 cATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCt                >  5:403242/1‑137 (MQ=255)
                                                                                                                                                   tCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACgcgc  >  2:298100/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGC  >  NC_000913/2231702‑2231997

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑