Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,505,785 A→T 6.6% I79N (ATT→AAT)  yhfS ← putative aminotransferase YhfS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,505,7850AT 6.6% 541.5 / 18.7 176I79N (ATT→AAT) yhfSputative aminotransferase YhfS
Reads supporting (aligned to +/- strand):  ref base A (83/80);  new base T (6/6);  total (90/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 182 total aligned reads displayed.

CCTGCTTCAGTGCCGACAGGTCATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGG  >  NC_000913/3505636‑3505934
                                                                                                                                                     |                                                                                                                                                     
ccTGCTTCAGTGCCGACAGGTCATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGaa                                                                                                                                                       >  2:2164950/1‑150 (MQ=255)
ccTGCTTCAGTGCCGACAGGTCATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGaa                                                                                                                                                       >  2:1055426/1‑150 (MQ=255)
    cTTCAGTGCCGACAGGTCATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATcgc                                                                                                                                                   >  2:2263001/1‑150 (MQ=255)
    cTTCAGTGCCGACAGGTCATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATcgc                                                                                                                                                   >  1:2048/1‑150 (MQ=255)
         gTGCCGACAGGTCATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGGCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAg                                                                                                                                              <  1:1111989/150‑1 (MQ=255)
          tGCCGACAGGTCATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGATTCGCGccggg                                                                                                                                             <  2:1792462/150‑4 (MQ=255)
           gCCGACAGGTCATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTc                                                                                                                                            >  1:1031045/1‑150 (MQ=255)
           gCCGACAGGTCATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTc                                                                                                                                            >  2:1011728/1‑150 (MQ=255)
           gCCGACAGGTCATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAACCGCGCCAGTc                                                                                                                                            <  1:2134227/150‑1 (MQ=255)
           gCCGACAGGTCATTGAAATCAACAGTAATAAGCGTCAGCCCCATATGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTc                                                                                                                                            <  2:796499/150‑1 (MQ=255)
                 aGGTCATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCcgcg                                                                                                                                      <  2:1267634/150‑1 (MQ=255)
                   gTCATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGcc                                                                                                                                    <  1:1765642/150‑1 (MQ=255)
                    tCATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGccc                                                                                                                                   >  2:978355/1‑150 (MQ=255)
                     cATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCt                                                                                                                                  <  2:154544/150‑1 (MQ=255)
                     cATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCt                                                                                                                                  >  1:380673/1‑150 (MQ=255)
                      aTTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTg                                                                                                                                 >  1:652893/1‑150 (MQ=255)
                       ttGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGc                                                                                                                                >  2:1947588/1‑150 (MQ=255)
                        tGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCa                                                                                                                               <  1:305487/150‑1 (MQ=255)
                           aaTCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCa                                                                                                                            <  2:2282576/150‑1 (MQ=255)
                             tCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAAGCCGTGTCGTCGGGAAAACAGCCCCGTCATGCACCAGAAGACGCTGCCCCGGTGTGAGCAAAGCCGCCACCCCCGCGCGACTCCCGCCAGACCCCGCGCCCTCCACCCgc                                                                                                                          >  2:1792122/1‑150 (MQ=255)
                                 cAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGcgccg                                                                                                                      <  1:823143/150‑1 (MQ=255)
                                  aGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGcgccgc                                                                                                                     >  2:2386596/1‑150 (MQ=255)
                                   gTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCa                                                                                                                    >  2:337922/1‑150 (MQ=255)
                                   gTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCa                                                                                                                    >  1:121641/1‑150 (MQ=255)
                                    taataaGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAg                                                                                                                   <  1:482631/150‑1 (MQ=255)
                                     aataaGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGc                                                                                                                  >  1:586530/1‑150 (MQ=255)
                                       taaGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCt                                                                                                                >  2:319777/1‑150 (MQ=255)
                                        aaGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCtg                                                                                                               >  2:1925612/1‑150 (MQ=255)
                                          gCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCtgtg                                                                                                             >  1:2095301/1‑150 (MQ=255)
                                          gCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCtgtg                                                                                                             <  2:275019/150‑1 (MQ=255)
                                          gCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCtgtg                                                                                                             <  2:503049/150‑1 (MQ=255)
                                            gTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCg                                                                                                           <  2:1099810/150‑1 (MQ=255)
                                            gTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCg                                                                                                           <  2:594793/150‑1 (MQ=255)
                                                  cccATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAAt                                                                                                     <  1:2074930/150‑1 (MQ=255)
                                                    cATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCCCCAGCGCCGCAGCCTGTGCGTGAAATGc                                                                                                   >  1:2084701/1‑150 (MQ=255)
                                                       gtGCTGAATAATATCCCGTGTGGTCGGGTAATCAGGCGCGGCATGCAGCTGAAGTCGCTGCCCCGGTTTGAGCAAAGGCGTGAGTCAGGTGCGAATCGGGCCAGTCGGCGCGGCCTGCACGAGCGCCGCAGGCTGTGCGTGATATGCAtg                                                                                                <  2:2439928/149‑2 (MQ=255)
                                                        tGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCCGCCTGTGCGTGAAATGCATca                                                                                               >  2:205221/1‑150 (MQ=255)
                                                         gCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATcag                                                                                              <  1:1180222/150‑1 (MQ=255)
                                                         gCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATcag                                                                                              <  1:1098281/150‑1 (MQ=255)
                                                           tCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATcagcc                                                                                            <  2:88781/150‑1 (MQ=255)
                                                            cAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATcagcca                                                                                           <  1:2488075/150‑1 (MQ=255)
                                                            cAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATcagcca                                                                                           >  1:1400143/1‑150 (MQ=255)
                                                            cAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATca                                                                                               >  1:134017/1‑146 (MQ=255)
                                                             aataatAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATcagccag                                                                                          <  1:339965/150‑1 (MQ=255)
                                                              ataataACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATcagccagc                                                                                         >  2:825716/1‑150 (MQ=255)
                                                               taataaCCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCa                                                                                        <  2:1076790/150‑1 (MQ=255)
                                                                 ataaCCCGTGTCGTCGGGTAAAGAGGCGCTTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCTCCAGCCCGGCGCGATTCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCAcc                                                                                      <  2:1326475/150‑1 (MQ=255)
                                                                 ataaCCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCAcc                                                                                      <  1:1453320/150‑1 (MQ=255)
                                                                     cccGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGtt                                                                                  >  1:1598892/1‑150 (MQ=255)
                                                                     cccGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATAAGCCAGCACCTGtt                                                                                  >  1:1542159/1‑150 (MQ=255)
                                                                        gtgtCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCca                                                                               <  2:1635551/150‑1 (MQ=255)
                                                                         tgtCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCcac                                                                              <  2:187714/150‑1 (MQ=255)
                                                                         tgtCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCcac                                                                              <  1:2021660/150‑1 (MQ=255)
                                                                            cgtcgGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGTACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACg                                                                           >  1:1674538/1‑150 (MQ=255)
                                                                             gtcgGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCGCCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGc                                                                          >  2:1179615/1‑150 (MQ=255)
                                                                              tcgGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCt                                                                         >  1:1823290/1‑150 (MQ=255)
                                                                               cgGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTg                                                                        <  2:380673/150‑1 (MQ=255)
                                                                               cgGGTAAACAGGCGCGTCATGCACCAGAACACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTg                                                                        >  2:780667/1‑150 (MQ=255)
                                                                                gggTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTgg                                                                       <  1:948166/150‑1 (MQ=255)
                                                                                gggTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTgg                                                                       >  2:2410776/1‑150 (MQ=255)
                                                                                gggTAAACAGGCGCGTCATGCACAAGAAGACGCTGCCCCGGATTGAGCAAAGCCGCCAGCCCGGCGCGAATCTCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATCCATCAGCCAGCACCTGTACCCCACGCTgg                                                                       >  2:1236585/1‑150 (MQ=255)
                                                                                  gTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGt                                                                     <  1:285022/150‑1 (MQ=255)
                                                                                     aaCAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAAc                                                                  <  2:1369853/150‑1 (MQ=255)
                                                                                     aaCAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAAc                                                                  >  2:174389/1‑150 (MQ=255)
                                                                                       cAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAAcgc                                                                <  1:2263001/150‑1 (MQ=255)
                                                                                         ggCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAAcgcgc                                                              >  1:677534/1‑150 (MQ=255)
                                                                                           cgcgTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCgg                                                            <  1:1011728/150‑1 (MQ=255)
                                                                                            gcgTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCTCGGt                                                           <  2:267596/150‑1 (MQ=255)
                                                                                            gcgTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGt                                                           <  2:2042669/150‑1 (MQ=255)
                                                                                             cttCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCGAGTCCCCGCTCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGtt                                                          <  2:721901/148‑1 (MQ=255)
                                                                                                    cACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAg                                                   <  1:978355/150‑1 (MQ=255)
                                                                                                    cACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCAGGCGCGAATCGCGCAAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAg                                                   >  2:2029908/1‑150 (MQ=255)
                                                                                                         gaagaCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTg                                              <  1:337922/150‑1 (MQ=255)
                                                                                                            gaCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCg                                           <  2:1821830/150‑1 (MQ=255)
                                                                                                            gaCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCg                                           <  1:2164950/150‑1 (MQ=255)
                                                                                                               gCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCCCCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTa                                        >  1:768573/1‑150 (MQ=255)
                                                                                                                   ccccGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGcc                                    <  1:872921/150‑1 (MQ=255)
                                                                                                                   ccccGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGcc                                    >  2:1843799/1‑150 (MQ=255)
                                                                                                                    cccGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCt                                   >  2:1445788/1‑150 (MQ=255)
                                                                                                                     ccGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTa                                  >  2:1674158/1‑150 (MQ=255)
                                                                                                                         tttGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATc                              <  2:1598892/150‑1 (MQ=255)
                                                                                                                          ttGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCa                             >  1:1209644/1‑150 (MQ=255)
                                                                                                                           tGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCAc                            >  1:1813665/1‑150 (MQ=255)
                                                                                                                              gCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCCCCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCCCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGc                         >  1:826071/1‑150 (MQ=255)
                                                                                                                               cAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGcc                        <  1:83359/150‑1 (MQ=255)
                                                                                                                                    ccgccAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTaa                   >  1:1949446/1‑150 (MQ=255)
                                                                                                                                    ccgccAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTACACACGCTGGGTAACGCGCGGTTGATTCACTCCTGGCGTTAAGCCTAAATCACCGCCAGTaa                   <  2:2066210/150‑1 (MQ=255)
                                                                                                                                        gggactggCGCGATTCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGaaa               >  1:128063/7‑150 (MQ=255)
                                                                                                                                         ggactggCGCGATTCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTATTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTaaa                                >  2:1527294/6‑132 (MQ=255)
                                                                                                                                          gactggCGCGATTCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACt             >  1:1522416/5‑150 (MQ=255)
                                                                                                                                           cccGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTc            <  1:1480347/150‑1 (MQ=255)
                                                                                                                                            ctggCGCGATTCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCg           >  2:1162953/3‑150 (MQ=255)
                                                                                                                                            ctggCGCGATTCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCACCAGTAAGAAACTCg           >  2:338368/3‑150 (MQ=255)
                                                                                                                                            ccGGCGCGAATCGCGCCAGTCCCCGCGCCCTGGACGAGCGCCGCAGCGTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTTAATCACCGCCAGTAAGAAACTCg           <  2:2006615/150‑1 (MQ=255)
                                                                                                                                            ccGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCg           <  2:96261/150‑1 (MQ=255)
                                                                                                                                                cgcgAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGc       >  1:2369600/1‑150 (MQ=255)
                                                                                                                                                  cgAATCGCGCCAGTCCCCGCGCCCTGCACCAGTGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACCCGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCg     >  1:2276908/1‑150 (MQ=255)
                                                                                                                                                   gAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCgg    >  2:1958804/1‑150 (MQ=255)
                                                                                                                                                   gAATCGCCCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATCCATCAGCCAGCACCTGTTCCACACGCTGCGTAACGCGCGGATCATTCAGTCCTGGCGATAAGCCTAAATCACCGCCAGTTAGCAACTCGCTGCCgg    >  2:1045987/1‑150 (MQ=255)
                                                                                                                                                     aTCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCgggg  <  2:1592774/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTGCTTCAGTGCCGACAGGTCATTGAAATCAACAGTAATAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGG  >  NC_000913/3505636‑3505934

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