Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,115 C→G 8.0% intergenic (+85/‑203) PP_0367 → / → PP_0368 hypothetical protein/acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947445,1150CG 8.0% 438.2 / 2.3 157intergenic (+85/‑203)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (75/64);  new base G (8/7);  total (86/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

TCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGAT  >  NC_002947/444966‑445263
                                                                                                                                                     |                                                                                                                                                    
 cGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                                                     <  2:3917913/150‑1 (MQ=255)
   gcgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTga                                                                                                                                                   <  1:953567/150‑2 (MQ=255)
     gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTc                                                                                                                                                 <  1:5745383/150‑1 (MQ=255)
       aaaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                                               <  2:2979074/150‑1 (MQ=255)
       aaaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                                               <  1:834112/150‑1 (MQ=255)
         aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                             <  2:3597597/150‑1 (MQ=255)
         aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                             <  1:1621827/150‑1 (MQ=255)
         aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGCCCGACCAGACAAAATTGTTGTCTGATCggac                                                                                                                                             <  2:6248864/150‑3 (MQ=37)
         aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGCCCGACCAGACAAAATTGTTGTCTGATCggac                                                                                                                                             <  1:5819972/150‑3 (MQ=37)
                 aaGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgc                                                                                                                                     <  1:6328727/150‑1 (MQ=255)
                  aGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgcg                                                                                                                                    <  1:5579661/150‑1 (MQ=255)
                   ggCATCCGACGACAGGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCGGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                   <  2:5728282/150‑1 (MQ=17)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCGTCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                   <  2:6197076/150‑1 (MQ=18)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGGGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                   <  1:6353561/150‑1 (MQ=18)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                   <  1:2231502/150‑1 (MQ=255)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                   <  1:6089021/150‑1 (MQ=255)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                   <  2:4642513/150‑1 (MQ=255)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                   <  2:2212471/150‑1 (MQ=255)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                   <  1:3433796/150‑1 (MQ=255)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                   <  2:3344718/150‑1 (MQ=255)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                   <  1:987461/150‑1 (MQ=255)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                   <  2:6308824/150‑1 (MQ=21)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                   <  1:3887943/150‑1 (MQ=18)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                   <  1:6434391/150‑1 (MQ=18)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                   <  2:3327848/150‑1 (MQ=18)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                   <  1:2265969/150‑1 (MQ=18)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                   <  2:6789821/150‑1 (MQ=18)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                   <  1:4899146/150‑1 (MQ=18)
                   ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                   <  2:5280060/150‑1 (MQ=18)
                   ggCATCCGACGAAATCCTTGGCCGGCGCATCGACACCTCCGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                   <  2:4664967/150‑1 (MQ=17)
                   ggCATCCGACGAAATACTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                   <  2:4665078/150‑1 (MQ=18)
                         cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGATCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                                             >  2:4314774/1‑150 (MQ=255)
                         cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGc                                                                                                                             <  1:6597807/150‑1 (MQ=17)
                          gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGcc                                                                                                                            <  1:1272922/150‑1 (MQ=17)
                           acgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTTCCCCCCACAGCGGCCCATCAAACAACAAATTTGTGTGGTGGGGCCGCGTGGCGGGCCAcccc                                                                                                                           >  2:4128693/1‑150 (MQ=255)
                             gaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTTCCCCCGCAAGCGGGCCGTCCGACAACCAATTTGTGTGGGGGGGGCGCGTAGCCGGCCAGccccc                                                                                                                         >  2:2416704/1‑148 (MQ=255)
                             gaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGc                                                                                                                         <  1:3342489/150‑1 (MQ=255)
                               aaTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGTTTCCCCCGAATGCGTCCGATAAGAAAACAATTTTGTGTGGTAGGGCCGCTTAGCGGGCAAGcccccgc                                                                                                                       >  2:5062924/1‑146 (MQ=255)
                               aaTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCGCGAATGCGTCCGATAAGACAACAAATTTGTCTGGTCGGGCCTCGTCGCGGGCAAGCCCGcgc                                                                                                                       >  2:5614083/1‑148 (MQ=255)
                                           gATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCCGCGCGTGCCTTGTCGGTGGGGGCTTGCCCGCGAATGCGTCCGATCCGACAACAATTTTGTCTGGTGGGGCCTCTTCGGGGGCAAGCCCGCTCCTACAAGGGttt                                                                                                           <  2:1639732/150‑1 (MQ=38)
                                             tcatcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCCGCGCCGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGAGGGCAAGCCCGCTCCTACAAGGGTTtgt                                                                                                         <  1:6475090/150‑1 (MQ=255)
                                                   ccacCTCCGCCTTAATCCTTGCGTCACCCCCTATTCCCGCGCCGGCCTTGTCGGAGCGGGCTGGCCCGCGAATGCGCCCTACCAGAAAACAATTTTGCCGGGTCGGCCCACTTCGCGGGAACCCCCGCTCCTACAAGGGTTTGtggcatg                                                                                                   <  2:1043694/149‑1 (MQ=255)
                                                       cTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCCCGAATTCGTGCGATAAGACAACAAATTTTTGTGGTCGGGGCTCTTCGCGGGCAAGCCCGCACCAACAAGGGTTTGTggcatggcat                                                                                               >  2:6154622/1‑150 (MQ=255)
                                                        tACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTTCCCGCGAAATCCTCCCCACAGACAACAAATTTTTGTGGGGGGGGCGCGTGGGGGGGAAGCCCCCGCCAAAAACGGGTTGTGggaatgcatg                                                                                              >  2:1551938/1‑150 (MQ=255)
                                                        cccGCCTTATTCCTGGCGTCAGCCCCTATTCCCGCGCCGGCCTTGTTGGTGGGGGGGTGCCCGCTAATGCGTCCCATCACACAACAATTTTGTGTGGTGGGGCCTCTTGGCGGGAACGCCCGCTCCTACAAGGGTTTGTGgcatggcatg                                                                                              <  2:915010/148‑1 (MQ=255)
                                                               taaTCCTTGCTCCACCCCCTTTTCCCGCGCGTGCTTTGTTGGCGGGGGCTCGCCGTCTAATGCGCCCTAAAAAAAAAAATTTTTGTGTGGGCGGCCCTCTTCGCGGGAACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAtt                                                                                       <  2:1971625/149‑1 (MQ=255)
                                                                aaTACTTGCTTCTTCCCCTTTTTCCGCGCCTGCCTTGTAGGTGCGGGCTGGCCCGCGATTGCGTCCGATCAGACAACAATTTTGTCTGGTGGGGCGGCTTCGCGGGCAAGCCAGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                      <  2:6548226/150‑1 (MQ=255)
                                                                aaTACTTGCGTCACCCCCAATTCCCGCGCCGGCCTTGTTGGCGGGGGCGGGCCCGCGAATGCGCCCGAACACAAAACATTTTTGGGGGGGGGGCCCTCTCCGGGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                      <  2:2495574/150‑1 (MQ=255)
                                                                               ccttttttCCGCGCCTGCCTGTTTGGAGGGGGTTGGCCGGCTACGGCTCCCGATAACCAAACATTTTTGTGGGGTGGGGCATCTTCGGGGGAAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                                       <  2:2980893/145‑1 (MQ=255)
                                                                               cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATCTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                                             <  1:1901925/96‑1 (MQ=255)
                                                                               cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATATTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                                             >  2:1901925/1‑96 (MQ=255)
                                                                               cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCGCGAATGCGTCCCATCAGACAACAAGTTTTTCTGGTCGGGCCACTTAGCGGGCAAGCCCCCACCAACAAGGGTTTGTGGCATGGCAATACTGATTTTCCAATTTTCGTGGc                                                                       >  2:6155046/1‑150 (MQ=255)
                                                                               cccTATTCCCGCCCCTGCCTGGGCGCAGGGGGGTGGCCCGCTAATGCGTCCTAACCGACAAAAATTTTGTCTGGTGGGGCCTCTCCGGGGGCTAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                                       <  2:1540046/150‑1 (MQ=255)
                                                                                ccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTGCCCCGCATTCGTCCGATCAGACAACAAATTTTTGAGGTCGGGCCACATAGCGGGCCAGCCCCCTCCAACAAGGCTTTGTTGGAAGGCCTGACTGATTTTACAATTTTAGTgggg                                                                      >  2:5521312/1‑148 (MQ=255)
                                                                                  tATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTTCCCCCCAATTCGTCCGATAACACAACAATGTTTTCTGGGCGGGCCCCATCGCCGGCAAGCCCCCGCCGACAAGGGGTTGTGGCATGGCCAGACTGATTTGCCAAATTTAGTGGCGTa                                                                    >  2:5125112/1‑150 (MQ=255)
                                                                                   aTTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                                               <  2:1037009/74‑1 (MQ=255)
                                                                                   aTTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                                               >  1:1037009/1‑74 (MQ=255)
                                                                                                     ttgttgggggTTTGCCCGCTATTTCGTCCTTTTTTACAAAATTTTTGTGTGGTGGGGCTTCTGGGGGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATc                                                 <  2:3039172/143‑1 (MQ=255)
                                                                                                       ggAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTGATc                                                                                                                                                 >  2:1673956/1‑52 (MQ=11)
                                                                                                       ggAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTGATc                                                                                                                                                 <  1:1673956/52‑1 (MQ=11)
                                                                                                         aGCGGGCTTGCCCGCGACTGGGTCCGATCAGACAACAATTTGGTGTGGTCGGGCCTCTTCGGGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTg                                             <  2:4414413/150‑1 (MQ=255)
                                                                                                                ttGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                      >  1:2404310/1‑150 (MQ=17)
                                                                                                                ttGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                      >  1:2404366/1‑150 (MQ=17)
                                                                                                                ttGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                      >  1:5089658/1‑150 (MQ=17)
                                                                                                                  gCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc                                    >  2:4254679/1‑150 (MQ=255)
                                                                                                                   cccGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                   >  2:3285808/1‑150 (MQ=18)
                                                                                                                   cccGCCATTGCGCCCTAACACACAACAATTTGGTCTGGGCGGCCCTCTGCGCGGCCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                   <  2:104155/150‑1 (MQ=255)
                                                                                                                    ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                  >  1:2768266/1‑150 (MQ=255)
                                                                                                                    ccGCGAATGCGTCCGATCAGACAACAATTATGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                  >  2:458198/1‑150 (MQ=255)
                                                                                                                           gaggcccGACCAGACAAAATTGTTGTCTGATCGGACCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTaaa                           >  2:6301695/6‑150 (MQ=37)
                                                                                                                              gcccGACCAGACAAAATTGTTGTATGATCGGACCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGAGTAAATTGAACCATATTACATCCTTGACAGAAAGCTGAGCTAAACGt                        >  2:5896096/3‑150 (MQ=25)
                                                                                                                               tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa                       >  1:1822345/1‑150 (MQ=255)
                                                                                                                                ccGATCAGACAACAATTTTGTCTGGTCGGGCCTCTGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa                      <  2:6799500/150‑1 (MQ=255)
                                                                                                                                ccGACCAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa                      >  1:4927402/1‑150 (MQ=255)
                                                                                                                                ccGACCAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa                      >  2:6045023/1‑150 (MQ=255)
                                                                                                                                ccGACCAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa                      >  1:4927417/1‑150 (MQ=255)
                                                                                                                                  gaccAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGt                    >  1:3714598/4‑150 (MQ=255)
                                                                                                                                   aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt                   >  2:887714/1‑150 (MQ=255)
                                                                                                                                     cAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc                 >  1:186560/1‑150 (MQ=255)
                                                                                                                                       gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa               <  2:4901911/150‑1 (MQ=255)
                                                                                                                                       gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTGTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa               >  2:6222994/1‑150 (MQ=255)
                                                                                                                                        acaacaATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa              >  1:5512340/1‑150 (MQ=255)
                                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGTCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca            >  2:2762821/1‑150 (MQ=255)
                                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca            >  1:681356/1‑150 (MQ=255)
                                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca            >  1:3850816/1‑150 (MQ=255)
                                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTGGTGGCATGGCATGACGGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCGTGTCAGAAAGCTGAGCTAGTCGTCAGTTTCattca            >  1:2824678/1‑146 (MQ=255)
                                                                                                                                           acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa           >  1:2155743/1‑150 (MQ=255)
                                                                                                                                            caaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac          >  2:855248/1‑150 (MQ=255)
                                                                                                                                            aattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac          >  1:5386736/6‑150 (MQ=38)
                                                                                                                                             aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt         >  1:1013076/1‑150 (MQ=255)
                                                                                                                                             aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt         >  1:1011017/1‑150 (MQ=255)
                                                                                                                                             aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt         >  1:3859012/1‑150 (MQ=255)
                                                                                                                                               ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGt       >  1:5652764/1‑150 (MQ=255)
                                                                                                                                               ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGt       >  2:2892421/1‑150 (MQ=255)
                                                                                                                                               ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGt       >  2:3420601/1‑150 (MQ=255)
                                                                                                                                                  tGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTg    <  1:741949/150‑1 (MQ=255)
                                                                                                                                                   gTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGa   >  1:324196/1‑150 (MQ=255)
                                                                                                                                                    tCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGAt  >  1:6809531/1‑150 (MQ=255)
                                                                                                                                                    tCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGAt  >  2:217932/1‑150 (MQ=255)
                                                                                                                                                    tCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGAt  >  2:5272247/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGAT  >  NC_002947/444966‑445263

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