Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,088,091 G→C 6.1% intergenic (+34/+16) yiiG → / ← frvR DUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,088,0910GC 6.1% 573.3 / 21.7 159intergenic (+34/+16)yiiG/frvRDUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain
Reads supporting (aligned to +/- strand):  ref base G (88/61);  new base C (5/4);  total (94/65)
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.42e-01

Only 100 of 168 total aligned reads displayed.

ACTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGC  >  NC_000913/4087942‑4088240
                                                                                                                                                     |                                                                                                                                                     
aCTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcgg                                                                                                                                                       <  1:1071774/150‑1 (MQ=255)
  tGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggcg                                                                                                                                                     >  1:68840/1‑150 (MQ=255)
  tGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggcg                                                                                                                                                     >  1:2603822/1‑150 (MQ=255)
     gATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGTAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCAt                                                                                                                                                  >  1:2244228/1‑150 (MQ=255)
     gATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCAt                                                                                                                                                  >  1:2657961/1‑150 (MQ=255)
      aTAAAGAGCAACTTCAGGACGCAAACTCAATCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATc                                                                                                                                                 <  2:392986/150‑1 (MQ=255)
      aTAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAACCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATc                                                                                                                                                 <  2:487283/150‑1 (MQ=255)
      aTAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATc                                                                                                                                                 <  2:392987/150‑1 (MQ=255)
      aTAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATc                                                                                                                                                 <  2:2847278/150‑1 (MQ=255)
       tAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATcc                                                                                                                                                <  1:3247361/150‑1 (MQ=255)
       tAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATcc                                                                                                                                                >  1:2499594/1‑150 (MQ=255)
         aaGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGa                                                                                                                                              <  2:1468080/150‑1 (MQ=255)
           gagCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGaca                                                                                                                                            <  2:899141/150‑1 (MQ=255)
            agCAACTTCAGGACGCAAACTCAAGCTGGAAGGTCGAAGACACTATCCCGAGAGCATTACGCGAGAACAACGAAATGGGTGATGACAATAAAAGCCTGCGTTAAACGATTTTGATGAACAATGCAGAGTCGGATGCGGCGCATCCCacac                                                                                                                                           >  2:810379/1‑150 (MQ=255)
                  ttCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCa                                                                                                                                     >  2:2809932/1‑150 (MQ=255)
                        aCGCAATCTCAAGCTGGATGGTCGAAGTCTCATTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccccaa                                                                                                                               >  2:2495279/1‑148 (MQ=255)
                        aCGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCtt                                                                                                                               >  1:2406826/1‑150 (MQ=255)
                        aCGAAAACTCAAGCTGGATGGTCTAAGACTCTTTCCCGAGAGCATTACGCGATTAAAACGAAATGGTTGATGACTATAATCGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTGGGATGCGGCGCATCCGACACAGTTCACCCCtt                                                                                                                               <  1:2125109/150‑1 (MQ=255)
                          gCAAACTCAAGCTGTATGGTCTAAGACTCTTTCCCGAGTGAATTACGCGAGTACAACGAATTGGTTGATGACTATAATCGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtgt                                                                                                                             <  1:2073444/150‑1 (MQ=255)
                          gCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtgt                                                                                                                             >  2:3089839/1‑150 (MQ=255)
                          gCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtgt                                                                                                                             >  1:3209246/1‑150 (MQ=255)
                            aaaCTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGa                                                                                                                           >  2:625619/1‑150 (MQ=255)
                            aaaCTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGa                                                                                                                           >  1:2190723/1‑150 (MQ=255)
                            aaaCTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGa                                                                                                                           >  1:581978/1‑150 (MQ=255)
                             aaCTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                                                                          <  2:2557054/150‑1 (MQ=255)
                             aaCTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCAGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                                                                          <  2:1777026/150‑1 (MQ=255)
                              aCTCATGCTGGATGGTCGAAGCCTCTTTCCCGAGTGCTTTACGCGTGTACAACGAAATGTTTGATGAATATAATAGCCTGATTTAAACGTTTTTGATTATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAg                                                                                                                         <  1:1817632/150‑1 (MQ=255)
                                  atggtggttggtCTACGACTCTTCCCCGAGAGCATTCCGCGATTACAACGAATTGGTTGTTGACTATAATCGCCTGCTTTAAAGTTTTTTTATGATAATGGCCGTGTCGGTTGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAAc                                                                                                                     <  1:1409077/146‑1 (MQ=255)
                                  aaGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAAc                                                                                                                     <  2:1359170/150‑1 (MQ=255)
                                    gCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTc                                                                                                                   >  2:2089163/1‑150 (MQ=255)
                                          ttgTCTAATACTCTTTCCCGAGAGCATTACGCGAGTACATATAATTGGTTGTTGAATATAATAGCCTGCGTTAAACGTTTTTTATTATAAATGCAGTGTCGTATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaaa                                                                                                             <  1:117237/148‑1 (MQ=255)
                                          tGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaaa                                                                                                             >  2:392482/1‑150 (MQ=255)
                                           ggTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAAt                                                                                                            <  2:1899344/150‑1 (MQ=255)
                                              cGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGAAAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                                                                         <  1:1376613/150‑1 (MQ=255)
                                                     tctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATc                                                                                                  <  1:2021097/150‑1 (MQ=255)
                                                      ctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACTTTTTTGATGATAAATGCTGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATcc                                                                                                 <  1:932443/150‑1 (MQ=255)
                                                       tttGCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTATAGGTTTTTGATGTTAAATGCAGTGTCGGTTGCGGCGCATCCGACACAGTTCACCCCTTGTGTAGTTACTCCATAAATGCCTCATATCCt                                                                                                >  2:1064405/1‑150 (MQ=255)
                                                       tttCCCGAGAGCATTACGCGAGTACAACGAATTGGTTGATGACTATAATAGCCTGAGTTAAACGTTTTTGATGATAAATGCAGTGTCGTATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCt                                                                                                <  2:619391/150‑1 (MQ=255)
                                                       tttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATTCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCt                                                                                                >  1:1699741/1‑150 (MQ=255)
                                                       tttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCCTCCGACACAGTTCACCCCATGTGAAGTAACGCCATAAATGCCTCATATccg                                                                                                >  1:676007/1‑149 (MQ=255)
                                                       tttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCt                                                                                                >  2:2315188/1‑150 (MQ=255)
                                                       tttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCt                                                                                                >  1:1170368/1‑150 (MQ=255)
                                                       tttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACACCATAAATGCCTCATATCCt                                                                                                >  1:3027288/1‑150 (MQ=255)
                                                       tttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGACGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATACt                                                                                                >  2:2315202/1‑150 (MQ=255)
                                                        ttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                                                                               >  2:2900907/1‑150 (MQ=255)
                                                        ttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGAGGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                                                                               <  1:2829050/150‑1 (MQ=255)
                                                           ccGAGAGCATTACGCGAGTACAAAGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTa                                                                                            >  1:2740000/1‑150 (MQ=255)
                                                                agCATTACGCGAGTACAACGACATGGTTGATGACTATAATAGCCTGCTTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCg                                                                                       <  1:3015510/150‑1 (MQ=255)
                                                                 gCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg                                                                                      <  1:486695/150‑1 (MQ=255)
                                                                 gCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg                                                                                      <  2:961152/150‑1 (MQ=255)
                                                                          gAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtct                                                                             >  2:1787074/1‑150 (MQ=255)
                                                                          gAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtca                                                                             >  2:1787282/1‑149 (MQ=255)
                                                                           aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtctc                                                                            <  2:1258748/150‑1 (MQ=255)
                                                                             tACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg                                                                          >  2:3237634/1‑150 (MQ=255)
                                                                                aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGTAGTAACTCCATAAATTGCTCATATGCTGTTTACCCGGCAATAATATccgccg                                                                       >  1:1885471/1‑150 (MQ=255)
                                                                                aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCAACCGACACAGTTAACCCCGTGTGTAGTAACTCCATAAATGCCTCAAATCCGGTTAAACCGGCAATAATATccgccg                                                                       >  1:1567286/1‑150 (MQ=255)
                                                                                 aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCgg                                                                      >  2:2132557/1‑150 (MQ=255)
                                                                                  cGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa                                                                     <  1:1349369/150‑1 (MQ=255)
                                                                                   gAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGAt                                                                    >  2:2520439/1‑150 (MQ=255)
                                                                                   gAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCATGTGAAGTAAATCCATAAAATCCTCATATCCTGTTAACCCGGCAAAAATCTCCGCCGGAt                                                                    >  1:2440447/1‑150 (MQ=255)
                                                                                    aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATACCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATg                                                                   >  2:157010/1‑150 (MQ=255)
                                                                                    aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATACCTCAAATCCTGTTAACCCGGCAATAATCTCCGCCGGATg                                                                   >  2:1688362/1‑150 (MQ=255)
                                                                                       tGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTg                                                                <  2:3149272/150‑1 (MQ=255)
                                                                                         gTTGATGACTATAATAGCCTGCGTTAACCTTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac                                                              <  1:3064063/150‑1 (MQ=255)
                                                                                         gTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac                                                              <  1:3064047/150‑1 (MQ=255)
                                                                                           tgatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGacac                                                            >  2:2448752/1‑150 (MQ=255)
                                                                                            gatgaCTATAATAGCCTGCGTTAATCGTTTTTGATGATAAATGCAGTGTCGGTTGCGGCGCATCCGTCACAGTTCACCCCTTGTGTAGTAACTCCATAAATGCGTCATATCGTGTTTACCCGGCCATAATGTCCGCCGGATGCTGAcaca                                                           >  1:1896478/1‑150 (MQ=255)
                                                                                             atgaCTATAATGGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACaa                                                          <  2:196500/150‑1 (MQ=255)
                                                                                             atgaCTATAATGGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACaa                                                          <  2:2500928/150‑1 (MQ=255)
                                                                                               gaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAata                                                        >  2:2417541/1‑150 (MQ=255)
                                                                                                    taataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCtt                                                   <  1:918183/150‑1 (MQ=255)
                                                                                                    taataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCACTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCtt                                                   >  2:1082568/1‑150 (MQ=255)
                                                                                                     aataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCGTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTg                                                  >  2:3199320/1‑150 (MQ=255)
                                                                                                          ccTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCCTAGCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGa                                             >  2:2671642/1‑150 (MQ=255)
                                                                                                           cTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGAt                                            <  2:2365758/150‑1 (MQ=255)
                                                                                                           cTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGAt                                            >  2:2679488/1‑150 (MQ=255)
                                                                                                            tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATa                                           <  1:1288972/150‑1 (MQ=255)
                                                                                                            tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATa                                           <  1:657287/150‑1 (MQ=255)
                                                                                                              cGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATACTTGCCAGATAGc                                         <  2:1607020/150‑1 (MQ=255)
                                                                                                               gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCt                                        >  2:814141/1‑150 (MQ=255)
                                                                                                                         tttGATGATAATCGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATAttttt                              <  1:886864/150‑1 (MQ=255)
                                                                                                                         tttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATAttttt                              >  2:2253430/1‑150 (MQ=255)
                                                                                                                         tttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCACATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATAttttt                              >  2:2253427/1‑150 (MQ=255)
                                                                                                                               gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCa                        >  2:1093645/1‑150 (MQ=255)
                                                                                                                               gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAAAATATTTAGCtat                        <  1:2404991/150‑3 (MQ=255)
                                                                                                                                atAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCAt                       >  2:2422681/1‑150 (MQ=255)
                                                                                                                                            tCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcgg           <  2:80711/150‑1 (MQ=255)
                                                                                                                                            tCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcgg           <  1:1018608/150‑1 (MQ=255)
                                                                                                                                            tCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcgg           >  2:1475997/1‑150 (MQ=255)
                                                                                                                                            tCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcgg           >  2:1966793/1‑150 (MQ=255)
                                                                                                                                            tCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcgg           >  2:321792/1‑150 (MQ=255)
                                                                                                                                            tCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcgg           >  2:3063705/1‑150 (MQ=255)
                                                                                                                                            tCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcgg           >  2:2181424/1‑150 (MQ=255)
                                                                                                                                               gATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcggcgg        <  2:2265239/150‑1 (MQ=255)
                                                                                                                                                 tGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcggcggcg      >  1:369675/1‑150 (MQ=255)
                                                                                                                                                   cggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATcc                                                                                                 <  2:2067853/57‑1 (MQ=255)
                                                                                                                                                   cggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATcc                                                                                                 >  1:2067853/1‑57 (MQ=255)
                                                                                                                                                   cggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcggcggcggt    >  1:2156615/1‑149 (MQ=255)
                                                                                                                                                     gcgcATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGgcgc  >  1:1734144/1‑150 (MQ=255)
                                                                                                                                                     gcgcATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGgcgc  >  1:3090806/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGC  >  NC_000913/4087942‑4088240

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: