Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,145 C→G 12.4% intergenic (+128/‑70) lldD → / → trmL L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1450CG12.4% 2600.8 / 130.3 868intergenic (+128/‑70)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (553/199);  new base G (2/109);  total (555/313)
Fisher's exact test for biased strand distribution p-value = 0.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 875 total aligned reads displayed.

GCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCT  >  NC_000913/3780997‑3781293
                                                                                                                                                     |                                                                                                                                                     
gCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGcc                                                                                                                                                       >  2:898534/1‑149 (MQ=255)
  gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                                                                                                                     <  1:614983/149‑1 (MQ=255)
  gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                                                                                                                     <  4:632862/149‑1 (MQ=255)
  gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                                                                                                                     >  3:1396324/1‑149 (MQ=255)
   aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                                                                                    >  4:1754981/1‑149 (MQ=255)
     aGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                                                                                                                  <  1:363163/149‑1 (MQ=255)
     aGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                                                                                                                  <  2:1402083/149‑1 (MQ=255)
     aGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                                                                                                                  <  2:83355/149‑1 (MQ=255)
      gggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTg                                                                                                                                                 >  3:463039/1‑149 (MQ=255)
        gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                                                                                                                               <  4:899544/149‑1 (MQ=255)
        gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                                                                                                                               <  2:589872/149‑1 (MQ=255)
        gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                                                                                                                               <  3:2057377/149‑1 (MQ=255)
        gAATGCAGCATGGTCGTTTGCGCCCCCTCACCCTAACCCGCTCCCTCAGGGAGAGGGGACCGTTTGGCGCTGTATGTACTCCCTCACTCTGAAACGGCACCGCACTCTTTTTTTCTCCCTCGCCCCTTCGGGGAGGGGGCCGGGGGGAg                                                                                                                                               >  2:1955739/1‑149 (MQ=255)
          aTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                                                                                                                             <  2:659929/149‑1 (MQ=255)
          aTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                                                                                                                             <  2:756343/149‑1 (MQ=255)
          aTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACCCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                                                                                                                             >  3:520862/1‑149 (MQ=255)
           tGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgggg                                                                                                                                            >  2:377790/1‑149 (MQ=255)
           tGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgggg                                                                                                                                            <  2:1016224/149‑1 (MQ=255)
                       gTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                                                                               >  1:2466198/1‑149 (MQ=255)
                       gTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                                                                               >  1:912832/1‑149 (MQ=255)
                        tttGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGccagatc                                                                                                                              >  3:1010014/1‑144 (MQ=255)
                        tttGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACt                                                                                                                              >  3:633313/1‑149 (MQ=255)
                             gcCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCa                                                                                                                         >  4:270527/1‑149 (MQ=255)
                                 ccccccccTACCCCCCCCCCCCGGGGGGGGGGGCCCGTCCGGCGTTTTAGGTTCCCCCCCCCGCTAAACCGCCCCCCCCTCTTTTTTTTCCCCCCCGCCCCGCGGGGGGGGGGGCGGGGTGGGGGGAAAAAGGCCGCACTGCTCACCaa                                                                                                                     <  3:2461164/144‑1 (MQ=255)
                                  cacACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGATGGGAAAAGGCCGCACTGCTCACCCAt                                                                                                                    >  1:2101332/3‑149 (MQ=255)
                                  cTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAt                                                                                                                    >  3:1565503/1‑149 (MQ=255)
                                   tCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCCCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtt                                                                                                                   >  4:222914/1‑149 (MQ=255)
                                     aCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                                                                                                                  >  1:817138/1‑117 (MQ=255)
                                      cccTACCCCTCTCCCTGAGGGAGGGGGGACCTTTGGGTTTTCTGTTTCCCCCCAAATTCCAAAAGGGCCCCGCACGTTTTATTCCCTGCCCCCCCCTCGGGGGGGGGGGCGGGGGTGAGGGGTAAAGGCCGCACTGCTCACCACTGTTc                                                                                                                <  4:1576733/149‑1 (MQ=255)
                                       ccTAACCCTCTGCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                                                               >  1:1315687/1‑149 (MQ=255)
                                       ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgggg                                                                                                                                                   >  1:2380279/1‑114 (MQ=255)
                                        caaaCCCCCCCCCTGAGGGAGGGGGGCCCTTCCGGCGCTTTATGTCCCCCCCCTCTCAGAACCGCCCCCCCCCTTTTTTTTCCCCCCTCCCCCCTCCGGGAAGGGGGCGGGGGAGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                                                                              <  3:1861811/147‑1 (MQ=255)
                                             ccTCTCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                                                                         >  3:410704/1‑149 (MQ=255)
                                             ccTCTCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                                                                         >  1:1601927/1‑149 (MQ=255)
                                             ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCTCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                                                                         >  4:1472048/1‑149 (MQ=255)
                                               tctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                                                                       >  1:26530/1‑149 (MQ=255)
                                                   ccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                                                                   >  3:1673327/1‑149 (MQ=255)
                                                   ccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                                                                   >  1:1006455/1‑149 (MQ=255)
                                                   ccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                                                                   >  2:1237046/1‑149 (MQ=255)
                                                    cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                                                              <  2:1593560/106‑1 (MQ=255)
                                                         ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGGTATCCTGCCCCCGCAc                                                                                             >  3:2119161/1‑149 (MQ=255)
                                                         ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                                                                             >  4:1391979/1‑149 (MQ=255)
                                                         ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                                                                             >  4:208158/1‑149 (MQ=255)
                                                             gggggACCGTGCGGCGTTGTATGTCCCCCCTCCCTCTGAAACGACCCCGCACTCTTTTTTTCCCCCTCCCCCCTCGGGGGAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg                                                                                         <  3:1257373/148‑1 (MQ=255)
                                                               gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                                                                       >  1:174386/1‑149 (MQ=255)
                                                                ggACCGTTCGGCGCTGTATTTCCCCCCCCTCTCTGAAACGACACCGCTCTCTTTTTTTCCCCCCCCCCCCCCGGGGGAGAGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                                                                                      <  2:176957/149‑1 (MQ=255)
                                                                ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                                                                                      >  3:594601/1‑149 (MQ=255)
                                                                 gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGCGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCACTTTTCTGCTATCCTGCCCCCGCACTAAggggg                                                                                     >  1:473899/1‑149 (MQ=255)
                                                                  aCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGTGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                                                                    >  3:1293166/1‑149 (MQ=255)
                                                                   cttttCGGCGTTGTTTGTCCCCCCCCCCTCTAAAACGCCACCCCTCTTTTTTTCCCCCCCCCCCCCCTCGGGGGGGGGGCGGGGGGTGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGcg                                                                                   <  1:111000/146‑2 (MQ=255)
                                                                     gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGt                                                                                 >  1:1412253/1‑149 (MQ=255)
                                                                      ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCACTTTTCTGCTATCCTGCCCCCGCACTCAGGGGGCAGTa                                                                                >  1:1971449/1‑149 (MQ=255)
                                                                       tCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                                                               >  3:1945858/1‑149 (MQ=255)
                                                                            gttgtatgtcCCCCCCCCCTCTAAAACGCCCCCCCACTTTTTTTTTCCCCCCCCCCCCCCCGGGGAGGGGGCGGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                                                                          <  2:283416/147‑1 (MQ=255)
                                                                            gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                                                                          >  3:1506304/1‑149 (MQ=255)
                                                                             cTGTATGTCCCCCCCCCCTTTAAACCGACCCCGCCTTTTTTTTTTCCCCCCCGCCCCTCCGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                                                                         <  4:455676/149‑1 (MQ=255)
                                                                              ggtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGGGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAAc                                                                        >  4:776388/2‑149 (MQ=255)
                                                                               gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa                                                                       >  3:862379/1‑149 (MQ=255)
                                                                                   gcaCTCCCTCACTCTAAAACTACCCCGCACTTTTTTTTTCTCCCTCGCCCCTGCGGGGAGGGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                                                                   <  1:1519526/147‑1 (MQ=255)
                                                                                        cncTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt                                                              >  3:326169/1‑149 (MQ=255)
                                                                                         ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa                                                             >  3:938975/1‑149 (MQ=255)
                                                                                             aCTCTGAAACGACACCGCGCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGaa                                                         >  3:2296624/1‑149 (MQ=255)
                                                                                                  gAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGa                                                    >  3:79967/1‑149 (MQ=255)
                                                                                                    aaCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGGGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGaaa                                                  <  3:1183471/149‑1 (MQ=255)
                                                                                                    aaCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGaaa                                                  >  2:322106/1‑149 (MQ=255)
                                                                                                       gACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                                                                                                                     >  2:185080/1‑48 (MQ=255)
                                                                                                       gACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                                                                                                                     <  1:185080/48‑1 (MQ=255)
                                                                                                       gACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTc                                               >  1:2174608/1‑149 (MQ=255)
                                                                                                          acCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGc                                            >  2:2227400/1‑149 (MQ=255)
                                                                                                             gCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCCATTTTCTGCTATCCTGCGCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCaa                                         >  3:730389/1‑149 (MQ=255)
                                                                                                              cACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCaaa                                        >  1:125516/1‑149 (MQ=255)
                                                                                                                   tttttttcCCCCCCCCCCCCTCCGGGGAGGGCGCCGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACt                                    <  2:2389592/148‑1 (MQ=255)
                                                                                                                    tttttttCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTTCCCCCGCCCTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTg                                   >  4:1109004/1‑149 (MQ=255)
                                                                                                                    tttttttCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTg                                   >  3:1717480/1‑149 (MQ=255)
                                                                                                                    tttttttCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTg                                   >  1:2273363/1‑149 (MQ=255)
                                                                                                                    tttttttCCCCCCCCCCCCCCCGGGGGGGGGGGCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACt                                    <  1:32006/149‑1 (MQ=255)
                                                                                                                      ttttcCCCCCCCCCCCCCCCGGGGGGGGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGc                                 <  2:2099288/149‑1 (MQ=255)
                                                                                                                        cccggcccTCTCCCCGGAGGGGCGAGGGCCGGGGAGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCaa                               >  4:1041643/6‑149 (MQ=255)
                                                                                                                          tctcCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAca                             >  4:1442523/1‑149 (MQ=255)
                                                                                                                            cccccCGCCCCCCCGGGGGGGGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAcatc                           <  4:1294769/148‑1 (MQ=255)
                                                                                                                             ccccccccccccGGGGGGGAGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAcatca                          <  4:1234953/143‑1 (MQ=255)
                                                                                                                              tccccccccTCGGGGGAGGGGCCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAACTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAcatcat                         <  2:1109527/144‑1 (MQ=255)
                                                                                                                              ccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACGTCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAcatcat                         >  1:742359/1‑149 (MQ=255)
                                                                                                                              ccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAcatcat                         >  4:847856/1‑149 (MQ=255)
                                                                                                                               cTCGCGCCTCCGGGGAGCGGGCCGGGGTGATGGGAAAAGGCCGAACTGCTCACCAATTTTCTGCTATTCTGCACCCGCCCTAAGATGGCATTATGCTAAACATCGTACTTTATGAACTAGAAATTCCTACAAAAACGGGGAAGatcatc                        >  4:1680440/1‑149 (MQ=255)
                                                                                                                               cTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg                                                                                         <  4:396734/84‑1 (MQ=255)
                                                                                                                                tCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATcc                       >  3:2236732/1‑149 (MQ=255)
                                                                                                                                tCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATcc                       >  2:1338541/1‑149 (MQ=255)
                                                                                                                                tCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATcc                       >  4:2428933/1‑149 (MQ=255)
                                                                                                                                    cccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                                                                            <  1:2044671/92‑1 (MQ=255)
                                                                                                                                    cccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCt                   >  2:2046609/1‑149 (MQ=255)
                                                                                                                                     ccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCTTAGTGCGGGGGCAgg                                                                                 <  3:107604/86‑2 (MQ=38)
                                                                                                                                      cTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTc                                               >  2:2285601/1‑119 (MQ=255)
                                                                                                                                      cTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCttt                 >  4:1879758/1‑149 (MQ=255)
                                                                                                                                      cTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCttt                 >  2:2490539/1‑149 (MQ=255)
                                                                                                                                       tCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTg                >  4:1986361/1‑149 (MQ=255)
                                                                                                                                        cgggggAGGGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTgc               <  4:52703/147‑1 (MQ=255)
                                                                                                                                             gggggggCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTaa          <  2:481854/145‑1 (MQ=255)
                                                                                                                                                agGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCACACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATAc       <  1:54985/149‑1 (MQ=255)
                                                                                                                                                    ccGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGc   <  2:1490987/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCT  >  NC_000913/3780997‑3781293

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: