Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,863,153 T→C 6.8% intergenic (‑137/+100) ligA ← / ← PP_4275 DNA ligase/cell division protein ZipA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,863,1530TC 6.8% 450.6 / 12.4 148intergenic (‑137/+100)ligA/PP_4275DNA ligase/cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base T (94/42);  new base C (3/7);  total (98/50)
Fisher's exact test for biased strand distribution p-value = 3.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

GCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGC  >  NC_002947/4863009‑4863298
                                                                                                                                                  |                                                                                                                                                  
 cggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGt                                                                                                                                              <  1:5245636/150‑1 (MQ=255)
   gTCATGGGATTTGTGTTCTCTTGCATTGCAATAGAGCAGCTTGGGGTGGTGGGTGGGGTTGCAGGTGTCTGTTGACTGTGCCAGCCTCTTGGCGGGTTTAGCCGCTGCGAGAAGGGGGTCAGTTTGCCTTGGGCGGTGTTGTTGCTGttg                                                                                                                                            >  2:2960619/1‑148 (MQ=255)
   gTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAg                                                                                                                                            <  1:5023431/150‑1 (MQ=255)
      aTGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCTACCACAACGGCGTCACATTCCCAAGGGCGGTGATGTACCAGTAGGAg                                                                                                                                         >  1:3517425/1‑150 (MQ=255)
       tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGc                                                                                                                                        >  1:2515065/1‑150 (MQ=255)
       tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGc                                                                                                                                        <  2:3610661/150‑1 (MQ=255)
        gggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTGCCACAAGGGCGTGACATGCCCTAGGGGGGTGATGTAGGTGTTGGAGCg                                                                                                                                      >  2:4632136/1‑150 (MQ=255)
        gggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCg                                                                                                                                      <  1:1472600/150‑1 (MQ=255)
         ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                                     <  2:2515065/150‑1 (MQ=255)
          gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCggg                                                                                                                                    <  1:3916244/150‑1 (MQ=14)
          gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                                    <  2:3848958/150‑1 (MQ=255)
          gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                                    <  2:4856083/150‑1 (MQ=255)
          gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCGCCCACAACGGCCTCACAAGCCCAAGGGCGGGGATGTACCGGTAGGAGCggg                                                                                                                                    >  1:3294550/1‑150 (MQ=255)
           aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                                                   <  1:138900/150‑1 (MQ=14)
           aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                                                   <  2:1807067/150‑1 (MQ=14)
           aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                                                   <  1:3403243/150‑1 (MQ=14)
           aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGt                                                                                                                                   <  2:690661/150‑1 (MQ=255)
             ttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGttt                                                                                                                                 <  2:1132315/150‑1 (MQ=255)
                tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTAcc                                                                                                                              >  1:4464072/1‑150 (MQ=14)
                tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTaac                                                                                                                              >  1:4489415/1‑148 (MQ=255)
                 cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTAccc                                                                                                                             <  2:1879393/150‑1 (MQ=14)
                 cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCGCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCGGGAGGAGCCGGTTTAccc                                                                                                                             >  2:938547/1‑150 (MQ=255)
                  ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACACGCGCACCCCACGGCGGCACAAGCCCACGGGCCGTGATTTAACCGTTGGATCCGGTTTAaccc                                                                                                                            >  1:1105550/1‑149 (MQ=255)
                   tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCcgc                                                                                                                           >  2:4950666/1‑150 (MQ=255)
                   tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCACCCACACGGGCGTCACATGCCCGAGGGCGGTGATGTAACAGTAGGAGCGGGGTTAccccc                                                                                                                           >  1:4596654/1‑148 (MQ=255)
                     tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGTATACCCCGCCCACAAGGGCGGGACAAGCCCGAGGGCGGGGGTGTGCCGGGTGGTGCGGGTGTTACCGCGa                                                                                                                         >  1:636720/1‑150 (MQ=255)
                     tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTAACCCGCGCCCCCAAGGGCGGCACATGCCCGAGGGCGGGGGTGTGCCGGGAGGAGCGGGGTTTAccccgc                                                                                                                          >  1:2476362/1‑147 (MQ=255)
                      ctTGCAAAGCAAAAGAGCAGCTGGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAGGCCCTAGGGCGGTGTTGTACCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                        <  2:3961019/150‑1 (MQ=255)
                       ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATTCCCTAGGGCGGTGAAGTACCTGTAGGAGCGGGTTTACCCGCGAag                                                                                                                       >  2:1119677/1‑150 (MQ=255)
                       ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAag                                                                                                                       <  2:3242638/150‑1 (MQ=255)
                       ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCGGTAGGAGCGGGTTTACCCGCGAag                                                                                                                       >  2:5344297/1‑150 (MQ=255)
                         gctttgcaagAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                                     >  2:2776503/6‑150 (MQ=255)
                          caaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCTGTCGGAGCGGGTTTACCCGCGAAGCgg                                                                                                                    >  2:3728179/1‑150 (MQ=255)
                           aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTAACCCCCGCCCACAACGGCGGCACAAGCCCAAAGGCGGGGATGTACCGGTAGGAGCGGGGTTAACCCCGAAGAGGc                                                                                                                   >  1:2834589/1‑150 (MQ=255)
                           aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGGAACCCCCTCCCCCAAGGGCGGGACAAGCCCAACGGCGGGGGTGTAACGGTAGGAGCCGGTTTAACCCCCAAGAGGc                                                                                                                   >  1:994536/1‑150 (MQ=255)
                              gcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCTGGTGTCTGTTGACTGGGCCAGCCTCTTCGCGGTTAAACCCGCTCCCACAAGGGCGTCACATGCCCTCGGGCGGTGATGTACCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCgg                                                                                                                <  2:4826177/150‑1 (MQ=38)
                              gcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCgg                                                                                                                >  2:809346/1‑150 (MQ=255)
                               caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTATGGCGGTGATGTACCGGTAGGAGCGGGTTTAACCGCGAAGAGGCCGGc                                                                                                               >  2:2924829/1‑150 (MQ=255)
                               caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGAAACCCCCGCCCCCAACGGCGTCAAAAGCCCAAGGGCGGGGATGTACCTGGAGGAGCGGGGTTTAACCCGAAGAGGCCCGc                                                                                                               >  1:4688467/1‑150 (MQ=255)
                                aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCTAGGGCGGTGATGTACCTGGAGGAGCGGGTTTACCCGCGAAGAGGCCGGca                                                                                                              >  2:4278432/1‑150 (MQ=255)
                                 aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcac                                                                                                             >  1:4432688/1‑150 (MQ=255)
                                 aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCGCCCACAAGGGCGTCACAAGCCCTAGGGCGGGGCTGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcac                                                                                                             >  1:4432581/1‑150 (MQ=255)
                                 aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCGCCCACAAGGGCGTGCCCAGCCCGTGGGGGGTGGTGTACCTGTAGGAGGGGGTTTTCCCCCGAAGAGGCCGCcac                                                                                                             >  1:4753578/1‑150 (MQ=255)
                                 aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCGCCCACAAGGGCGTCACAAGCCCAAGGGCGGTGCTGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcac                                                                                                             >  1:3946110/1‑150 (MQ=38)
                                  aaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCTGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcaca                                                                                                            >  2:5013574/1‑150 (MQ=255)
                                     agcagcTGGGGGTGGTGCGTGGGGTGCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGCGGGGGTTTACCCGCGAAGAGGCCGGCACAGGc                                                                                                         <  2:3932986/150‑1 (MQ=255)
                                                gtggtgCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATaa                                                                                              >  2:1112186/1‑149 (MQ=255)
                                                            gTTCCAGGTGTCTGTTGACTGGGCCAGCCTCTGCGCGGGTAACCCCGCTCCCACACGGGCGTCACAGGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCa                                                                                  <  2:5353530/150‑1 (MQ=38)
                                                             ttCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAg                                                                                 >  2:380411/1‑150 (MQ=255)
                                                             ttCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGGAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAg                                                                                 >  2:380400/1‑150 (MQ=255)
                                                                cAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCACCCACAAGGGCCTCACAAGCCCTAGGGCCGTGATGTACCGGTAGGAGCGGGTTTACCCCCGAAGAGGCCGGCACAGGCAACCACGATACCCGAACccgcccgccc                                                                              >  2:3408642/1‑148 (MQ=37)
                                                                       gtGTTGAGTGGGCCAGCCTCTTCGGGGTTAAACCCGCGCCCACCCGGGCGCCACTGGCCCTCGGGGGTTGTTGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCa                                                                       <  2:5182292/149‑1 (MQ=25)
                                                                       gtGTTCACTGGGCCCGCCTCTGCGGGGTTAAACCCGCCCCCACAGGGGCTACACATGCCCTCGGGGGGTTTTGTCCGTGTTGGCGGGGTTTTCCCCGCCAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCa                                                                       <  2:3789027/149‑1 (MQ=255)
                                                                        tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCCGCACGCCCTCAg                                                                      >  2:3586521/1‑150 (MQ=255)
                                                                        tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAg                                                                      >  2:1262574/1‑150 (MQ=255)
                                                                        tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCGCCCACAAGGGCGGCACATGCCCGAGGGCGGGGGTGTACCGGGAGGAGCGGGGTTACCCGCGCAGAGGCCGGCCCAGGCCACGCCGCAACCCGAACCGGCCAGCACGCCCGCAg                                                                      >  2:1098691/1‑150 (MQ=25)
                                                                        tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCACCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGCACAGGCCCGCCCCGGCAACGACGATACCCGAACCGGCCAGCACGCCCACAg                                                                      >  2:4707107/1‑150 (MQ=38)
                                                                          ttGACTGTCCCCCCCTCTTCGCGGTTAAACCCGCACCCAAACGGGGGAAACTTCCCCTTGGGGGTTTTTTTACCTGTTGGAGGGTTTTAACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGgc                                                                    <  2:2041170/150‑1 (MQ=255)
                                                                                   ccAGCCTCTTCGCGGGTAAACCCCCGCCCACAACGGCCTCACATGCCCCAGGGCGGGGATGTACCGGGAGGCGCCGGGGTACCCGCGAAGAGGCCCGCCCCCGCAACCCCCAACACCGCACCGGCCCGCCCCCCCACAggcggcgcgtgc                                                           >  2:5138286/1‑144 (MQ=255)
                                                                                         tctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACCATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTgg                                                     >  1:2229900/1‑150 (MQ=255)
                                                                                         tctGCGCGGGAACACCCGCTCCCACACGGGCGTCACATCCCCTTGGGGGGTGTGGTACGTGTTGGCGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTgg                                                     <  2:2834589/150‑1 (MQ=37)
                                                                                           ttCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTATGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTAGTCGGGAt                                                   >  2:2490477/1‑150 (MQ=255)
                                                                                              gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTg                                                >  1:1140630/1‑150 (MQ=255)
                                                                                              gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCGGGATCTg                                                >  1:4183843/1‑150 (MQ=255)
                                                                                                  ggaacccccGCTCCCACCCGGGCGCCACGGGCCCGGGGGGGGTGTTGCCCGTGTCGGCGGGGTTTCCCCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                                             <  2:656113/144‑1 (MQ=255)
                                                                                                  ggTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                                             <  2:1770828/150‑1 (MQ=255)
                                                                                                   gTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAg                                            >  2:1846280/1‑150 (MQ=255)
                                                                                                   gTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGAAACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAg                                            >  1:1758815/1‑150 (MQ=255)
                                                                                                    tAAACCCGCTCCCACCCGGGCGTCACCTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGAGACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgc                                           <  2:5361721/150‑1 (MQ=255)
                                                                                                       aCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTGTCGTCTGGATCTGTCAGCGCt                                        >  1:52675/1‑150 (MQ=14)
                                                                                                       aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                        >  1:284984/1‑150 (MQ=255)
                                                                                                        cccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCtt                                       >  2:3098937/1‑150 (MQ=255)
                                                                                                         ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                      >  1:1210814/1‑150 (MQ=255)
                                                                                                         ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                      >  1:1099513/1‑150 (MQ=255)
                                                                                                         ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                      >  1:3082241/1‑150 (MQ=255)
                                                                                                           gCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTg                                    >  1:1256361/1‑150 (MQ=255)
                                                                                                            cgcccACCAGGGCGTCACCTGCCCTTGGGGGGTGATGTACGTGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTgg                                   <  2:4308116/148‑1 (MQ=38)
                                                                                                             tCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg                                  >  2:5280309/1‑150 (MQ=255)
                                                                                                             tCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg                                  >  1:2696610/1‑150 (MQ=255)
                                                                                                             tCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg                                  >  2:2160229/1‑150 (MQ=255)
                                                                                                                 acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg                              >  1:509860/1‑150 (MQ=255)
                                                                                                                  caAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgc                             <  1:3098937/150‑1 (MQ=255)
                                                                                                                       gCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACg                        <  2:2182762/150‑1 (MQ=255)
                                                                                                                          tCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGgcg                     >  1:2739465/1‑150 (MQ=255)
                                                                                                                           cactggCCCTCGGGCGGTGTTTTACCTTTAGGCGCGGGTTAACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGgcgc                    <  2:2313100/145‑1 (MQ=38)
                                                                                                                           cactTGCCCTCGGGGGGTGTTTTCCGTGTCGGCGGGGTTTTCCCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGgcgc                    <  2:2313096/150‑1 (MQ=255)
                                                                                                                             ccggCCCTTGGGGGGTGTTGTACGTTTTGGCGGGGGTTAACCCGCTAAGTGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTc                  <  1:3502077/147‑1 (MQ=255)
                                                                                                                             caTGCCCTGGGGCGGTGATGTACCTGTCGGTGCGGGTTTCCCCGCGACGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTc                  <  1:453055/150‑1 (MQ=255)
                                                                                                                             caTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTc                  >  2:4369572/1‑150 (MQ=255)
                                                                                                                                 cccTTGGGCGGTGTTGTACCTGTAGGAGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAAc              <  2:1638807/150‑1 (MQ=255)
                                                                                                                                  ccTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACt             >  1:605997/1‑150 (MQ=255)
                                                                                                                                   cTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTc            >  2:3628532/1‑150 (MQ=255)
                                                                                                                                   cTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTc            >  2:1543492/1‑150 (MQ=255)
                                                                                                                                   cTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTc            >  2:3628639/1‑150 (MQ=255)
                                                                                                                                     aGGGCGGTGTTGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcga          <  2:447695/150‑1 (MQ=255)
                                                                                                                                      gggCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcgac         >  2:4221583/1‑150 (MQ=255)
                                                                                                                                           ggTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATa     >  2:3120505/1‑150 (MQ=255)
                                                                                                                                             tGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCt                                                       >  2:472814/1‑98 (MQ=255)
                                                                                                                                             tGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCt                                                       <  1:472814/98‑1 (MQ=255)
                                                                                                                                              gATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGc  <  1:4221583/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGC  >  NC_002947/4863009‑4863298

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