Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,125 C→G 13.4% intergenic (+95/‑193) 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,1250CG13.4% 475.6 / 59.8 172intergenic (+95/‑193)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (92/57);  new base G (10/13);  total (102/70)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-01

Only 100 of 177 total aligned reads displayed.

GTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCGCC  >  NC_002947/444977‑445274
                                                                                                                                                     |                                                                                                                                                     
gtcgtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCct                                                                                                                                                      <  2:2989247/150‑1 (MQ=255)
   gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCtt                                                                                                                                                   <  2:7671043/150‑1 (MQ=255)
    tcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTGCCCCCGAAAGCGTCCGGGCGGACAACAAATATTTGTGGTGGGGccgcgtc                                                                                                                                                  >  1:6800816/1‑145 (MQ=255)
        tgCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                              <  2:8511761/149‑1 (MQ=255)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCTCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgcgt                                                                                                                                              <  2:2562934/150‑2 (MQ=255)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                              <  2:3904066/150‑1 (MQ=255)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                              <  2:3465918/150‑1 (MQ=255)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                              <  2:3351809/150‑1 (MQ=255)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                              <  1:7031319/150‑1 (MQ=18)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                              <  2:3488369/150‑1 (MQ=18)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                              <  2:7880979/150‑1 (MQ=18)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                              <  1:4546502/150‑1 (MQ=18)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                              <  1:9099582/150‑1 (MQ=18)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACATACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                              <  2:2055880/150‑1 (MQ=255)
              cgacgaAATGCTTGGCCGGATCTTCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                                                        <  1:41504/150‑1 (MQ=255)
              cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGc                                                                                                                                        <  1:3619271/150‑1 (MQ=17)
              cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGc                                                                                                                                        <  1:7711284/150‑1 (MQ=17)
                          tGGCCGGACCATCGACACCTACGCCTGAATCCTGGCTTCATCCCCTATTTCCGCGCCGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGACCGATCAGACAACAATTTTGTGTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                                            <  2:7632886/150‑1 (MQ=37)
                                  acatcGACACCTACGCCTCAATACTTGCGACATCCCCTATCCCCGCGCCGGCCTGGTCGGAGCGGGCTTGCCCGCGACTGCGCCCGATCACACAACAATTTTGTCTGGTGGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTtgt                                                                                                                    <  2:3768562/149‑1 (MQ=37)
                                   catcGACACCTACGCCTGAATACTTGCGTCATCCCCTCTTTCAGCGCCTGCCTTGTAGTCGCGGGCTTGCCCGCGAATGCGTCGGACCACACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGATtgtg                                                                                                                   <  2:8205472/150‑1 (MQ=38)
                                       gACACCACCGCGTTAATAGTTGGGACAGCCCCTATTCCCGCGCGGGCTTTGTAGGCGGGGGGTGGCCGGCTATGGGGCCCGATCACACAAAAATTTGGTGGGGTGGGGCCTCTGCGGGGCCACCCCCGCTCCTACAAGGGTTTGtggcat                                                                                                               <  2:8730419/150‑1 (MQ=255)
                                         cacCTCCGCGTTATTCCTGGCGTCAGCCCCTTTTCCCGCGCCTGCCTTGTCGGCGGGGGCTGGCCCGCTACTGCGCCCGAACACAAAAAAATTTGGTCTGGGCGGGCCACTGCGGGGGAACGCCCGCTCCTACAAGGGTTTGtggcatgg                                                                                                             <  2:8073700/150‑1 (MQ=255)
                                                 ccTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCCCCAATGCGTCCCCTCAGACAAAAAATTTTTGTTGTGGGGCCGCGTCGCGGGCAACCCCCCGCCAAAAAGGGGTTTTTGGAAGGCATTAcag                                                                                                     >  2:7209568/1‑148 (MQ=255)
                                                    gAATCCTTGCGTCATCCCCTATTTCCGCGCCTGCCTGGTAGGCGCGGGCTTGCCCGCGAATGCGTCCGACCAGACAACAATTTTGTCTGGTGGGGCCTCTGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAtt                                                                                                  <  2:6640805/150‑1 (MQ=38)
                                                                 aTCCCCTATTCCCGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGACTGCGTCCGATCAGACACCAATTTTGTGTGGTCGGGCCTCTGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGt                                                                                     <  2:8194288/150‑1 (MQ=38)
                                                                  cccccAATTCCCGCGCGGGCCTTGTAGGCGGGGGCTTGCCCGCGATTGCTGCCGAACCCCAAACATTGTTGTGTGGTCGGGCCTCTCCGGGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTg                                                                                    <  2:8689215/149‑1 (MQ=255)
                                                                    cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTGCCCCCCAATTCGGGCGCTCAGACAAAAAATTTTTGTGGTCGGGGCGCGTCCCGGGCCATCCCCCACCAACAATGGTTTTTGGCATGGCATGAAAGATTTTCCAATTTTGGTGGc                                                                                  >  2:1414537/1‑150 (MQ=255)
                                                                     ccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAATTTGTCTGGTCGGGCCTCGTAGCGGGCAATCCCGCTCCTACAAGGGTTTGTTGCATGGCATGACTGATTTGCCAATTTTCGTTGCg                                                                                 >  2:9155378/1‑150 (MQ=38)
                                                                     ccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCATTCCCGCGAAAGCCTACCATAAGACAACAATTATGTCTGGTCGGGCCGCTTCGCCGGCAAGACCCCACCGACAAGGGTTTTTGGCAAGGCATGACTGATTTTCCAATTTTCGTGGCg                                                                                 >  2:7524415/1‑150 (MQ=255)
                                                                         ttcccgcgcCTGCCTTTTCGGGGCGGGCGGCCCGGCTACTGCGCCCCCCCCCACACCAAGTGTGTCGGGTGGGCCCTCGACGGGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                                             <  2:1922313/145‑1 (MQ=255)
                                                                                   tgcttTGTCGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGAAAACATTTTGGTCTGGTCGGGCCTCTTTGCGGGAAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATa                                                                   <  1:2474849/150‑1 (MQ=255)
                                                                                   gcctttgttgGCGCGGGCGGCCCCGCGACTGCGTCCCAACCAACACCATTTTGGGGTGGTGGGGCCTCCCCGCGGGCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATa                                                                   <  1:7537879/146‑1 (MQ=255)
                                                                                           tggAGCGGGGTTGCCCGCTAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCATGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATcc                                                           <  2:6740798/149‑1 (MQ=255)
                                                                                             gAGCGGGCTTGCCCGCTAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTCCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                                         <  2:4564600/150‑1 (MQ=255)
                                                                                              gcgggggTTGCCCGGGTATTGCGTCCTATCACACAAAATTTTTGTGGGGCGGGGCCTCTGCGCGGCCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                                         <  1:2986595/147‑1 (MQ=255)
                                                                                               gCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTTTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTTCCATTTTTCGTGGCGTAAATTGCACCATAGTCCATCATTGa                                                       >  1:7998886/1‑150 (MQ=255)
                                                                                               gCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGGTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCCTAGTCCATCCTTgc                                                       >  2:5801179/1‑149 (MQ=255)
                                                                                               gCGGGCTTGCCCGCGAAAGCGTCCGAACAGACAACAATTTTGTCTGGGCGGGCCTCTACGCGGGCAAGCCCCCTCCTACAAGGGGTTGTGGCATGGCATGACTGATTTGCCAATTTTCGTGGCGTAAATGGCACCAAAGTCCATACTTGa                                                       >  1:8819618/1‑150 (MQ=255)
                                                                                                   gCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                                                   >  2:694331/1‑150 (MQ=255)
                                                                                                     ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                                 >  1:7281357/1‑150 (MQ=255)
                                                                                                        cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                              >  2:7230839/1‑150 (MQ=255)
                                                                                                        cccGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                              >  2:2994741/1‑150 (MQ=18)
                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGGGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                             >  1:4993116/1‑150 (MQ=255)
                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                                            >  2:7185886/1‑71 (MQ=255)
                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                                            <  1:7185886/71‑1 (MQ=255)
                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGTCGTAAAt                                                                            >  2:7172161/1‑119 (MQ=255)
                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGTCGTAAAt                                                                            <  1:7172161/119‑1 (MQ=255)
                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                             >  2:662023/1‑150 (MQ=255)
                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                             >  2:6175340/1‑150 (MQ=255)
                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                             >  1:1901012/1‑150 (MQ=255)
                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                             >  1:444634/1‑150 (MQ=255)
                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                             >  2:3094776/1‑150 (MQ=255)
                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                             >  1:6503236/1‑150 (MQ=255)
                                                                                                          cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                            >  2:1368129/1‑150 (MQ=255)
                                                                                                                    gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt                                   >  2:400243/1‑150 (MQ=255)
                                                                                                                     tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa                                  >  1:7493975/1‑150 (MQ=255)
                                                                                                                      ccGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                            <  1:1263837/139‑1 (MQ=255)
                                                                                                                      ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                                 >  2:7055818/1‑86 (MQ=18)
                                                                                                                      ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                                 <  1:7055673/86‑1 (MQ=18)
                                                                                                                       cGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAg                                >  2:1011676/1‑150 (MQ=255)
                                                                                                                         aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt                              >  1:6246004/1‑150 (MQ=255)
                                                                                                                         aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt                              >  1:5121193/1‑150 (MQ=255)
                                                                                                                          tCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGttt                             >  2:4988145/1‑150 (MQ=255)
                                                                                                                             gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                          >  2:979310/1‑150 (MQ=255)
                                                                                                                             gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                          >  1:1677932/1‑150 (MQ=255)
                                                                                                                             gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                          >  1:6678504/1‑150 (MQ=255)
                                                                                                                             gACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                          <  1:2994741/150‑1 (MQ=37)
                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcat                                                                                                               <  2:4220457/62‑1 (MQ=255)
                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                       >  2:511333/1‑150 (MQ=255)
                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                       >  2:7732837/1‑150 (MQ=255)
                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGAAAGTTTCAaaca                       >  2:7671977/1‑150 (MQ=255)
                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAAAGTAAGTTTCAaaca                       >  1:6986146/1‑150 (MQ=255)
                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGAGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                       >  1:5587022/1‑150 (MQ=255)
                                                                                                                                aaaattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                       >  1:1189921/8‑150 (MQ=37)
                                                                                                                                 acatTTTGGTGTGGTGGGGCCTCTGCGCGGGAACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa                      <  1:8139219/150‑1 (MQ=255)
                                                                                                                                 acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa                      >  1:808252/1‑150 (MQ=255)
                                                                                                                                 aaatttttGTAGGGTCGGGCCTCCCCGGGGCCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGATAGCTGAGCTAAACGTAAGTTTCAaacaa                      <  1:3649191/146‑1 (MQ=255)
                                                                                                                                  caaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac                     >  1:5812002/1‑150 (MQ=255)
                                                                                                                                  aattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac                     >  1:4694327/6‑150 (MQ=38)
                                                                                                                                   aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt                    >  1:7083304/1‑150 (MQ=255)
                                                                                                                                   aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCAACCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt                    >  2:6075372/1‑150 (MQ=255)
                                                                                                                                    aTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTg                   >  2:1555301/1‑150 (MQ=255)
                                                                                                                                    aTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTg                   >  2:9091468/1‑150 (MQ=255)
                                                                                                                                     ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACGGt                  >  2:7853607/1‑150 (MQ=255)
                                                                                                                                      tttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGtt                 >  1:1407600/1‑150 (MQ=255)
                                                                                                                                      ttggtcgggtcggGCCTCTGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGtt                 <  2:5074692/147‑1 (MQ=255)
                                                                                                                                          tCTGGTGGGGCCTCTGCGGGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGAt             <  2:8587860/150‑1 (MQ=255)
                                                                                                                                          tCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGAt             >  2:7918406/1‑150 (MQ=255)
                                                                                                                                           cTGGTCGGGCCTCTTCGGGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATc            <  2:2120061/150‑1 (MQ=255)
                                                                                                                                           cTGGTCGGGCCTCTCCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTATCAAACAGTTGTTTGaaa            <  2:6706519/150‑3 (MQ=255)
                                                                                                                                            tGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCg           >  1:4381200/1‑150 (MQ=255)
                                                                                                                                              gTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAAATGTTTGATCGGc         >  1:8459563/1‑150 (MQ=255)
                                                                                                                                               tCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGc         >  2:4840708/1‑149 (MQ=255)
                                                                                                                                               tCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCt        >  1:1291141/1‑150 (MQ=255)
                                                                                                                                                cGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAATCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCtt       >  1:2342324/1‑150 (MQ=255)
                                                                                                                                                 gggCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTc      >  1:5598410/1‑150 (MQ=255)
                                                                                                                                                 gggCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTc      >  2:2051739/1‑150 (MQ=255)
                                                                                                                                                  ggCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTcc     >  2:8276243/1‑150 (MQ=255)
                                                                                                                                                   gCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCg    >  2:3721748/1‑150 (MQ=255)
                                                                                                                                                   accTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCg    >  2:2281443/2‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCGCC  >  NC_002947/444977‑445274

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: