Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,657,937 T→G 8.5% intergenic (‑38/+123) PP_4961 ← / ← PP_4962 lipoprotein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,657,9370TG 8.5% 811.2 / 42.5 266intergenic (‑38/+123)PP_4961/PP_4962lipoprotein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (99/142);  new base G (5/19);  total (105/161)
Fisher's exact test for biased strand distribution p-value = 7.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 277 total aligned reads displayed.

GGGTGTCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCT  >  NC_002947/5657788‑5658085
                                                                                                                                                     |                                                                                                                                                    
 ggTGTCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTtctc                                                                                                                                                     <  2:4165897/150‑1 (MQ=255)
   tgtCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCCAAAGCGCTTTCATTCCGTTGCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCgg                                                                                                                                                   >  2:5999103/1‑150 (MQ=255)
   tgtCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCgg                                                                                                                                                   <  1:1948306/150‑1 (MQ=255)
    gtCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCggg                                                                                                                                                  <  2:4156497/150‑1 (MQ=255)
    gtCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCggg                                                                                                                                                  <  2:1749278/150‑1 (MQ=255)
    gtCGCAGACCCAACGGGTCCAGTTATCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCggg                                                                                                                                                  <  2:4627922/150‑1 (MQ=255)
     tCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTtctc                                                                                                                                                     >  2:2813250/1‑146 (MQ=255)
     tCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGc                                                                                                                                                 <  1:5451252/150‑1 (MQ=255)
     tCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGc                                                                                                                                                 <  2:7183736/150‑1 (MQ=255)
     tCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGc                                                                                                                                                 <  2:5924234/150‑1 (MQ=255)
           aCCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGccc                                                                                                                                           >  2:3200589/1‑150 (MQ=255)
           aCCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCGCGGGCTCGccc                                                                                                                                           <  1:2067905/150‑1 (MQ=255)
           aCCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCGCGGGCTCGccc                                                                                                                                           >  1:3529992/1‑150 (MQ=255)
             ccAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCGGCCCTTCGCGGGCTCGCCCGc                                                                                                                                         >  2:326715/1‑150 (MQ=255)
              caccGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCGAGCCCgag                                                                                                                                        <  1:6676259/147‑3 (MQ=255)
               aaCGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTc                                                                                                                                       >  2:822764/1‑150 (MQ=255)
                aCGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTcc                                                                                                                                      <  1:3570778/150‑1 (MQ=255)
                aCGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTcc                                                                                                                                      <  1:7278311/150‑1 (MQ=255)
                  gggTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCca                                                                                                                                    <  2:4100791/150‑1 (MQ=255)
                   ggTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCcac                                                                                                                                   <  2:2138281/150‑1 (MQ=255)
                   ggTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCcac                                                                                                                                   <  1:4972668/150‑1 (MQ=255)
                    gTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCcaca                                                                                                                                  >  2:1501959/1‑150 (MQ=255)
                     tCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGACTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAg                                                                                                                                 <  1:1303853/150‑1 (MQ=255)
                      ccAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAgg                                                                                                                                <  1:6459458/150‑1 (MQ=255)
                      ccAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAga                                                                                                                                <  1:4907090/150‑2 (MQ=255)
                       cAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTAACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGt                                                                                                                               >  1:1418582/1‑150 (MQ=255)
                        aGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTa                                                                                                                              <  1:4705330/150‑1 (MQ=255)
                          ttGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACa                                                                                                                            <  2:634900/150‑1 (MQ=255)
                           tGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAg                                                                                                                           <  1:2708270/150‑1 (MQ=255)
                           tGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAg                                                                                                                           <  1:2843557/150‑1 (MQ=255)
                            gTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGt                                                                                                                          <  2:4539948/150‑1 (MQ=255)
                            gTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGt                                                                                                                          <  2:105611/150‑1 (MQ=255)
                               tGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTCCAGTGCa                                                                                                                       >  2:5212571/1‑150 (MQ=255)
                               tGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCa                                                                                                                       <  2:5083299/150‑1 (MQ=255)
                               tGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCa                                                                                                                       <  2:6552533/150‑1 (MQ=255)
                                 ccGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGt                                                                                                                     >  2:1195857/1‑150 (MQ=255)
                                       ccGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCtgt                                                                                                               <  1:3484917/150‑1 (MQ=255)
                                       ccGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCtgt                                                                                                               <  1:5759628/150‑1 (MQ=255)
                                       ccGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCtgt                                                                                                               <  1:5270752/150‑1 (MQ=255)
                                         ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                             <  2:2658180/150‑1 (MQ=255)
                                         ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                             <  2:1308927/150‑1 (MQ=255)
                                         ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                             <  2:3529992/150‑1 (MQ=255)
                                           ttGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGa                                                                                                           <  2:5524622/150‑1 (MQ=255)
                                            tGCGCCGGCTGCAGCTGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAg                                                                                                          <  2:7014726/150‑1 (MQ=255)
                                            tGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCCAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGACTCTGTTGGAg                                                                                                          >  1:4549418/1‑150 (MQ=255)
                                              cgcCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTGTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCGGCCCTTCGCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCg                                                                                                        <  2:7116871/150‑1 (MQ=21)
                                               gcCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCgg                                                                                                       >  2:3426986/1‑150 (MQ=255)
                                               gcCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCgg                                                                                                       <  1:5060122/150‑1 (MQ=25)
                                               gcCGGCAGCAGCAGGGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGGTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCgg                                                                                                       <  1:5650352/150‑1 (MQ=25)
                                                   gCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGa                                                                                                   <  1:4304186/150‑1 (MQ=255)
                                                   gCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGa                                                                                                   <  2:904666/150‑1 (MQ=255)
                                                          cagTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                             >  2:3176648/1‑133 (MQ=255)
                                                           agTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGa                                                                                           >  1:6073664/1‑150 (MQ=255)
                                                             tGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCgg                                                                                                       <  1:2226983/136‑1 (MQ=25)
                                                             tGACCATCCGGCCCGTGTCGCCAGGGCCGCCCCGGCAAAACGCGCTTTCATTCCGTTTCCCCGATTTGGGGCTGCTGTGCCGCCCTTCTCGGGCTCGCCCGCCCCCACAGGTACCGGGCTGTCTGGGTGGGAGCGGGCGAGCCCGCGAAg                                                                                         <  2:3720962/150‑1 (MQ=255)
                                                                 cATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCa                                                                                                                       >  1:2417040/1‑116 (MQ=255)
                                                                 cATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCa                                                                                                                       <  2:2417040/116‑1 (MQ=255)
                                                                 cATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGcc                                                                                     >  1:7287768/1‑150 (MQ=255)
                                                                  aTCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCg                                                                                    >  1:3894750/1‑150 (MQ=255)
                                                                     cTGCCAGTGTCGCCCGGGCCGTCACGGCCAAACGCGCTTTCATGCAGTTTCGTGGCTGGGGGGCTGCTGTGCAGCCCTGCTCGGGCGCGCCCGCTCCCCCCGGTACAGTGCCGTCTCTGTGGGCGGGGGCGAGCCCGCGAAGGGCCGcac                                                                                 <  2:6464880/150‑1 (MQ=255)
                                                                      tGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGcacc                                                                                >  1:3237851/1‑149 (MQ=255)
                                                                      tGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGcaca                                                                                >  1:3206934/1‑150 (MQ=255)
                                                                      tGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCCCGCCCACAGGTACACTGCACTCACAGGGGGAGCGGGCGAGCCCGCGAAGGGCCGcaca                                                                                >  1:2159587/1‑150 (MQ=255)
                                                                       gCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAg                                                                               >  1:6783406/1‑150 (MQ=255)
                                                                          aGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCgg                                                                            >  1:1801062/1‑150 (MQ=255)
                                                                           gtgtCGCCAGGGCCGTCACGGCCAAAAGCGCTTGCATTCAGTTCCCGCGCTTTGGGGCTGCTGTGCCGCCCTTCGCGGGCTCGCCCGCCCCCACCGGTACCGTGCCGTCTCGGTGGGAGGGGGCGAGCCCGCGAAGGGCCGCACAGCGGc                                                                           <  2:2173717/150‑1 (MQ=25)
                                                                             gtCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGa                                                                                                           <  1:434418/116‑1 (MQ=255)
                                                                              tCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGcccc                                                                        >  2:2776196/1‑150 (MQ=255)
                                                                                 ttttggCCGTGACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCcttc                                                                     >  1:3767983/5‑150 (MQ=255)
                                                                                        cGTCACGGCCCAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGa                                                              >  2:427077/1‑150 (MQ=255)
                                                                                        cGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCCGtct                                                                                                                   >  1:4462596/1‑97 (MQ=255)
                                                                                        cGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGtct                                                                                                                   >  1:748330/1‑97 (MQ=255)
                                                                                        cGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGtct                                                                                                                   <  2:748330/97‑1 (MQ=255)
                                                                                        cGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGTGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGa                                                              <  2:3608029/150‑1 (MQ=255)
                                                                                              ggCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACaa                                                        >  2:4401959/1‑150 (MQ=255)
                                                                                                       cgcTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACCAATCAACGtt                                               >  2:6791794/1‑150 (MQ=255)
                                                                                                         ctttatttcaTTTTCCTCGTTTGGGGGCGGCTGTGCAGCCATTCTGGGGCTCGCCCGCTCCCACCGGTACAGTGCAGTCTCTGTGGGAGGGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCa                                             <  1:1184140/144‑1 (MQ=255)
                                                                                                             ctttcaGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCGCGGGCTCGCCCGCTCCCACCGGTACAGTGCTGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGtt                                         <  2:6316759/148‑1 (MQ=255)
                                                                                                              attcaGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCCCAACCACACGTACACTGCCGTCTCTGGGGGAGCGGGCGAGCCCCCCCAGGGCCGCCCCGCGGCCCCTTCGTCAAGACCACACATCAACGGTAAGGTTg                                        >  1:6250796/1‑150 (MQ=255)
                                                                                                               ttcaGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGc                                       >  1:593474/1‑150 (MQ=255)
                                                                                                                 caGTTTCCTCGATTTGGGGCTGCTGTGCCGCCCTTCTCGGGCTCGCCCGCCCCCACAGGTACAGTGCTGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCt                                     <  2:1944359/150‑1 (MQ=255)
                                                                                                                    tttCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTAtt                                  >  2:220183/1‑150 (MQ=255)
                                                                                                                        cTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCAGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAccc                              >  2:5708089/1‑150 (MQ=255)
                                                                                                                         tCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACCGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCGTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCg                             >  2:5428727/1‑150 (MQ=255)
                                                                                                                            gagaagggCTGCACAGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAAACAACGTTCAGGTTGCCTATTACCCGATc                          >  1:5332049/6‑150 (MQ=38)
                                                                                                                            aTTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCCCAAATCAACGTTCAGGTTGCCTATTACCCGATc                          >  2:3935287/1‑150 (MQ=255)
                                                                                                                               tGGGGCTGCGGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACCGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAg                       <  2:3877560/150‑1 (MQ=255)
                                                                                                                                ggggCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAgg                      >  2:4852713/1‑150 (MQ=255)
                                                                                                                                ggggCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAgg                      >  2:1190396/1‑150 (MQ=255)
                                                                                                                                ggggCTGCTGTGCAGCCCTTCGCGGGCGCCCCCGCCCCCCCCGGTACCGTGCTGTCTGTGTGGGCGGGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAgg                      <  1:822764/150‑1 (MQ=255)
                                                                                                                                 gggCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGTTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAggg                     <  2:993574/150‑1 (MQ=255)
                                                                                                                                   gctgctGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCa                   >  2:29236/1‑150 (MQ=255)
                                                                                                                                      ggtgtgCAGCCCTTCTGGGGCTCGCCCGCTCCCACTGTTACCGTGCTGTCTCTGTGGGAGGGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAAcc                <  1:1501959/148‑1 (MQ=255)
                                                                                                                                        tgtgCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCGGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCg              >  1:1757709/1‑150 (MQ=255)
                                                                                                                                            cAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGgcgc          <  2:5487536/150‑1 (MQ=255)
                                                                                                                                               cccTTCTCGGGCTCGCCCGCTCCCACAGGTACTGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCtt       <  1:6029617/150‑1 (MQ=255)
                                                                                                                                               cccTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAc                                <  2:2363166/125‑1 (MQ=255)
                                                                                                                                                ccTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTg      <  2:1115846/150‑1 (MQ=255)
                                                                                                                                                 cTTCTGGGGCTCGCCCGCGCCCCCTGGTACAGTGCAGTCTCGGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTgc     <  1:4098975/150‑1 (MQ=255)
                                                                                                                                                  ttCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCGTgcg    >  2:2537296/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGTGTCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCT  >  NC_002947/5657788‑5658085

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