Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,108 A→T 22.3% intergenic (+78/‑210) 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,1080AT22.3% 365.9 / 90.2 172intergenic (+78/‑210)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (58/74);  new base T (7/32);  total (65/107)
Fisher's exact test for biased strand distribution p-value = 4.42e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 197 total aligned reads displayed.

GTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTG  >  NC_002947/444959‑445257
                                                                                                                                                      |                                                                                                                                                      
gTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGAcaacaa                                                                                                                                                        <  1:508795/150‑1 (MQ=255)
gTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGAcaacaa                                                                                                                                                        <  1:3919186/150‑1 (MQ=255)
  cGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtt                                                                                                                                                      <  2:8367447/150‑1 (MQ=255)
  cGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtt                                                                                                                                                      <  2:4803577/150‑1 (MQ=255)
  cGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtt                                                                                                                                                      <  1:9126814/150‑1 (MQ=255)
    aGCTCGGCGCACAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACTCCTACTCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGTCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                    <  2:6943357/150‑1 (MQ=255)
    aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                    <  1:5238578/150‑1 (MQ=255)
       tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTc                                                                                                                                                <  1:3579079/150‑1 (MQ=255)
       tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTc                                                                                                                                                <  1:556594/150‑1 (MQ=255)
         ggCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTg                                                                                                                                              <  2:3847882/150‑1 (MQ=255)
         ggCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTg                                                                                                                                              >  2:7537879/1‑150 (MQ=255)
          gcgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTga                                                                                                                                             <  2:2390159/150‑2 (MQ=255)
          gcgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTga                                                                                                                                             <  2:386023/150‑2 (MQ=25)
           cgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTGTGGt                                                                                                                                            >  2:5225636/1‑150 (MQ=255)
            gctaaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATc                                                                                                                                           <  2:6678371/147‑1 (MQ=25)
            gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTc                                                                                                                                           <  2:7992943/150‑1 (MQ=255)
            gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGACCAGACAAAATTGTTGTCTGATc                                                                                                                                           <  1:4890969/150‑1 (MQ=38)
               aaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCggg                                                                                                                                        <  2:1108682/150‑1 (MQ=255)
                aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                       <  1:8444266/150‑1 (MQ=255)
                aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                       <  1:7269408/150‑1 (MQ=255)
                aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                       <  1:7717578/150‑1 (MQ=255)
                aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                       <  1:5265480/150‑1 (MQ=255)
                aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGACCAGACAAAATTGTTGTCTGATCggac                                                                                                                                       <  1:9028645/150‑3 (MQ=38)
                aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGCCCGACCAGACAAAATTGTTGTCTGATCggac                                                                                                                                       <  1:784197/150‑3 (MQ=37)
                    cgtcAAGGCATCCGCCGCAATGCGGGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTAATTCAGCGCCGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCctct                                                                                                                                   <  1:6496812/150‑1 (MQ=255)
                      tcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTGCCCCCGAAAGCGTCCGGGCGGACAACAAATATTTGTGGTGGGGccgcgtc                                                                                                                                 >  1:6800816/1‑145 (MQ=255)
                        aaGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTgg                                                                                                                                             >  1:6700098/1‑136 (MQ=255)
                          tgCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                             <  2:8511761/149‑1 (MQ=255)
                          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCTCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgcgt                                                                                                                             <  2:2562934/150‑2 (MQ=255)
                          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                             <  2:3465918/150‑1 (MQ=255)
                          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                             <  1:7031647/150‑1 (MQ=18)
                          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                             <  1:4546444/150‑1 (MQ=18)
                          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                             <  2:7880979/150‑1 (MQ=18)
                          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                             <  1:9099582/150‑1 (MQ=18)
                          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                             <  1:2370291/150‑1 (MQ=18)
                          ggCATCCGACGAAATGCTTGGCCGGATCATCGACACATACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                             <  2:2055880/150‑1 (MQ=255)
                           gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCggg                                                                                                                            <  2:4886023/150‑1 (MQ=255)
                           gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCggg                                                                                                                            <  1:5705929/150‑1 (MQ=18)
                             aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCa                                                                                                                          <  2:7224138/150‑1 (MQ=255)
                                cgacgaAATGCTTGGCCGGATCTTCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                                       <  1:41504/150‑1 (MQ=255)
                                cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGc                                                                                                                       <  1:3619271/150‑1 (MQ=17)
                                      aaTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTTCCCCCGAATGCGGCCGATCAGAAAACAATTTTTTGTTGGCGGGGCGCGTCGCGGGCCAGcccccac                                                                                                                 >  2:5089305/1‑146 (MQ=255)
                                            tGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                    <  1:1511749/110‑1 (MQ=255)
                                            tGGCCGGACCATCGACACCTACGCCTGAATCCTGGCTTCATCCCCTATTTCCGCGCCGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGACCGATCAGACAACAATTTTGTGTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                           <  2:7632886/150‑1 (MQ=37)
                                                     catcGACACCTACGCCTGAATACTTGCGTCATCCCCTCTTTCAGCGCCTGCCTTGTAGTCGCGGGCTTGCCCGCGAATGCGTCGGACCACACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGATtgtg                                                                                                  <  2:8205472/150‑1 (MQ=38)
                                                       tcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCCCGCAAGCGGCCCCTCAGACAACAAAGTTTGCGGGGCGGGCCGCGTCGCGGGCCACCCCCCGCCAACAACGGGGTGtggc                                                                                                >  2:3748259/1‑150 (MQ=255)
                                                          acacCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                       >  1:286698/1‑108 (MQ=255)
                                                          acacCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                       <  2:286698/108‑1 (MQ=255)
                                                           cacCTCCGCGTTATTCCTGGCGTCAGCCCCTTTTCCCGCGCCTGCCTTGTCGGCGGGGGCTGGCCCGCTACTGCGCCCGAACACAAAAAAATTTGGTCTGGGCGGGCCACTGCGGGGGAACGCCCGCTCCTACAAGGGTTTGtggcatgg                                                                                            <  2:8073700/150‑1 (MQ=255)
                                                                      gAATCCTTGCGTCATCCCCTATTTCCGCGCCTGCCTGGTAGGCGCGGGCTTGCCCGCGAATGCGTCCGACCAGACAACAATTTTGTCTGGTGGGGCCTCTGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAtt                                                                                 <  2:6640805/150‑1 (MQ=38)
                                                                                   aTCCCCTATTCCCGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGACTGCGTCCGATCAGACACCAATTTTGTGTGGTCGGGCCTCTGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGt                                                                    <  2:8194288/150‑1 (MQ=38)
                                                                                    cccccAATTCCCGCGCGGGCCTTGTAGGCGGGGGCTTGCCCGCGATTGCTGCCGAACCCCAAACATTGTTGTGTGGTCGGGCCTCTCCGGGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTg                                                                   <  2:8689215/149‑1 (MQ=255)
                                                                                       ccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAATTTGTCTGGTCGGGCCTCGTAGCGGGCAATCCCGCTCCTACAAGGGTTTGTTGCATGGCATGACTGATTTGCCAATTTTCGTTGCg                                                                >  2:9155378/1‑150 (MQ=38)
                                                                                       ccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCCAAATCGTCCCATCAGTCAACAAATTTGTCTGGGAGGGCCGCGTGGCGGGCAAGCCCCCTCCGACAAGGGTGTTTGGCAAGGCATGAATGATTTGGCAATTTTAGTGGCg                                                                >  2:4554678/1‑150 (MQ=255)
                                                                                       ccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCCCGAATGCGGACGCTCCCACAACAAATTTTTGTGGGGGGGCCTCGTCGCGGGCCAGCCCCCTCCCACAAAGGGGTGTGGGCTGGCATGTCTGAATTGCCAATTTTCCTGGCg                                                                >  2:3649191/1‑150 (MQ=255)
                                                                                           ttcccgcgcCTGCCTTTTCGGGGCGGGCGGCCCGGCTACTGCGCCCCCCCCCACACCAAGTGTGTCGGGTGGGCCCTCGACGGGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                            <  2:1922313/145‑1 (MQ=255)
                                                                                            ttCCGCGCCGGCCTTGTCGGCGGGGGCTGCCCCGCCACTGCGCCCGAACACACAAAATTTTGTTCGGGTCGGGCCTCTCCGGGGCCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAAt                                                           <  2:3296616/150‑1 (MQ=255)
                                                                                                     gcctttgttgGCGCGGGCGGCCCCGCGACTGCGTCCCAACCAACACCATTTTGGGGTGGTGGGGCCTCCCCGCGGGCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATa                                                  <  1:7537879/146‑1 (MQ=255)
                                                                                                             tggAGCGGGGTTGCCCGCTAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCATGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATcc                                          <  2:6740798/149‑1 (MQ=255)
                                                                                                                gcgggggTTGCCCGGGTATTGCGTCCTATCACACAAAATTTTTGTGGGGCGGGGCCTCTGCGCGGCCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                        <  1:2986595/147‑1 (MQ=255)
                                                                                                                 gCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTTTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTTCCATTTTTCGTGGCGTAAATTGCACCATAGTCCATCATTGa                                      >  1:7998886/1‑150 (MQ=255)
                                                                                                                 gCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGGTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCCTAGTCCATCCTTgc                                      >  2:5801179/1‑149 (MQ=255)
                                                                                                                       ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                >  1:7281445/1‑150 (MQ=255)
                                                                                                                          cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                             >  2:8784295/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)
                                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                            >  1:444634/1‑150 (MQ=255)
                                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                            >  1:6504071/1‑150 (MQ=255)
                                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                            >  2:662023/1‑150 (MQ=255)
                                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                            >  1:1901012/1‑150 (MQ=255)
                                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                            >  2:3094776/1‑150 (MQ=255)
                                                                                                                            cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                           >  2:1368129/1‑150 (MQ=255)
                                                                                                                               gAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCt                        >  2:1437650/1‑150 (MQ=255)
                                                                                                                                      gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt                  >  2:400243/1‑150 (MQ=255)
                                                                                                                                       tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa                 >  2:9122252/1‑150 (MQ=255)
                                                                                                                                       tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa                 >  1:140473/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:7055818/86‑1 (MQ=18)
                                                                                                                                           aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt             >  1:6246004/1‑150 (MQ=255)
                                                                                                                                            tCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGttt            >  2:4988145/1‑150 (MQ=255)
                                                                                                                                               gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa         >  2:6478689/1‑150 (MQ=255)
                                                                                                                                               gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa         >  2:979310/1‑150 (MQ=255)
                                                                                                                                               gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa         >  1:1677932/1‑150 (MQ=255)
                                                                                                                                                  aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca      >  2:7732837/1‑150 (MQ=255)
                                                                                                                                                  aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca      >  1:5587028/1‑150 (MQ=255)
                                                                                                                                                  aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca      >  2:511333/1‑150 (MQ=255)
                                                                                                                                                  aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca      >  2:7286965/1‑150 (MQ=255)
                                                                                                                                                  aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAAAGTAAGTTTCAaaca      >  1:6986146/1‑150 (MQ=255)
                                                                                                                                                  aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTAGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca      >  2:7671976/1‑150 (MQ=255)
                                                                                                                                                  aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGAAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca      >  2:7672240/1‑150 (MQ=255)
                                                                                                                                                   acatTTTGGTGTGGTGGGGCCTCTGCGCGGGAACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa     <  1:8139219/150‑1 (MQ=255)
                                                                                                                                                   acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa     >  1:808252/1‑150 (MQ=255)
                                                                                                                                                     aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGATGAGCTAAACGTAAGTTTCAAACAACt   >  1:7083362/1‑150 (MQ=255)
                                                                                                                                                     aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCAACCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt   >  2:6075372/1‑150 (MQ=255)
                                                                                                                                                      aTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTg  >  2:1555301/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTG  >  NC_002947/444959‑445257

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: