Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,863,154 A→C 7.5% intergenic (‑138/+99) 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,1540AC 7.5% 405.4 / 16.8 148intergenic (‑138/+99)ligA/PP_4275DNA ligase/cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base A (92/41);  new base C (3/9);  total (98/50)
Fisher's exact test for biased strand distribution p-value = 3.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 149 total aligned reads displayed.

GCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGCTGGCG  >  NC_002947/4863009‑4863303
                                                                                                                                                   |                                                                                                                                                      
   gTCATGGGATTTGTGTTCTCTTGCATTGCAATAGAGCAGCTTGGGGTGGTGGGTGGGGTTGCAGGTGTCTGTTGACTGTGCCAGCCTCTTGGCGGGTTTAGCCGCTGCGAGAAGGGGGTCAGTTTGCCTTGGGCGGTGTTGTTGCTGttg                                                                                                                                                 >  2:2960619/1‑148 (MQ=255)
   gTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAg                                                                                                                                                 <  1:5170034/150‑1 (MQ=255)
   gTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAg                                                                                                                                                 >  2:2960647/1‑150 (MQ=255)
      aTGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCTACCACAACGGCGTCACATTCCCAAGGGCGGTGATGTACCAGTAGGAg                                                                                                                                              >  1:3517425/1‑150 (MQ=255)
       tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGc                                                                                                                                             <  2:3845764/150‑1 (MQ=255)
       tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGc                                                                                                                                             <  2:3610699/150‑1 (MQ=255)
       tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGc                                                                                                                                             <  2:3610661/150‑1 (MQ=255)
       tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGc                                                                                                                                             >  1:2515065/1‑150 (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:4856083/150‑1 (MQ=255)
          gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCGCCCACAACGGCCTCACAAGCCCAAGGGCGGGGATGTACCGGTAGGAGCggg                                                                                                                                         >  1:3294550/1‑150 (MQ=255)
           aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                                                        <  1:3403243/150‑1 (MQ=14)
           aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                                                        <  1:138900/150‑1 (MQ=14)
           aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                                                        <  2:1807067/150‑1 (MQ=14)
           aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGt                                                                                                                                        <  2:3221918/150‑1 (MQ=255)
           aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGt                                                                                                                                        <  2:690661/150‑1 (MQ=255)
                tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTaac                                                                                                                                   >  1:4489415/1‑148 (MQ=255)
                 cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCGCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCGGGAGGAGCCGGTTTAccc                                                                                                                                  >  2:938547/1‑150 (MQ=255)
                  ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACACGCGCACCCCACGGCGGCACAAGCCCACGGGCCGTGATTTAACCGTTGGATCCGGTTTAaccc                                                                                                                                 >  1:1105550/1‑149 (MQ=255)
                   tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCcgc                                                                                                                                >  2:4950666/1‑150 (MQ=255)
                   tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCGAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCcgc                                                                                                                                >  2:1843698/1‑150 (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:2468055/1‑150 (MQ=255)
                       ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCGGTAGGAGCGGGTTTACCCGCGAag                                                                                                                            >  2:5344297/1‑150 (MQ=255)
                         gcaaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCCCGAagag                                                                                                                          >  1:4482330/1‑150 (MQ=255)
                           aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTAACCCCCGCCCACAACGGCGGCACAAGCCCAAAGGCGGGGATGTACCGGTAGGAGCGGGGTTAACCCCGAAGAGGc                                                                                                                        >  1:2834589/1‑150 (MQ=255)
                              gcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCTGGTGTCTGTTGACTGGGCCAGCCTCTTCGCGGTTAAACCCGCTCCCACAAGGGCGTCACATGCCCTCGGGCGGTGATGTACCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCgg                                                                                                                     <  2:4826177/150‑1 (MQ=38)
                              gcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCgg                                                                                                                     >  2:809346/1‑150 (MQ=255)
                               caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGAAACCCCCGCCCCCAACGGCGTCAAAAGCCCAAGGGCGGGGATGTACCTGGAGGAGCGGGGTTTAACCCGAAGAGGCCCGc                                                                                                                    >  1:4688467/1‑150 (MQ=255)
                                naaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCCCACGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGca                                                                                                                   >  1:5380945/2‑150 (MQ=255)
                                aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTGAACCCCCCCCCCCCAGGGGGGCACAAGCCCCCGGGCGGGGGTGGACCCGTGGGGGCGGGGGTGACCGCCAAGAGggcggcc                                                                                                                   >  1:656113/1‑149 (MQ=255)
                                aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGca                                                                                                                   >  2:2071994/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)
                                          gtgggggTGGTGGGGGGGTTGCCTGGTGTCTGTTGAGTGGGCCCGCCTCTTCGGGGTTAAACCCGCCCCAACCGGGGCGACACTGGCCCTGGGGGGGTTTTGTACGTTTGGGCGGGGTTTTCCCCGCGAAGAGGCCGGCACAGGCAacga                                                                                                         <  2:1105550/147‑1 (MQ=255)
                                             ggggTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAcc                                                                                                                                   <  1:4939437/121‑1 (MQ=255)
                                                  gggggggggggTCCCTGGTTTTTGTTGACTGTGCCAGCATGTCCGGGGTTAAACCCGCGCCCACCCGGGGGTAACTGCCCCTTGGGGGTTTTTTAACGTTTTGGAGCGTTTTTACCCGCGAAGAGGCCGGAACAGGCAACGACGATAccc                                                                                                 <  2:2518033/143‑1 (MQ=255)
                                                            gTTCCAGGTGTCTGTTGACTGGGCCAGCCTCTGCGCGGGTAACCCCGCTCCCACACGGGCGTCACAGGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCa                                                                                       <  2:5353530/150‑1 (MQ=38)
                                                             ttCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAg                                                                                      >  2:380257/1‑150 (MQ=255)
                                                             ttCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAg                                                                                      >  2:380411/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)
                                                                        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)
                                                                                   cccgCCTTTTCGCGTTTAACCCCGCTCCCCAAGGGGCGTCCCTGGCCCTCGGGCGTGGTTGTACCTGTCGGAGGGGTTTTACCGGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTc                                                                <  1:4532465/147‑1 (MQ=255)
                                                                                   ccAGCCTCTTCGCGGGTAAACCCCCGCCCACAACGGCCTCACATGCCCCAGGGCGGGGATGTACCGGGAGGCGCCGGGGTACCCGCGAAGAGGCCCGCCCCCGCAACCCCCAACACCGCACCGGCCCGCCCCCCCACAggcggcgcgtgc                                                                >  2:5138286/1‑144 (MQ=255)
                                                                                         tctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATTCCCGAGGGCGGTGTTGTACCTGTAGGAGCGGGTTTACCCCCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCACAGGCTGCTCGTTCGTCTgg                                                          >  2:1492684/1‑150 (MQ=38)
                                                                                         tctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTgg                                                          >  2:2213076/1‑150 (MQ=255)
                                                                                         tctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTgg                                                          >  2:3631964/1‑150 (MQ=255)
                                                                                         tctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACCATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTgg                                                          >  1:2229900/1‑150 (MQ=255)
                                                                                         tctGCGCGGGAACACCCGCTCCCACACGGGCGTCACATCCCCTTGGGGGGTGTGGTACGTGTTGGCGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTgg                                                          <  2:2834589/150‑1 (MQ=37)
                                                                                              gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTg                                                     >  1:1140630/1‑150 (MQ=255)
                                                                                              gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCGGGATCTg                                                     >  1:4183843/1‑150 (MQ=255)
                                                                                                 gggTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTc                                                   >  2:1668817/1‑150 (MQ=255)
                                                                                                  ggaacccccGCTCCCACCCGGGCGCCACGGGCCCGGGGGGGGTGTTGCCCGTGTCGGCGGGGTTTCCCCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                                                  <  2:656113/144‑1 (MQ=255)
                                                                                                  ggTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                                                  >  2:34345/1‑150 (MQ=255)
                                                                                                   gTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAg                                                 >  2:1846280/1‑150 (MQ=255)
                                                                                                   gTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGAAACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAg                                                 >  1:1758815/1‑150 (MQ=255)
                                                                                                    tAAACCCGCTCCCACCCGGGCGTCACCTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGAGACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgc                                                <  2:5361721/150‑1 (MQ=255)
                                                                                                       aCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCTGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                             >  2:1119734/1‑150 (MQ=14)
                                                                                                       aCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTGTCGTCTGGATCTGTCAGCGCt                                             >  1:52675/1‑150 (MQ=14)
                                                                                                       aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                             >  1:284984/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                                       >  1:2696610/1‑150 (MQ=255)
                                                                                                             tCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg                                       >  2:5280309/1‑150 (MQ=255)
                                                                                                                 acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg                                   >  1:509860/1‑150 (MQ=255)
                                                                                                                  caAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgc                                  <  1:3098937/150‑1 (MQ=255)
                                                                                                                       gCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACg                             <  2:284984/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)
                                                                                                                            acaTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCt                        >  1:472182/1‑150 (MQ=255)
                                                                                                                             ccggCCCTTGGGGGGTGTTGTACGTTTTGGCGGGGGTTAACCCGCTAAGTGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTc                       <  1:3502077/147‑1 (MQ=255)
                                                                                                                             caTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTc                       >  2:4369572/1‑150 (MQ=255)
                                                                                                                                 cccTTGGGCGGTGTTGTACCTGTAGGAGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAAc                   <  2:1638807/150‑1 (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)
                                                                                                                                              gATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGc       <  1:4221583/150‑1 (MQ=255)
                                                                                                                                                   aCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGCTGGCg  >  1:4575183/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGCTGGCG  >  NC_002947/4863009‑4863303

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: