Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,123 G→A 19.9% intergenic (+93/‑195) 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,1230GA19.9% 340.0 / 63.5 150intergenic (+93/‑195)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (79/39);  new base A (10/17);  total (89/61)
Fisher's exact test for biased strand distribution p-value = 7.63e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 154 total aligned reads displayed.

AGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTC  >  NC_002947/444975‑445270
                                                                                                                                                    |                                                                                                                                                   
aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                                    <  1:1621827/150‑1 (MQ=255)
   tcgtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCctc                                                                                                                                                 <  1:3335383/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)
          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGGGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                          <  1:3226214/150‑1 (MQ=18)
          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                          <  2:3344718/150‑1 (MQ=255)
          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                          <  2:4642513/150‑1 (MQ=255)
          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                          <  1:2231502/150‑1 (MQ=255)
          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                          <  1:6674281/150‑1 (MQ=18)
          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                          <  2:5280060/150‑1 (MQ=18)
          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                          <  1:2265969/150‑1 (MQ=18)
          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                          <  2:6789821/150‑1 (MQ=18)
          ggCATCCGACGAAATACTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                          <  2:4665078/150‑1 (MQ=18)
                cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGATCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                                                    >  2:4314774/1‑150 (MQ=255)
                 gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGTTTCCCGCGAATTCGTGCGATAAGACAAAAAATTTTTGTTGTGGGGGCTCTTAGCGGGGAAGcc                                                                                                                                   >  2:4646372/1‑150 (MQ=255)
                 gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGcc                                                                                                                                   <  1:1272922/150‑1 (MQ=17)
                 gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGCAGGAGCGGGCGTGCCCGCGGATGCGTCCGATGAGACAACAATTTTGTATGGTGGGGGCGCGTGGCGGGCAAGcc                                                                                                                                   >  2:3342643/1‑150 (MQ=255)
                  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)
                           ttGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCCCGAAAGCGTCCGAACACACAACAATTTTGTATGGTAGGGCCTCGACGCGGGCAAGCCCCCTCCAACa                                                                                                                         >  2:6200567/1‑150 (MQ=255)
                                  gATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCCGCGCGTGCCTTGTCGGTGGGGGCTTGCCCGCGAATGCGTCCGATCCGACAACAATTTTGTCTGGTGGGGCCTCTTCGGGGGCAAGCCCGCTCCTACAAGGGttt                                                                                                                  <  2:1639732/150‑1 (MQ=38)
                                    tcatcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCCGCGCCGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGAGGGCAAGCCCGCTCCTACAAGGGTTtgt                                                                                                                <  1:6475090/150‑1 (MQ=255)
                                     aatcGACACCCACGCCTCAATACTTGCGACACCCCCTATTCCAGCGCCTGCCTTGTTGGCGCGGGCTTGCCCGCTAATGCGTCCGATCACAAAAAATTTTTGTCTGGTCGGGCCTCTGCGCGGGCACGCCCGCTCCTACAAGGGTTtgtg                                                                                                               <  2:1007000/149‑1 (MQ=255)
                                          ccacCTCCGCCTTAATCCTTGCGTCACCCCCTATTCCCGCGCCGGCCTTGTCGGAGCGGGCTGGCCCGCGAATGCGCCCTACCAGAAAACAATTTTGCCGGGTCGGCCCACTTCGCGGGAACCCCCGCTCCTACAAGGGTTTGtggcatg                                                                                                          <  2:1043694/149‑1 (MQ=255)
                                              cTACGCCTTATTACTTGCTTCATCCCCTATTTCTGCGCCGGCTTTGTTGGTGGGGGATTGCCCGCGATTGCGTCCTATCATAAAACATTTTTGTCTGGGGGGGCTTCTTCGGGGGCATGCCCGCTCCTACAAGGGTTTGTggcatggcat                                                                                                      <  2:1095043/150‑1 (MQ=255)
                                               tACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTTCCCGCGAAATCCTCCCCACAGACAACAAATTTTTGTGGGGGGGGCGCGTGGGGGGGAAGCCCCCGCCAAAAACGGGTTGTGggaatgcatg                                                                                                     >  2:1551938/1‑150 (MQ=255)
                                               cccGCCTTATTCCTGGCGTCAGCCCCTATTCCCGCGCCGGCCTTGTTGGTGGGGGGGTGCCCGCTAATGCGTCCCATCACACAACAATTTTGTGTGGTGGGGCCTCTTGGCGGGAACGCCCGCTCCTACAAGGGTTTGTGgcatggcatg                                                                                                     <  2:915010/148‑1 (MQ=255)
                                                       aaTACTTGCTTCTTCCCCTTTTTCCGCGCCTGCCTTGTAGGTGCGGGCTGGCCCGCGATTGCGTCCGATCAGACAACAATTTTGTCTGGTGGGGCGGCTTCGCGGGCAAGCCAGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                             <  2:6548226/150‑1 (MQ=255)
                                                                      ccttttttCCGCGCCTGCCTGTTTGGAGGGGGTTGGCCGGCTACGGCTCCCGATAACCAAACATTTTTGTGGGGTGGGGCATCTTCGGGGGAAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                                              <  2:2980893/145‑1 (MQ=255)
                                                                      cccTCTTTCCGCGCCTGCCTTGTCGGCGCGGGCTTGCCCGCGAATGCGGCCGATCTGACACCAATTTTGTCTGGGCGGGCCTGTCCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                                              <  2:4450460/150‑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)
                                                                        cTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCCCGGATGCGTGCCGTAACACAACAAGTTTGTCTGGGGCGGCCTGATCGCGGGCAAGCCCGCGCCAACAAGGGTTTGTGGCATGGCAAGACTGATTTTCCAATTTTAGGGGCGt                                                                            >  2:2663887/1‑150 (MQ=255)
                                                                         tATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTTCCCCCCAATTCGTCCGATAACACAACAATGTTTTCTGGGCGGGCCCCATCGCCGGCAAGCCCCCGCCGACAAGGGGTTGTGGCATGGCCAGACTGATTTGCCAAATTTAGTGGCGTa                                                                           >  2:5125112/1‑150 (MQ=255)
                                                                                gcgcCGGCCGTGTTGGCGGGGGCTTGCCCGCGAATGCGTCCGACCAGACACAAATTTTGTCGGGTGGGGCCTCTGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCa                                                                    <  2:3782753/150‑1 (MQ=38)
                                                                                            ttgttgggggTTTGCCCGCTATTTCGTCCTTTTTTACAAAATTTTTGTGTGGTGGGGCTTCTGGGGGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATc                                                        <  2:3039172/143‑1 (MQ=255)
                                                                                                aGCGGGCTTGCCCGCGACTGGGTCCGATCAGACAACAATTTGGTGTGGTCGGGCCTCTTCGGGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTg                                                    <  2:4414413/150‑1 (MQ=255)
                                                                                                       ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                             >  1:2485132/1‑150 (MQ=255)
                                                                                                       ttGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                             >  1:2404310/1‑150 (MQ=17)
                                                                                                       ttGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                             >  1:5089658/1‑150 (MQ=17)
                                                                                                         gCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc                                           >  2:4254679/1‑150 (MQ=255)
                                                                                                          cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGGCGGGCCTCTTCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                          <  2:103982/150‑1 (MQ=255)
                                                                                                          cccGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                          >  2:3285808/1‑150 (MQ=18)
                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                         >  1:6219391/1‑150 (MQ=255)
                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                         >  1:4788429/1‑150 (MQ=255)
                                                                                                            cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                        >  1:914217/1‑150 (MQ=255)
                                                                                                                  gaggcccGACCAGACAAAATTGTTGTCTGATCGGACCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTaaa                                  >  2:6301695/6‑150 (MQ=37)
                                                                                                                      tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa                              >  1:1822345/1‑150 (MQ=255)
                                                                                                                       ccGATCAGACAACAATTTTGTCTGGTCGGGCCTCTGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa                             <  2:6799500/150‑1 (MQ=255)
                                                                                                                       ccGACCAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa                             >  1:4927417/1‑150 (MQ=255)
                                                                                                                       ccGACCAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa                             >  1:4927402/1‑150 (MQ=255)
                                                                                                                       ccGACCAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa                             >  2:6045023/1‑150 (MQ=255)
                                                                                                                       ccGACCAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCATACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTATTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa                             >  2:6045062/1‑150 (MQ=38)
                                                                                                                         gaccAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGt                           >  1:3714598/4‑150 (MQ=255)
                                                                                                                         gATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGt                           >  2:2911738/1‑150 (MQ=255)
                                                                                                                          accAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt                          >  1:3482458/3‑150 (MQ=255)
                                                                                                                          aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt                          >  2:887714/1‑150 (MQ=255)
                                                                                                                            cAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc                        >  1:186560/1‑150 (MQ=255)
                                                                                                                              gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                      <  2:4901911/150‑1 (MQ=255)
                                                                                                                              gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                      >  2:6278295/1‑150 (MQ=255)
                                                                                                                              gACAACAATATAGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAAAGTAAGTTTCaa                      >  2:220326/1‑150 (MQ=255)
                                                                                                                               acaacaATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa                     >  1:5512339/1‑150 (MQ=255)
                                                                                                                               acaacaATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa                     >  1:5512340/1‑150 (MQ=255)
                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGTCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                   >  2:2762821/1‑150 (MQ=255)
                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                   >  2:4886736/1‑150 (MQ=255)
                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                   >  1:6386098/1‑150 (MQ=255)
                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                   >  1:3850985/1‑150 (MQ=255)
                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                   >  1:6522327/1‑150 (MQ=255)
                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                   >  2:3454717/1‑150 (MQ=255)
                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                   >  1:681356/1‑150 (MQ=255)
                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTGGTGGCATGGCATGACGGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCGTGTCAGAAAGCTGAGCTAGTCGTCAGTTTCattca                   >  1:2824678/1‑146 (MQ=255)
                                                                                                                                  acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa                  >  1:2155743/1‑150 (MQ=255)
                                                                                                                                  acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa                  >  1:5927914/1‑150 (MQ=255)
                                                                                                                                   caaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac                 >  2:855248/1‑150 (MQ=255)
                                                                                                                                    attgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt                >  2:191524/5‑150 (MQ=38)
                                                                                                                                    aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt                >  1:3859012/1‑150 (MQ=255)
                                                                                                                                    aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt                >  1:1013076/1‑150 (MQ=255)
                                                                                                                                    aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt                >  1:6836307/1‑150 (MQ=255)
                                                                                                                                    aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt                >  1:1011017/1‑150 (MQ=255)
                                                                                                                                      ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGt              >  1:5652764/1‑150 (MQ=255)
                                                                                                                                      ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGt              >  2:3224247/1‑150 (MQ=255)
                                                                                                                                        ttGTCTGGTCGGCCCTCTCCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGttt            <  2:2383739/150‑1 (MQ=255)
                                                                                                                                           tCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGAt         >  2:217932/1‑150 (MQ=255)
                                                                                                                                           tCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGAt         >  2:653334/1‑150 (MQ=255)
                                                                                                                                           tCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGAt         >  2:5272247/1‑150 (MQ=255)
                                                                                                                                             tGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCg       >  2:5752405/1‑150 (MQ=255)
                                                                                                                                             tGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCg       >  1:5963018/1‑150 (MQ=255)
                                                                                                                                              ggTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCgg      >  2:3475224/1‑150 (MQ=255)
                                                                                                                                               gTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGc     >  1:941112/1‑150 (MQ=255)
                                                                                                                                               gTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGc     >  1:1074141/1‑150 (MQ=255)
                                                                                                                                                tggggCCTCTCCGCGGGCAACCCCGCTCCTACTAGGGTTTGTGGCATGGCATGACTGTTTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCATTGTCAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCt    <  1:6200567/148‑1 (MQ=255)
                                                                                                                                                  gggCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTc  >  1:4567522/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTC  >  NC_002947/444975‑445270

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: