Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 53,272 T→G 5.6% M49L (ATG→CTG)  pdxA ← 4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091353,2720TG 5.6% 337.5 / 8.3 119M49L (ATG→CTG) pdxA4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base T (51/61);  new base G (5/2);  total (56/63)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

Only 100 of 121 total aligned reads displayed.

CACCACATAATGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTT  >  NC_000913/53123‑53421
                                                                                                                                                     |                                                                                                                                                     
ccccacATAATGCCCATTTTCCACCGCTACCTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCAt                                                                                                                                                       <  2:29566/148‑1 (MQ=255)
    acaTAATGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCt                                                                                                                                                   >  1:1778271/1‑150 (MQ=255)
      aTAATGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGc                                                                                                                                                 <  1:931570/150‑1 (MQ=255)
        aaTGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGccc                                                                                                                                               <  1:1467250/150‑1 (MQ=255)
            cccATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGTAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGtt                                                                                                                                           <  2:1261856/150‑1 (MQ=255)
             ccATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTg                                                                                                                                          <  2:448528/150‑1 (MQ=255)
             ccATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTg                                                                                                                                          >  1:1245987/1‑150 (MQ=255)
             ccATTTTCAACCGCAAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTg                                                                                                                                          <  2:238224/150‑1 (MQ=255)
              cATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTgg                                                                                                                                         >  2:901271/1‑150 (MQ=255)
               aTTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGt                                                                                                                                        >  2:216231/1‑150 (MQ=255)
                ttttCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGGGCAGCGATGCTCGGtttgcc                                                                                                                                       <  2:565162/150‑5 (MQ=255)
                ttttCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTGTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGTGGGTGAGCGGCAACCCGAGCATCGCTGCCCGGTGGGTa                                                                                                                                       >  2:340091/1‑150 (MQ=255)
                tgcacaaccgcaaactgcgggcacattaacgctacttcctgtcgcgctacgtgcacctgtcactgcggggcagttagcggttgTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTa                                                                                                                                       >  1:43364‑M2/84‑150 (MQ=255)
                    CAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGga                                                                                                                                   <  2:1558644‑M1/150‑1 (MQ=255)
                     AACCGCTACCTGCCCCGCAGAGACAGGTGCACGTAGCGCGACAGGAACTAGCGTAAACGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGCTGGGGGAATAAGGGCGGCGGGGGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGgag                                                                                                                                  >  2:1195151‑M1/1‑150 (MQ=255)
                      ACCGCTAACTGCCCCGCAGTGACAGGTGCACGTCGCGCGCCAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTGGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGCGCGGCAAACCGCGCCTCGCTGCCCGGTTGGTAAGgaga                                                                                                                                 <  1:1509674‑M1/150‑1 (MQ=255)
                      ACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGgaga                                                                                                                                 <  1:17353‑M1/150‑1 (MQ=255)
                        CGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGt                                                                                                                               <  1:418356‑M1/150‑1 (MQ=255)
                          CTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGGGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTgg                                                                                                                             >  2:1095220‑M1/1‑150 (MQ=255)
                            aaCTGCCCCGCAGTGACAGGTGCACGTAGCGCGACCGGAAGTAGCGTTAATGTGCCCGCAGTTCGCGGCTGTGCAGGGGAGGTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCa                                                                                                                           <  2:910700/150‑1 (MQ=255)
                              cTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATc                                                                                                                         >  1:1187180/1‑150 (MQ=255)
                                gCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCgg                                                                                                                       >  1:1207034/1‑150 (MQ=255)
                                  cccGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGca                                                                                                                     >  1:667310/1‑150 (MQ=255)
                                    cGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGGTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGcaca                                                                                                                   <  1:441979/150‑1 (MQ=255)
                                        gTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCacaacca                                                                                                               >  2:51889/1‑150 (MQ=255)
                                          gACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAaacaac                                                                                                             >  2:321773/1‑150 (MQ=255)
                                           aCAGGTGCACGTAGCGCGACAGGACGTAGCGGTAATGTGGCCGCAGTTTGCGGGTGTGCAGGGGAGGTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAAcc                                                                                                            <  2:237177/150‑1 (MQ=255)
                                           aCAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAAcc                                                                                                            >  1:1478563/1‑150 (MQ=255)
                                              ggTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGt                                                                                                         >  1:802140/1‑150 (MQ=255)
                                               gTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGtt                                                                                                        <  2:1384215/150‑1 (MQ=255)
                                               gTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGCCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGtt                                                                                                        >  1:1210360/1‑150 (MQ=255)
                                                tGCACGTAGCGCGCCAGGCAGTAGCGTTAATGTGCCCGCCGTTTGCGGGTGTGCAGGCGAGTTGGGGGAATAAGGGCGGAGGGGGAGCGGCAAACCGAGCATCGCTGCCCGGGAGGTAAGGAGCGTGCCATCGGCACAAACACCCAGtgc                                                                                                       >  2:290298/1‑148 (MQ=255)
                                                     gTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGGAGTTTGCGGGTGTGCAGGGGAGGTGGGGGAATTAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCg                                                                                                  <  1:946862/150‑1 (MQ=255)
                                                     gTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCg                                                                                                  >  2:719450/1‑150 (MQ=255)
                                                      tAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCgg                                                                                                 >  2:842479/1‑150 (MQ=255)
                                                       aGCGCGACAGGAAGTAGCGTGAATGTGCCCGCAGTTGGCGGTTGTGCAGGGGAGTTGGGGGAATAATGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGCGGCATCGGCACAAACAACCAGTTCGACCGGc                                                                                                <  1:339134/150‑1 (MQ=255)
                                                        gcgcGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGcc                                                                                               <  1:1058552/150‑1 (MQ=255)
                                                         cgcgACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCa                                                                                              <  1:210375/150‑1 (MQ=255)
                                                         cgcgACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCa                                                                                              >  2:441751/1‑150 (MQ=255)
                                                           cgACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACt                                                                                            >  2:1640991/1‑150 (MQ=255)
                                                               aGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACg                                                                                        >  2:627518/1‑150 (MQ=255)
                                                                ggAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGCCCGGCCACTCACGc                                                                                       >  2:366722/1‑150 (MQ=255)
                                                                 gAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCt                                                                                      <  2:461947/150‑1 (MQ=255)
                                                                   aGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAACGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCtgt                                                                                    <  1:291796/150‑1 (MQ=255)
                                                                    ttatcGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGTGCGGATGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGTCATCGGCACAAACAACCAGTTCGACCTGCCACTCACGCtgtg                                                                                   <  1:513508/146‑1 (MQ=255)
                                                                    gTAGCGTTACTCTGCCCGCAGTTTGCGGTTGTGCAGGCGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCGGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCtgtg                                                                                   <  1:293284/150‑1 (MQ=255)
                                                                    gTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCtgtg                                                                                   <  1:547853/150‑1 (MQ=255)
                                                                    gTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCtgtg                                                                                   <  2:1778271/150‑1 (MQ=255)
                                                                    gTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCtgtg                                                                                   <  1:1781916/150‑1 (MQ=255)
                                                                         gTTAATGTGCCCGCTGTTGGCGGTTGTGCAGGGGAGGTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTTGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGc                                                                              <  1:176440/150‑1 (MQ=255)
                                                                           tactgtgCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGAAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTg                                                                            <  1:1630059/147‑1 (MQ=255)
                                                                           tAATGTGCCCGCAGGTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTg                                                                            <  1:1029312/150‑1 (MQ=255)
                                                                            aaTGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAAAAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGGTTCGGCACAAACAACCAGTTCGACCGGCCAATCACGCTGTGCAAGCTgg                                                                           >  1:66169/1‑150 (MQ=255)
                                                                               gtgCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACa                                                                        >  2:1778617/1‑150 (MQ=255)
                                                                                tgCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACaa                                                                       >  2:911011/1‑150 (MQ=255)
                                                                                   ccGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTa                                                                    <  1:1772058/150‑1 (MQ=255)
                                                                                     gCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGACACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAg                                                                  <  2:1267168/150‑1 (MQ=255)
                                                                                        gTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTcc                                                               >  1:1311586/1‑150 (MQ=255)
                                                                                          ttGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCCTCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCgg                                                             <  2:126755/150‑1 (MQ=255)
                                                                                          ttGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAACCAACAAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCgg                                                             >  1:273507/1‑150 (MQ=255)
                                                                                          ttGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCgg                                                             <  1:1808065/150‑1 (MQ=255)
                                                                                           tGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCggg                                                            >  1:1399867/1‑150 (MQ=255)
                                                                                            gcggttgTGCAGGGGAGTAGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGc                                                           <  1:1177725‑M2/143‑1 (MQ=255)
                                                                                             cggttgTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGcc                                                          <  2:418356‑M2/144‑1 (MQ=255)
                                                                                             cggttgTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGcc                                                          >  1:13598‑M2/7‑150 (MQ=255)
                                                                                                    gCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCcg                                                   >  1:478300/1‑150 (MQ=255)
                                                                                                      aGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAGCCAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCcggc                                                 >  2:841009/1‑150 (MQ=255)
                                                                                                       ggggAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTTGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCcggcg                                                <  1:1733415/150‑1 (MQ=255)
                                                                                                       ggggAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCcggcg                                                <  1:1499050/150‑1 (MQ=255)
                                                                                                       ggggAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGGGCAAGCTGGACAACTAAGTCCGGGCCAATCCcggcg                                                <  1:1332914/150‑1 (MQ=255)
                                                                                                          gcgTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGTCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGc                                             <  1:1124777/148‑1 (MQ=255)
                                                                                                                ggggAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCg                                       >  1:1439162/1‑150 (MQ=255)
                                                                                                                  ggAGTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACATCTAAGTCCGGGCCAATCCCGGCGAGCCCGCCggg                                     <  2:1403798/150‑1 (MQ=255)
                                                                                                                  ggAATAAGTGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTGCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCggg                                     <  2:857812/150‑1 (MQ=255)
                                                                                                                    aaTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAg                                   <  1:49407/150‑1 (MQ=255)
                                                                                                                     aTAAGGGCGGAGGGGGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGCCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCGGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCGGGatt                                  <  2:874792/150‑3 (MQ=255)
                                                                                                                      tAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCAGCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCGCGCCGGGAGTg                                 <  2:395243/150‑1 (MQ=255)
                                                                                                                       aaGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGa                                <  1:321773/150‑1 (MQ=255)
                                                                                                                         gggCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGGTGGTAAGGAGAGTGGCATCGGCACAAACAACCCGATCGACCGGCATCTCACGCTGGGCAAGCTGGACAACTAAGTCCGTGCCAATCCCGGCGGGCTCGCCGGGAGTGATc                              <  1:1211811/150‑1 (MQ=255)
                                                                                                                           gCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTACGGAGAGTCGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCGGAGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATcac                            <  1:298735/150‑1 (MQ=255)
                                                                                                                              gAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCacaac                         >  2:1876323/1‑150 (MQ=255)
                                                                                                                                    gAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACGCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTg                   <  2:1382789/150‑1 (MQ=255)
                                                                                                                                     aGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTgg                  <  2:667310/150‑1 (MQ=255)
                                                                                                                                      gCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTggg                 >  2:1683516/1‑150 (MQ=255)
                                                                                                                                       cGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGt                >  2:331366/1‑150 (MQ=255)
                                                                                                                                       cGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGt                >  2:1561144/1‑150 (MQ=255)
                                                                                                                                          ccaaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGtttt             >  1:1288983/3‑150 (MQ=255)
                                                                                                                                          ccaaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGGCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGtttt             >  1:69346/3‑150 (MQ=255)
                                                                                                                                          cAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGtttt             >  1:717668/1‑150 (MQ=255)
                                                                                                                                           caaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTa            >  2:1060763/2‑150 (MQ=255)
                                                                                                                                             aCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAc          >  2:1549513/1‑150 (MQ=255)
                                                                                                                                             aCCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAc          >  2:1713660/1‑150 (MQ=255)
                                                                                                                                             aCCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAc          >  2:847911/1‑150 (MQ=255)
                                                                                                                                              ccGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAcc         <  2:1245987/150‑1 (MQ=255)
                                                                                                                                              ccGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCGCGCCGGGAGTGATCACAACACGTTGGGTTTTAAcc         <  2:1447034/150‑1 (MQ=255)
                                                                                                                                               cGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCa        >  2:919923/1‑150 (MQ=255)
                                                                                                                                                   cATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTTACCATTAg    <  1:1640991/150‑1 (MQ=255)
                                                                                                                                                   cATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGATCGCCGGGAGTGATCACACCACGTTGGGGTTTAACCATTAg    >  2:1105780/1‑150 (MQ=255)
                                                                                                                                                   cATCGCTGCCCGGTTGGTAAGGAGAGTGGCAGCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAg    <  2:1861208/150‑1 (MQ=255)
                                                                                                                                                    aTCGCTGCCCGGTTTGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGt   <  1:1573872/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACCACATAATGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTT  >  NC_000913/53123‑53421

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