Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,642,335 C→T 8.7% intergenic (+53/+13) 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,642,3350CT 8.7% 955.7 / 56.6 321intergenic (+53/+13)pdxJ/PP_1437pyridoxine 5'‑phosphate synthase/heavy metal sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (145/148);  new base T (20/8);  total (165/156)
Fisher's exact test for biased strand distribution p-value = 2.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 328 total aligned reads displayed.

GAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAAC  >  NC_002947/1642186‑1642484
                                                                                                                                                     |                                                                                                                                                     
 aaCATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGca                                                                                                                                                      >  2:6879/1‑150 (MQ=255)
   cATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGcaca                                                                                                                                                    <  1:2456838/150‑1 (MQ=255)
     ttGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAgg                                                                                                                                                  <  2:4028273/150‑1 (MQ=255)
     ttGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAgg                                                                                                                                                  <  2:1380447/150‑1 (MQ=38)
      tGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAggc                                                                                                                                                 <  2:2884173/150‑1 (MQ=37)
        gCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAggcgg                                                                                                                                               >  1:6909768/1‑150 (MQ=255)
        gCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTAGCGGGCAAGACCGCGACGAGGCCCGCACAggcgg                                                                                                                                               >  1:3203136/1‑150 (MQ=255)
         ccATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAggcggc                                                                                                                                              >  1:1647580/1‑150 (MQ=255)
         ccATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAggcggc                                                                                                                                              <  1:6600276/150‑1 (MQ=255)
          cATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGGGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGcc                                                                                                                                             <  2:8636946/150‑1 (MQ=255)
              cgcTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATc                                                                                                                                         >  2:8393912/1‑150 (MQ=255)
               gcTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCa                                                                                                                                        <  2:4573859/150‑1 (MQ=255)
                cTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTCGGCCGCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGAccaacaat                                                                                                                                       <  1:6230656/150‑7 (MQ=18)
                      gcgcATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCa                                                                                                                                 <  2:2630788/150‑1 (MQ=255)
                       cgcATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCATTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAg                                                                                                                                <  2:4231194/150‑1 (MQ=255)
                        gcATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAgg                                                                                                                               >  2:1805054/1‑150 (MQ=255)
                         cATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGa                                                                                                                              >  1:7619414/1‑150 (MQ=255)
                         cATGCCTCGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGa                                                                                                                              <  2:6474176/150‑1 (MQ=255)
                              cTTGTTGGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAacac                                                                                                                         <  1:6053674/150‑1 (MQ=255)
                                tgtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCa                                                                                                                       >  1:7896076/1‑150 (MQ=255)
                                tgtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCa                                                                                                                       <  1:6376426/150‑1 (MQ=255)
                                tgtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCa                                                                                                                       <  1:4954689/150‑1 (MQ=255)
                                  tttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAgg                                                                                                                     <  1:4533617/150‑1 (MQ=255)
                                  tttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAgg                                                                                                                     <  1:6723242/150‑1 (MQ=255)
                                  tttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAgg                                                                                                                     <  2:2546551/150‑1 (MQ=255)
                                      tCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTaa                                                                                                                 >  1:4019419/1‑150 (MQ=255)
                                      tCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTaa                                                                                                                 >  2:1736357/1‑150 (MQ=255)
                                        ggCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAAc                                                                                                               <  2:2045703/150‑1 (MQ=255)
                                         gCTTCAAGGCTACTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAAcc                                                                                                              >  1:7361510/1‑150 (MQ=255)
                                            tCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCggg                                                                                                           <  1:4128488/150‑1 (MQ=255)
                                              aaGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTa                                                                                                         >  2:5634537/1‑150 (MQ=255)
                                                ggCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAgg                                                                                                       <  2:4520895/150‑1 (MQ=255)
                                                 gctgctGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGc                                                                                                      >  1:3808949/1‑150 (MQ=255)
                                                 gctgctGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGc                                                                                                      <  1:612603/150‑1 (MQ=255)
                                                      tgtgGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGTGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAg                                                                                                 >  1:7703888/1‑150 (MQ=255)
                                                      tgtgGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAg                                                                                                 <  1:8750151/150‑1 (MQ=255)
                                                           ccGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGc                                                                                            <  1:711020/150‑1 (MQ=255)
                                                           ccGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGc                                                                                            <  1:1864112/150‑1 (MQ=255)
                                                           ccGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGc                                                                                            >  1:1614923/1‑150 (MQ=255)
                                                            cGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCa                                                                                           >  2:8752062/1‑150 (MQ=255)
                                                               gaTGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgct                                                                                        <  2:624217/150‑1 (MQ=255)
                                                                 tGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctgc                                                                                      <  1:3400324/150‑1 (MQ=255)
                                                                     ggCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGAc                                                                                  >  2:7922322/1‑150 (MQ=255)
                                                                      gcgcTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCCGCAGCTGAcc                                                                                 >  1:4792044/1‑150 (MQ=255)
                                                                      gcgcTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGAcc                                                                                 >  2:3585040/1‑150 (MQ=255)
                                                                       cgcTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCt                                                                                <  2:3905204/150‑1 (MQ=255)
                                                                        gcTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTc                                                                               <  1:4456058/150‑1 (MQ=255)
                                                                           gATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAc                                                                            <  1:2724826/150‑1 (MQ=255)
                                                                             tCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCACGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAccc                                                                          <  2:7245574/150‑1 (MQ=255)
                                                                             tCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCACGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAccc                                                                          <  2:118923/150‑1 (MQ=255)
                                                                                 ggcggcTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAcccgccc                                                                      <  2:4991444/150‑1 (MQ=255)
                                                                                  gcggcTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAcccgcccg                                                                     >  1:5856382/1‑150 (MQ=255)
                                                                                            gcTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACt                                                           >  1:1364579/1‑150 (MQ=255)
                                                                                            gcTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACt                                                           >  1:5440486/1‑150 (MQ=255)
                                                                                            gcTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACt                                                           >  1:508200/1‑150 (MQ=255)
                                                                                            gcTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCCCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTTGGCCCCAGGGGGCCGCTGCTGACCGCGACCCGCCCGCCAAGCCAca                                                           >  1:1011305/1‑149 (MQ=255)
                                                                                              tGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTAGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCGGCTGACCTCGACCCGCCCGCCATGCAACTGc                                                         >  1:1740230/1‑150 (MQ=255)
                                                                                                   gttgGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAGTTGCATggc                                                    <  1:1253672/150‑3 (MQ=255)
                                                                                                    ttgGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATc                                                   <  2:1651411/150‑1 (MQ=255)
                                                                                                    ttgGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATc                                                   <  2:2133236/150‑1 (MQ=255)
                                                                                                    ttgGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATc                                                   <  2:3203136/150‑1 (MQ=255)
                                                                                                       gccGCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGACCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCCTGCAACTGCATGATCGAc                                                >  1:7866449/3‑150 (MQ=25)
                                                                                                          gCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTc                                             >  1:2140664/1‑150 (MQ=255)
                                                                                                          gCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCAAGCAACAGCATGATCGACTTc                                             >  1:8162970/1‑150 (MQ=255)
                                                                                                           ccTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCa                                            <  2:7506461/150‑1 (MQ=255)
                                                                                                           ccTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCa                                            <  1:5709151/150‑1 (MQ=255)
                                                                                                             tgtgctggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGACCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACg                                          >  2:5307650/7‑150 (MQ=37)
                                                                                                              gTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGa                                         <  1:2886898/150‑1 (MQ=255)
                                                                                                                 ccGGCCTTTTCGCGGGCCAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACGCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCCGGGGGCAGCTGCTGACCGCGACCCGCCCGCCATGCAACTGCATGCTCGACTTCACGAGCg                                      >  2:4720353/1‑150 (MQ=255)
                                                                                                                 ccGGCCTTTTCGCGGGCAAGCCCGCGAATAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCg                                      <  1:559070/150‑1 (MQ=255)
                                                                                                                  cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCACCTGCATGATCGACTTCACGATCGc                                     >  2:5579557/1‑150 (MQ=255)
                                                                                                                  cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGc                                     <  1:7316875/150‑1 (MQ=255)
                                                                                                                  cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGc                                     >  2:2832995/1‑150 (MQ=255)
                                                                                                                  cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCAAAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGc                                     >  2:4854726/1‑150 (MQ=255)
                                                                                                                   ggCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGcc                                    >  2:913949/1‑150 (MQ=255)
                                                                                                                     ccTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaa                                  >  2:2245105/1‑150 (MQ=255)
                                                                                                                      ctcttCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaac                                 >  1:1148984/4‑150 (MQ=39)
                                                                                                                      ctcttCGCGGGCTTGCCCGCGAAAAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaac                                 >  1:2606115/4‑150 (MQ=255)
                                                                                                                       ttttCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaacc                                >  2:3641326/1‑150 (MQ=255)
                                                                                                                       ttttCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaacc                                >  2:7388281/1‑150 (MQ=255)
                                                                                                                       tcttCGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaacc                                >  2:4412101/3‑150 (MQ=255)
                                                                                                                       tcttCGCGGGCTTGCCCGCGAACAGGCCGGTACCGGCGACCGATCACTCGGCAGGACACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGcccacc                                >  2:50484/3‑150 (MQ=255)
                                                                                                                            gcgGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaaccccaac                           >  1:2998991/1‑150 (MQ=255)
                                                                                                                            gcgGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaaccccaac                           >  1:977416/1‑150 (MQ=255)
                                                                                                                                gCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCgg                       <  1:7808160/150‑1 (MQ=255)
                                                                                                                                      ccGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGc                 >  1:3538030/1‑150 (MQ=255)
                                                                                                                                       cgcgAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAACCACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGcc                >  1:2082422/1‑150 (MQ=255)
                                                                                                                                       cgcgAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGcc                >  2:7488273/1‑150 (MQ=255)
                                                                                                                                       cgcgAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACACCAACCCGGTCCCGcc                >  2:7047018/1‑150 (MQ=255)
                                                                                                                                        gcgAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTGGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGccc               >  2:3994051/1‑150 (MQ=255)
                                                                                                                                        gcgAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGccc               >  2:2994444/1‑150 (MQ=255)
                                                                                                                                         cgAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCt              <  2:6909768/150‑1 (MQ=255)
                                                                                                                                          gAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTc             <  1:8393912/150‑1 (MQ=255)
                                                                                                                                          gAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTc             >  2:2009307/1‑150 (MQ=255)
                                                                                                                                                gCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAggg                                                                                               >  2:7885248/1‑62 (MQ=255)
                                                                                                                                                gCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCg                                                                              <  2:1259540/79‑1 (MQ=255)
                                                                                                                                                gCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCcg       >  2:1543930/1‑150 (MQ=255)
                                                                                                                                                gCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCcg       >  2:5501245/1‑150 (MQ=255)
                                                                                                                                                  cAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCcgcg     >  2:7586796/1‑150 (MQ=255)
                                                                                                                                                    gCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGaa   <  2:6092551/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAAC  >  NC_002947/1642186‑1642484

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