Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 308,740 C→T 7.0% intergenic (+1754/‑116) hslO → / → PP_0254 molecular chaperone HslO/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947308,7400CT 7.0% 251.5 / 7.0 85intergenic (+1754/‑116)hslO/PP_0254molecular chaperone HslO/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/40);  new base T (2/4);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

GGCTGACAAAGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGACCGAAGCCCCTCAGCGCACTGCCTTCTCC  >  NC_002947/308593‑308888
                                                                                                                                                   |                                                                                                                                                    
ggCTGACAAAGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAg                                                                                                                                                    >  1:117125/1‑150 (MQ=255)
 gCTGACAAAGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAgg                                                                                                                                                   >  2:370360/1‑150 (MQ=255)
    gACAAAGCGGCCTACGACGCGGCTGCCAAAGGTCTGGCCAAGCTGTTCAGCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGggc                                                                                                                                                <  2:424673/150‑1 (MQ=255)
      cAAAGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAgggg                                                                                                                                                 <  1:300814/147‑1 (MQ=255)
      cAAAGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAgggg                                                                                                                                                 >  2:300814/1‑147 (MQ=255)
      cAAAGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGggcgg                                                                                                                                              <  2:501793/150‑1 (MQ=255)
       aaaGCGGCCTACGACGAGGCTGCCAACCTGCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGggcggc                                                                                                                                             >  1:1322552/1‑150 (MQ=255)
       aaaGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGggcggc                                                                                                                                             >  2:845211/1‑150 (MQ=255)
       aaaGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCAGCAAAGAGAAAGCCGACCCCAGGggcggc                                                                                                                                             >  1:1366727/1‑150 (MQ=255)
        aaGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCTGGGGGCGGCt                                                                                                                                            <  2:274177/150‑1 (MQ=255)
         aGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCtt                                                                                                                                           <  2:1555155/150‑1 (MQ=255)
          gCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCTGGGGGCGGCttt                                                                                                                                          <  1:1936886/150‑1 (MQ=255)
           cGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCtttt                                                                                                                                         <  2:117125/150‑1 (MQ=255)
           cGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCtttt                                                                                                                                         <  2:1097560/150‑1 (MQ=255)
               cTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTtgt                                                                                                                                     <  2:163454/150‑1 (MQ=255)
               cTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTtgt                                                                                                                                     <  1:599823/150‑1 (MQ=255)
                     cgaGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCACGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGTCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTa                                                                                                                               <  2:115345/150‑1 (MQ=255)
                      gaGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAg                                                                                                                              <  2:839917/150‑1 (MQ=255)
                       aGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAgg                                                                                                                             <  1:845211/150‑1 (MQ=255)
                        ggCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAggg                                                                                                                            >  1:343525/1‑150 (MQ=255)
                             ccAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCtt                                                                                                                         <  2:772140/148‑1 (MQ=255)
                             ccAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCtt                                                                                                                         >  1:772140/1‑148 (MQ=255)
                               aaaGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGccc                                                                                                                     >  1:1162199/1‑150 (MQ=255)
                                   gTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGa                                                                                                                 >  1:509361/1‑150 (MQ=255)
                                     cTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAg                                                                                                               <  2:1366727/150‑1 (MQ=255)
                                        gCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAgggg                                                                                                            >  2:460141/1‑150 (MQ=255)
                                          cAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGcc                                                                                                          >  1:687106/1‑150 (MQ=255)
                                            aGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCg                                                                                                        >  2:360810/1‑150 (MQ=255)
                                             gCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGggc                                                                                                                                                <  2:1691951/109‑1 (MQ=255)
                                             gCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCGCAGCCGCAACGAGAAAGCCGCCCCCAGGggc                                                                                                                                                >  1:1691951/1‑109 (MQ=255)
                                                    aTCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCa                                                                                                <  1:460141/150‑1 (MQ=255)
                                                     tCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAg                                                                                               <  2:687106/150‑1 (MQ=255)
                                                       gAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCg                                                                                             >  2:194554/1‑150 (MQ=255)
                                                             ttCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTc                                                                                       >  2:314586/1‑150 (MQ=255)
                                                             ttCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTc                                                                                       >  1:694792/1‑150 (MQ=255)
                                                              tCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCt                                                                                      <  2:1162199/150‑1 (MQ=255)
                                                               cAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTg                                                                                     <  1:370360/150‑1 (MQ=255)
                                                               cAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTg                                                                                     >  2:1002249/1‑150 (MQ=255)
                                                               cAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTg                                                                                     <  2:694792/150‑1 (MQ=255)
                                                                  gaagTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATCGGCAGCGAAAGTCTGGtt                                                                                  >  2:1563735/1‑150 (MQ=255)
                                                                       tCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCt                                                                             <  2:757278/150‑1 (MQ=255)
                                                                        cGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTgg                                                                                                                                   >  2:1800610/1‑95 (MQ=255)
                                                                        cGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTgg                                                                                                                                   <  1:1800610/95‑1 (MQ=255)
                                                                          aaGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGtt                                                                          <  2:1322552/150‑1 (MQ=255)
                                                                          aaGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGAGATCTTAATATCGGACGTGTGTTGTCCGCTAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGtt                                                                          >  1:228383/1‑150 (MQ=255)
                                                                                ccGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTg                                                                    <  1:682983/150‑1 (MQ=255)
                                                                                  gACGCCATCAAGGCCGCTTGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTTAGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCACAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTc                                                                  <  1:1106350/150‑1 (MQ=255)
                                                                                              gCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGaaa                                                                                          >  2:1091986/1‑114 (MQ=255)
                                                                                              gCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGaaa                                                                                          <  1:1091986/114‑1 (MQ=255)
                                                                                               ccGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATcgcg                                                     <  2:1740975/150‑1 (MQ=255)
                                                                                                 gCTGGCCCGCAGCTGTAAGCCGAGCTCAGGCGAACAGCAGAGGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGaa                                                   <  2:155100/150‑1 (MQ=255)
                                                                                                 gCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGaa                                                   >  2:1665507/1‑150 (MQ=255)
                                                                                                   tGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGa                                                 >  1:218791/1‑150 (MQ=255)
                                                                                                    ggCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGAc                                                <  1:1283431/150‑1 (MQ=255)
                                                                                                       ccGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGaa                                                                                                                >  1:232127/1‑83 (MQ=255)
                                                                                                       ccGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGaa                                                                                                                <  2:232127/83‑1 (MQ=255)
                                                                                                       ccGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGc                                                                                              >  2:831334/1‑101 (MQ=255)
                                                                                                       ccGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGc                                                                                              <  1:831334/101‑1 (MQ=255)
                                                                                                          cAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGagc                                          <  2:218791/150‑1 (MQ=255)
                                                                                                           aGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGagcc                                         >  1:381584/1‑150 (MQ=255)
                                                                                                              tGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTtgt                                                                                                                                     <  1:1294246/55‑1 (MQ=255)
                                                                                                              tGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTtgt                                                                                                                                     >  2:1294246/1‑55 (MQ=255)
                                                                                                               gTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGTGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGagccagcc                                     >  1:1845907/1‑150 (MQ=255)
                                                                                                                  aGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGagc                                          <  1:792105/142‑1 (MQ=255)
                                                                                                                  aGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGagc                                          >  2:792105/1‑142 (MQ=255)
                                                                                                                  aGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGa                                  >  1:637778/1‑150 (MQ=255)
                                                                                                                   gCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGAt                                 <  1:1665507/150‑1 (MQ=255)
                                                                                                                     cGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGtt                                                                          >  2:83406/1‑107 (MQ=255)
                                                                                                                     cGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGtt                                                                          <  1:83406/107‑1 (MQ=255)
                                                                                                                         cTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCg                                                                                             <  2:299656/84‑1 (MQ=255)
                                                                                                                         cTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCg                                                                                             >  1:299656/1‑84 (MQ=255)
                                                                                                                           cAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGACCGAAg                         <  1:1002249/150‑1 (MQ=255)
                                                                                                                                gCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGACCGAAGCCCCt                    <  1:360810/150‑1 (MQ=255)
                                                                                                                                  agcGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGACCGAAGCCCCTCa                  >  2:1837244/3‑150 (MQ=255)
                                                                                                                                   acaaaaaaGCCGCCCCTGGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCa                                                                 >  2:563133/6‑102 (MQ=38)
                                                                                                                                   acaaaaaaGCCGCCCCTGGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCa                                                                 <  1:563133/97‑1 (MQ=38)
                                                                                                                                     gagaAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGCGCCCGCATAGGCAGCGAAAGTCTGTTTTTTCTATTTCACTGTTATCATCCATTGCGGAGACAAGAGCCAGCCAGATGACCGAAGCCTCTCAgcg               >  2:1154407/1‑150 (MQ=255)
                                                                                                                                     aaaaaaGCCGCCCCTGGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGACCGAAGCCCCTCAg                 <  2:694200/145‑1 (MQ=255)
                                                                                                                                     aaaaaaGCCGCCCCTGGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGACCGAAGCCCCTCAg                 >  1:694200/4‑148 (MQ=255)
                                                                                                                                          aGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATa                                                                                                    >  2:1598620/1‑60 (MQ=255)
                                                                                                                                          aGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATa                                                                                                    <  1:1598620/60‑1 (MQ=255)
                                                                                                                                          aGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTc                                                                              >  2:29824/1‑82 (MQ=255)
                                                                                                                                          aGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTc                                                                              <  1:29824/82‑1 (MQ=255)
                                                                                                                                          aGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGACCGAAGCCCCTCAGCGCACTg          >  2:567451/1‑150 (MQ=255)
                                                                                                                                            ccgccCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGACCGAAGCCCCTCAGCGCACTGcc        >  2:1457140/1‑150 (MQ=255)
                                                                                                                                                 cccAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGACCGAAGCCCCTCAGCGCACTGCCTtctc   <  2:228383/150‑1 (MQ=255)
                                                                                                                                                  ccAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGACCGAAGCCCCTCAGCGCACTGCCTTCTcc  <  1:194554/150‑1 (MQ=255)
                                                                                                                                                  ccAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGACCGAAGCCCCTCAGCGCACTGCCTTCTcc  >  1:635886/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCTGACAAAGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGACCGAAGCCCCTCAGCGCACTGCCTTCTCC  >  NC_002947/308593‑308888

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: