Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 53,272 T→G 8.9% 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 8.9% 306.1 / 16.7 118M49L (ATG→CTG) pdxA4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base T (55/52);  new base G (8/3);  total (63/55)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

Only 100 of 122 total aligned reads displayed.

CACCACATAATGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTT  >  NC_000913/53123‑53421
                                                                                                                                                     |                                                                                                                                                     
caccacATAATGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCAt                                                                                                                                                       <  1:378537/150‑1 (MQ=255)
caccacATAATGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCAt                                                                                                                                                       >  2:502607/1‑150 (MQ=255)
caccacATAATGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCAt                                                                                                                                                       <  2:728458/150‑1 (MQ=255)
 accacATAATGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGCCAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGGTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGCTGAGCGGCAAACCGAGCATc                                                                                                                                                      >  1:285682/1‑150 (MQ=255)
   cacaTAATGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGc                                                                                                                                                    <  2:550259/150‑1 (MQ=255)
       tAATGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGcc                                                                                                                                                <  2:1642338/150‑1 (MQ=255)
            cccATTTTCATCCGCTAACTGCCCCGCTGTGACAGGTGCACGTATCGCGACAGGACGTAGCGTTAATGCGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATATGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGtt                                                                                                                                           <  1:1287267/150‑1 (MQ=255)
             ccATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTg                                                                                                                                          >  2:1067756/1‑150 (MQ=255)
              cATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTgg                                                                                                                                         <  2:975200/150‑1 (MQ=255)
               aTTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTACCGTGCATGTGCCCGCAGGTTGCGGTTGTGCAGTGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGt                                                                                                                                        <  2:284376/150‑1 (MQ=255)
                ttttCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTa                                                                                                                                       >  1:155817/1‑150 (MQ=255)
                   tCAACCGCTACCTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAgg                                                                                                                                    <  1:1068231/150‑1 (MQ=255)
                   tCAACAGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACTTGAAGTAGCGTTATTGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAgg                                                                                                                                    <  2:957062/150‑1 (MQ=255)
                    cAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGga                                                                                                                                   >  2:416599/1‑150 (MQ=255)
                    cAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGga                                                                                                                                   >  2:228569/1‑150 (MQ=255)
                      aCCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGgaga                                                                                                                                 >  1:508510/1‑150 (MQ=255)
                        cgcataCTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGt                                                                                                                               <  2:199282/145‑1 (MQ=255)
                        cGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGt                                                                                                                               >  2:530835/1‑150 (MQ=255)
                        cGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGt                                                                                                                               <  2:1393095/150‑1 (MQ=255)
                          cTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGATTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTgg                                                                                                                             >  2:1315827/1‑150 (MQ=255)
                               tcgcccGCAGTGACAGGTGCACGTAGCGCGACAGGTCTTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCTGATGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGTAGAGTGGCAGCg                                                                                                                        <  2:915859/147‑1 (MQ=255)
                                 cccaGCAGTGACAGGTGCACGTAGCGCGACAGGACTTAGCGGTCATGTGCCAGCAGGTTGCGGTTGTGCAGGTGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGc                                                                                                                      <  2:482860/150‑1 (MQ=255)
                                  cccGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGCGTGGCATCGGcc                                                                                                                     >  2:152303/1‑149 (MQ=255)
                                   ccGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGcac                                                                                                                    >  1:604984/1‑150 (MQ=255)
                                   ccGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGcac                                                                                                                    >  1:543174/1‑150 (MQ=255)
                                     gCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCacaa                                                                                                                  <  2:607263/150‑1 (MQ=255)
                                        gTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCacaaaca                                                                                                               >  2:117325/1‑150 (MQ=255)
                                           aCAGGTGCACGTAGCGCGACAGGCAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGCGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCAGCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAAcc                                                                                                            <  1:883333/150‑1 (MQ=255)
                                                tGCTCGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTc                                                                                                       <  2:157426/150‑1 (MQ=255)
                                                tGCACGGAGCGCGACAGGAAGGAGCTTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGATCATCGCTGCCCGGTTGGTAAGGAGAGTGGAATCGGCACAAACAACCAGTTc                                                                                                       <  2:1043339/150‑1 (MQ=255)
                                                 gCACGTAGCGCGACAGGACGTAGCGTTAATGAGCCCGCAGTAGGCGGTTGTGCAGGGGAGTGGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCGGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCg                                                                                                      <  1:838182/150‑1 (MQ=255)
                                                 gCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCg                                                                                                      <  2:1488854/150‑1 (MQ=255)
                                                  cACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGa                                                                                                     >  1:1724119/1‑150 (MQ=255)
                                                    cGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTTGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGccc                                                                                                   >  2:1113897/1‑150 (MQ=255)
                                                       aGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGc                                                                                                <  2:560781/150‑1 (MQ=255)
                                                                ggAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGc                                                                                       >  2:1677070/1‑150 (MQ=255)
                                                                   actAGCGTTAATGTGCCCGCAGGTTGCGGGTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCtgt                                                                                    <  1:543605/148‑1 (MQ=255)
                                                                      aGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCa                                                                                 >  2:1140/1‑150 (MQ=255)
                                                                      aGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCa                                                                                 >  2:685204/1‑150 (MQ=255)
                                                                        cGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAg                                                                               >  1:95563/1‑150 (MQ=255)
                                                                        cGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAg                                                                               >  2:1302242/1‑150 (MQ=255)
                                                                            aaTGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTgg                                                                           >  1:565561/1‑150 (MQ=255)
                                                                              tgtgCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGAc                                                                         <  1:630231/150‑1 (MQ=255)
                                                                              tgtgCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGAc                                                                         >  1:186688/1‑150 (MQ=255)
                                                                                 gCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAAc                                                                      >  1:1355332/1‑150 (MQ=255)
                                                                                 gCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAAc                                                                      >  2:620613/1‑150 (MQ=255)
                                                                                 gCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAAc                                                                      >  1:1326203/1‑150 (MQ=255)
                                                                                  cccGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACt                                                                     >  1:1326987/1‑150 (MQ=255)
                                                                                        gTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTcc                                                               >  1:449551/1‑150 (MQ=255)
                                                                                        gTTTGCGGTTGTGCAGGGGAGTTGAGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTcc                                                               >  2:1040952/1‑150 (MQ=255)
                                                                                          ttGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCgg                                                             <  1:502607/150‑1 (MQ=255)
                                                                                            gCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGc                                                           >  1:489004/1‑150 (MQ=255)
                                                                                             cGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGcc                                                          <  1:936022/150‑1 (MQ=255)
                                                                                             cGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGcc                                                          >  2:1545271/1‑150 (MQ=255)
                                                                                             cGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGc                                                           >  1:1493713/1‑149 (MQ=255)
                                                                                              ggTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCa                                                         >  1:1076525/1‑150 (MQ=255)
                                                                                                 tgtgCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATc                                                      <  2:1699928/150‑1 (MQ=255)
                                                                                                  gtgCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGGTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCACTcc                                                     <  2:693695/150‑1 (MQ=255)
                                                                                                   tgCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATccc                                                    >  2:886473/1‑150 (MQ=255)
                                                                                                    gCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCcg                                                   >  1:421179/1‑150 (MQ=255)
                                                                                                     cAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCcgg                                                  >  2:809367/1‑150 (MQ=255)
                                                                                                      aGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCcggc                                                 <  2:1326987/150‑1 (MQ=255)
                                                                                                       ggggAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCcggcg                                                >  2:999360/1‑150 (MQ=255)
                                                                                                         ggAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCggg                                              >  1:1273449/1‑150 (MQ=255)
                                                                                                         ggAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGATCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCggg                                              <  2:573798/150‑1 (MQ=255)
                                                                                                               gggggAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGcc                                        <  1:1240261/150‑1 (MQ=255)
                                                                                                                 gggAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCACGGTTGGTAAGGTGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCGACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCgg                                      <  2:1116140/150‑1 (MQ=255)
                                                                                                                  ggAATAAGGGCGGAGGGAGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCggg                                     <  2:449664/150‑1 (MQ=255)
                                                                                                                    aaTAAGGGCGGAGGGTGAGCGGCAAACCGCGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTGCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCGGGat                                   <  1:152303/150‑2 (MQ=255)
                                                                                                                    aaTAAGGGCGGAGGGTGAGCGGCAAACCGAGCAGCGCTGCCCGTTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAg                                   <  1:1060829/150‑1 (MQ=255)
                                                                                                                     aTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCGGGatt                                  <  2:1466680/150‑3 (MQ=255)
                                                                                                                     aTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCGGGatt                                  <  2:421179/150‑3 (MQ=255)
                                                                                                                       aaGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGa                                >  2:1545079/1‑150 (MQ=255)
                                                                                                                           gCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATcac                            <  1:1372901/150‑1 (MQ=255)
                                                                                                                           gCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATcac                            <  1:694980/150‑1 (MQ=255)
                                                                                                                            cGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATcaca                           >  2:1556818/1‑150 (MQ=255)
                                                                                                                            cGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATcaca                           <  1:886473/150‑1 (MQ=255)
                                                                                                                            cGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATcaca                           <  1:730672/150‑1 (MQ=255)
                                                                                                                                       cGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGt                <  1:1023655/150‑1 (MQ=255)
                                                                                                                                         gCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGttt              >  1:1130674/1‑150 (MQ=255)
                                                                                                                                          ccaaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACCCGTTGGGtttt             >  2:837110/3‑150 (MQ=255)
                                                                                                                                          ccaaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGtttt             >  2:1625120/3‑150 (MQ=255)
                                                                                                                                          cAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGtttt             <  2:977841/150‑1 (MQ=255)
                                                                                                                                          cAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGtttt             <  2:1719702/150‑1 (MQ=255)
                                                                                                                                           caaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTa            >  1:1594450/2‑150 (MQ=255)
                                                                                                                                           aaaCCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACCCGTTGGGTTTTa            >  2:1469931/1‑150 (MQ=255)
                                                                                                                                            aaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCGGGGAGTGATCACAACACGATGGGGTTTaa           >  1:1643196/1‑150 (MQ=255)
                                                                                                                                            aaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTaa           >  1:1534378/1‑150 (MQ=255)
                                                                                                                                            aaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTaa           >  2:1365249/1‑150 (MQ=255)
                                                                                                                                             aCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACATCCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAc          >  2:144106/1‑150 (MQ=255)
                                                                                                                                             aCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCCCGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACCCGTTGGGTTTTAAc          >  1:486188/1‑150 (MQ=255)
                                                                                                                                             aCCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACCCGTTGGGGTTTAAc          >  1:353202/1‑150 (MQ=255)
                                                                                                                                              ccGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAcc         <  1:1545271/150‑1 (MQ=255)
                                                                                                                                              ccGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGGGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAcc         <  2:1055705/150‑1 (MQ=255)
                                                                                                                                                gAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCAt       >  1:387611/1‑150 (MQ=255)
                                                                                                                                                 aGCATCGCTGCCCGGTTGGTATGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCAtt      <  2:218324/150‑1 (MQ=255)
                                                                                                                                                 aGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCAtt      >  1:1177451/1‑150 (MQ=255)
                                                                                                                                                   cATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAg    >  2:1158551/1‑150 (MQ=255)
                                                                                                                                                   cATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAg    >  1:1064943/1‑150 (MQ=255)
                                                                                                                                                     tCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGtt  >  1:22290/1‑150 (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 < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: