Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,637,357 T→C 6.9% intergenic (+24/+42) pdxJ → / ← PP_1437 pyridoxine 5'‑phosphate synthase/heavy metal sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,637,3570TC 6.9% 252.8 / 4.2 85intergenic (+24/+42)pdxJ/PP_1437pyridoxine 5'‑phosphate synthase/heavy metal sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (43/36);  new base C (4/2);  total (47/38)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GATCAAGGGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACC  >  NC_002947/1637216‑1637505
                                                                                                                                             |                                                                                                                                                    
gATCAAGGGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTcgcg                                                                                                                                                >  1:466472/1‑148 (MQ=255)
gATCAAGGGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTcgcg                                                                                                                                                <  2:466472/148‑1 (MQ=255)
      gggTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGcc                                                                                                                                        >  1:1315645/1‑150 (MQ=255)
          aTCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCcgcg                                                                                                                                    <  1:178453/150‑1 (MQ=255)
           tCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGa                                                                                                                                   >  1:514776/1‑150 (MQ=255)
           tCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGa                                                                                                                                   >  1:1595061/1‑150 (MQ=255)
            cAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGaa                                                                                                                                  <  2:514776/150‑1 (MQ=255)
             aaTGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGTCTTTTCGCGGGCAAGCCCGCGAag                                                                                                                                 >  1:968478/1‑150 (MQ=255)
             aaTGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCAGGCAAGCCCGCGaaa                                                                                                                                 >  1:143899/1‑149 (MQ=255)
               tGAGTTGAACATTGGCCCTGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCAGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAagag                                                                                                                               <  1:33725/150‑1 (MQ=255)
               tGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAagag                                                                                                                               <  2:96250/150‑1 (MQ=255)
               tGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAagag                                                                                                                               <  2:1595061/150‑1 (MQ=255)
                 aGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGc                                                                                                                             <  1:147058/150‑1 (MQ=255)
                  gTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGcc                                                                                                                            <  2:618454/150‑1 (MQ=255)
                  gTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGcc                                                                                                                            >  2:1512366/1‑150 (MQ=255)
                   ttGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGTCTTTTCGCGGGCAAGCCCGCGAAGAGGCCa                                                                                                                           <  2:968478/150‑1 (MQ=255)
                   ttGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTAGCGGGCAAGCCCGCGAAGAGGCCa                                                                                                                           >  2:512587/1‑150 (MQ=255)
                    tGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAg                                                                                                                          <  2:1398197/150‑1 (MQ=255)
                      aaCATTGGCCATGCACTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGc                                                                                                                         <  2:1802978/149‑1 (MQ=255)
                      aaCATTGGCCATGCACTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGc                                                                                                                         >  1:1802978/1‑149 (MQ=255)
                       aCATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGcac                                                                                                                       >  2:922556/1‑150 (MQ=255)
                        cATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAaga                                                                                                                                <  2:1642349/140‑1 (MQ=255)
                        cATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAaga                                                                                                                                >  1:1642349/1‑140 (MQ=255)
                        cATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGcaca                                                                                                                      >  1:750760/1‑150 (MQ=255)
                          ttGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAgg                                                                                                                    >  2:1571866/1‑150 (MQ=255)
                          ttGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAgg                                                                                                                    >  2:1075242/1‑150 (MQ=255)
                           tGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAggc                                                                                                                   <  2:942001/150‑1 (MQ=255)
                             gCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAggcgg                                                                                                                 >  2:1773156/1‑150 (MQ=255)
                              ccATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGc                                                                                                                             <  2:931977/137‑1 (MQ=255)
                              ccATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGc                                                                                                                             >  1:931977/1‑137 (MQ=39)
                                 tGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCCGCACAGGCGGCCGa                                                                                                             >  1:933478/1‑150 (MQ=255)
                                   cgcTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCTCTTTCTATCCCGCGAAGAGGCCAGCACAGGCGGCCGATc                                                                                                           >  2:547621/1‑150 (MQ=255)
                                       ggtggCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCTGCCGATCACTc                                                                                                       <  2:323896/150‑1 (MQ=255)
                                                    ttgtttgtCGGCGTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACAcc                                                                                          >  2:505666/1‑150 (MQ=255)
                                                    tggcttgtCGGCTTCAAGGCTGCGGGGGGCGAGATGAAGGCGCTGATCGTGGCGGCGTCGCGCGGAGTGTGGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACCGGCGGCCGATCACCCGGCCGGAAACAcc                                                                                          >  2:115194/5‑150 (MQ=38)
                                                      gtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAg                                                                                        <  2:1785682/150‑1 (MQ=255)
                                                       tttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAgg                                                                                       <  1:1512366/150‑1 (MQ=255)
                                                       tttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAgg                                                                                       <  1:922556/150‑1 (MQ=255)
                                                         tgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTcgcg                                                                                                                                                >  1:792685/1‑91 (MQ=255)
                                                         tgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTcgcg                                                                                                                                                <  2:792685/91‑1 (MQ=255)
                                                              gCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAAcc                                                                                >  1:668381/1‑150 (MQ=255)
                                                                 tCAAGGCTGCTGTGGCCGAGATGAAGGCTCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCggg                                                                             <  1:1100575/150‑1 (MQ=255)
                                                                 tCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCggg                                                                             <  1:1757955/150‑1 (MQ=255)
                                                                   aaGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGTTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTa                                                                           >  2:167203/1‑150 (MQ=255)
                                                                      gctgctGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGc                                                                        <  1:547621/150‑1 (MQ=255)
                                                                        tgctgTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCcgc                                                                                                                                     <  2:1121419/87‑1 (MQ=255)
                                                                        tgctgTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCcgc                                                                                                                                     >  1:1121419/1‑87 (MQ=255)
                                                                           tgtgGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAggcg                                                                                                                  >  2:1752906/1‑103 (MQ=255)
                                                                           tgtgGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAggcg                                                                                                                  <  1:1752906/103‑1 (MQ=255)
                                                                           tgtgGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAg                                                                   <  2:668381/150‑1 (MQ=255)
                                                                            gtgGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAgg                                                                  <  1:1669490/150‑1 (MQ=255)
                                                                               gCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAggggg                                                               >  2:441785/1‑150 (MQ=255)
                                                                                 cGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCa                                                             <  2:750760/150‑1 (MQ=255)
                                                                                  gagaTGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATAACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAg                                                            <  2:1709994/150‑1 (MQ=255)
                                                                                     aTGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctg                                                         <  1:441785/150‑1 (MQ=255)
                                                                                       gAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctgct                                                       >  1:1108138/1‑150 (MQ=255)
                                                                                             gcTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTc                                                 >  1:873691/1‑150 (MQ=255)
                                                                                                      ggcggcTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAggg                                                                 >  1:175665/1‑125 (MQ=255)
                                                                                                      ggcggcTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAggg                                                                 <  2:175665/125‑1 (MQ=255)
                                                                                                        cggcTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCt                                                  <  1:445066/138‑1 (MQ=255)
                                                                                                        cggcTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCt                                                  >  2:445066/1‑138 (MQ=255)
                                                                                                            ttCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATg                                  <  1:920949/150‑1 (MQ=255)
                                                                                                            ttCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACACGACCGCCATg                                  >  1:505151/1‑150 (MQ=255)
                                                                                                               gcgcTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGt                                                                            <  2:835687/105‑1 (MQ=255)
                                                                                                               gcgcTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGt                                                                            >  1:835687/1‑105 (MQ=255)
                                                                                                                     aTTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCggg                                                                             >  1:596781/1‑98 (MQ=255)
                                                                                                                     aTTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCggg                                                                             <  2:596781/98‑1 (MQ=255)
                                                                                                                     aTTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCAt                         >  1:1194560/1‑150 (MQ=255)
                                                                                                                         ttgGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATc                     >  1:1230496/1‑150 (MQ=255)
                                                                                                                          tgGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCg                    >  1:1056982/1‑150 (MQ=255)
                                                                                                                                 ctgtgctggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGACCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCAc             >  1:135675/8‑150 (MQ=25)
                                                                                                                                 cTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCAc             >  1:1813942/1‑150 (MQ=255)
                                                                                                                                     cccGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATc         <  2:1315645/149‑1 (MQ=255)
                                                                                                                                      ctggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGACCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCaa                               >  1:875693/3‑127 (MQ=25)
                                                                                                                                      ctggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGACCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCaa                               <  2:875693/125‑1 (MQ=25)
                                                                                                                                      ccGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATACACTTCAccagcg        >  2:1573557/1‑145 (MQ=255)
                                                                                                                                       cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGa                                                     <  1:856505/104‑1 (MQ=255)
                                                                                                                                       cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGa                                                     >  2:856505/1‑104 (MQ=255)
                                                                                                                                       cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGCTCGc       >  2:139534/1‑150 (MQ=255)
                                                                                                                                       cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGc       >  2:698439/1‑150 (MQ=255)
                                                                                                                                       cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGc       >  1:1760703/1‑150 (MQ=255)
                                                                                                                                        ggCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGa                       <  2:853955/133‑1 (MQ=255)
                                                                                                                                        ggCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGa                       >  1:853955/1‑133 (MQ=255)
                                                                                                                                         gCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGcca     >  2:642297/1‑150 (MQ=255)
                                                                                                                                          ccTCTTCGCGGGCTTGCCCGCGAAAAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCAt                         >  1:800830/1‑129 (MQ=39)
                                                                                                                                            ttttCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCg                                                >  2:733807/1‑104 (MQ=255)
                                                                                                                                            ttttCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCg                                                <  1:733807/104‑1 (MQ=255)
                                                                                                                                            ttttCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaacc  >  2:449818/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GATCAAGGGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACC  >  NC_002947/1637216‑1637505

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