Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,126 T→C 13.9% intergenic (+96/‑192) 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,1260TC13.9% 417.7 / 49.3 173intergenic (+96/‑192)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (87/54);  new base C (10/13);  total (103/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 = 5.66e-02

Only 100 of 185 total aligned reads displayed.

GTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCGCCA  >  NC_002947/444977‑445275
                                                                                                                                                      |                                                                                                                                                     
gtcgtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCct                                                                                                                                                       <  2:2989247/150‑1 (MQ=255)
   gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCtt                                                                                                                                                    >  2:8354663/1‑150 (MQ=255)
   gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCtt                                                                                                                                                    <  2:7671043/150‑1 (MQ=255)
      aaGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTTCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgc                                                                                                                                                 <  1:7924046/150‑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)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                               <  2:7180481/150‑1 (MQ=255)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                               <  1:2370291/150‑1 (MQ=18)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                               <  1:7031647/150‑1 (MQ=18)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                               <  2:7880979/150‑1 (MQ=18)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                               <  1:4546502/150‑1 (MQ=18)
        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACATACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                               <  2:2055880/150‑1 (MQ=255)
         gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCggg                                                                                                                                              <  1:5705929/150‑1 (MQ=18)
              cgacgaAATGCTTGGCCGGATCTTCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                                                         <  1:41504/150‑1 (MQ=255)
               gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGcc                                                                                                                                        <  2:7059876/150‑1 (MQ=255)
                    aaTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTTCCCCCGAATGCGGCCGATCAGAAAACAATTTTTTGTTGGCGGGGCGCGTCGCGGGCCAGcccccac                                                                                                                                   >  2:5089305/1‑146 (MQ=255)
                                   catcGACACCTACGCCTGAATACTTGCGTCATCCCCTCTTTCAGCGCCTGCCTTGTAGTCGCGGGCTTGCCCGCGAATGCGTCGGACCACACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGATtgtg                                                                                                                    <  2:8205472/150‑1 (MQ=38)
                                     tcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCCCGCAAGCGGCCCCTCAGACAACAAAGTTTGCGGGGCGGGCCGCGTCGCGGGCCACCCCCCGCCAACAACGGGGTGtggc                                                                                                                  >  2:3748259/1‑150 (MQ=255)
                                       gACACCACCGCGTTAATAGTTGGGACAGCCCCTATTCCCGCGCGGGCTTTGTAGGCGGGGGGTGGCCGGCTATGGGGCCCGATCACACAAAAATTTGGTGGGGTGGGGCCTCTGCGGGGCCACCCCCGCTCCTACAAGGGTTTGtggcat                                                                                                                <  2:8730419/150‑1 (MQ=255)
                                         cacCTCCGCGTTATTCCTGGCGTCAGCCCCTTTTCCCGCGCCTGCCTTGTCGGCGGGGGCTGGCCCGCTACTGCGCCCGAACACAAAAAAATTTGGTCTGGGCGGGCCACTGCGGGGGAACGCCCGCTCCTACAAGGGTTTGtggcatgg                                                                                                              <  2:8073700/150‑1 (MQ=255)
                                                 ccTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCCCCAATGCGTCCCCTCAGACAAAAAATTTTTGTTGTGGGGCCGCGTCGCGGGCAACCCCCCGCCAAAAAGGGGTTTTTGGAAGGCATTAcag                                                                                                      >  2:7209568/1‑148 (MQ=255)
                                                    gAATACTTGCGACAGCCCCTTTTTCCGCGCCGGCCTTGTAGGTGGGGGCTTGCCCGCGACTGCGTCCTAACAGACAACAATTTGGTCTGGTCGGGCCTGTTCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAtt                                                                                                   <  2:6640801/150‑1 (MQ=255)
                                                               tCATCCCCTATTCCCGCGCCTGCCTTGTCGGAGCGGGCTTGCCCGCGAATGCGTCCGCTCAGAAAACAATTTTGTGTGGTCGGGCCTCTGCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTc                                                                                        <  2:7551977/150‑1 (MQ=255)
                                                                 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)
                                                                     ccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCCCGAATGCGGACGCTCCCACAACAAATTTTTGTGGGGGGGCCTCGTCGCGGGCCAGCCCCCTCCCACAAAGGGGTGTGGGCTGGCATGTCTGAATTGCCAATTTTCCTGGCg                                                                                  >  2:3649191/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)
                                                                                           tggAGCGGGCTTGCGCTCTAAGGCTTCCTCTCCTCAAACATTTTGGTCGGGGGGGGCCTCTTGGCGGGAATGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATcc                                                            <  2:1968839/149‑1 (MQ=255)
                                                                                             gAGCGGGCTTGCCCGCTAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTCCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                                          <  2:4564600/150‑1 (MQ=255)
                                                                                              gcgggggTTGCCCGGGTATTGCGTCCTATCACACAAAATTTTTGTGGGGCGGGGCCTCTGCGCGGCCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                                          <  1:2986595/147‑1 (MQ=255)
                                                                                               gCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGGTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCCTAGTCCATCCTTgc                                                        >  2:5801179/1‑149 (MQ=255)
                                                                                               gCGGGCTTGCCCGCGAAAGCGTCCGAACAGACAACAATTTTGTCTGGGCGGGCCTCTACGCGGGCAAGCCCCCTCCTACAAGGGGTTGTGGCATGGCATGACTGATTTGCCAATTTTCGTGGCGTAAATGGCACCAAAGTCCATACTTGa                                                        >  1:8819618/1‑150 (MQ=255)
                                                                                                     ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                                  >  1:7281357/1‑150 (MQ=255)
                                                                                                        cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                               >  2:7230839/1‑150 (MQ=255)
                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                                             <  1:7185886/71‑1 (MQ=255)
                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                              >  2:3094776/1‑150 (MQ=255)
                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                              >  1:444634/1‑150 (MQ=255)
                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                              >  1:6503236/1‑150 (MQ=255)
                                                                                                                agaggcccGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTaa                                        >  1:9042438/7‑150 (MQ=21)
                                                                                                                    gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt                                    >  1:1923021/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)
                                                                                                                      ccGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                             >  2:1263837/1‑139 (MQ=255)
                                                                                                                      ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                                  >  2:7055673/1‑86 (MQ=18)
                                                                                                                      ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                                  <  1:7055818/86‑1 (MQ=18)
                                                                                                                        gATCAGACAACAATTTTGTCTGGTCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGt                                >  2:300651/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)
                                                                                                                           accacaacaTTTTGGTCTGGGCGCGCCACTCCGGGGCCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc                             <  1:6807959/147‑1 (MQ=255)
                                                                                                                             gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                           >  2:6478689/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)
                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcat                                                                                                                <  2:4220457/62‑1 (MQ=255)
                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAtaca                        >  1:3635340/1‑150 (MQ=255)
                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                        >  2:7732837/1‑150 (MQ=255)
                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                        >  2:7286965/1‑150 (MQ=255)
                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                        >  2:511333/1‑150 (MQ=255)
                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                        >  1:5587028/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)
                                                                                                                                aaaattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                        >  1:4642264/8‑150 (MQ=37)
                                                                                                                                 acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGTCGGAATGCTGAGCTAAAGGTAAGTTTCaatctc                       >  1:8122393/1‑146 (MQ=255)
                                                                                                                                  aattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac                      >  2:8556896/6‑150 (MQ=38)
                                                                                                                                   aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt                     >  1:8377065/1‑150 (MQ=255)
                                                                                                                                   aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGATGAGCTAAACGTAAGTTTCAAACAACt                     >  1:7083362/1‑150 (MQ=255)
                                                                                                                                   aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCAACCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt                     >  2:6075372/1‑150 (MQ=255)
                                                                                                                                    aTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTg                    >  2:9091468/1‑150 (MQ=255)
                                                                                                                                    aTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTg                    >  2:1555301/1‑150 (MQ=255)
                                                                                                                                     ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACGGt                   >  2:7853607/1‑150 (MQ=255)
                                                                                                                                      tttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGtt                  >  1:1407600/1‑150 (MQ=255)
                                                                                                                                      tttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGtt                  <  2:6768158/150‑1 (MQ=255)
                                                                                                                                      tttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGtt                  <  1:979310/150‑1 (MQ=255)
                                                                                                                                          tCTGGTGGGGCCTCTGCGGGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGAt              <  2:8587860/150‑1 (MQ=255)
                                                                                                                                          tCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGAt              >  2:7918406/1‑150 (MQ=255)
                                                                                                                                           cTGGTCGGGCCTCTTCGGGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATc             <  2:2120061/150‑1 (MQ=255)
                                                                                                                                           cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATc             >  1:8031584/1‑150 (MQ=255)
                                                                                                                                           cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTATCAAACAGTTGTTTGaaa             <  2:6706467/150‑3 (MQ=255)
                                                                                                                                           cTGGTCGGGCCTCTCCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTATCAAACAGTTGTTTGaaa             <  2:6706519/150‑3 (MQ=255)
                                                                                                                                            tGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCg            >  1:4381200/1‑150 (MQ=255)
                                                                                                                                              gTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGc          >  1:3131552/1‑150 (MQ=255)
                                                                                                                                               tCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGc          >  2:4840708/1‑149 (MQ=255)
                                                                                                                                               tCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGc          <  1:4840708/149‑1 (MQ=255)
                                                                                                                                               tCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCt         >  1:1291141/1‑150 (MQ=255)
                                                                                                                                                   accTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCg     >  2:2281443/2‑150 (MQ=255)
                                                                                                                                                      tctTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCGCCa  >  1:4009359/1‑150 (MQ=255)
                                                                                                                                                      tctTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCGCCa  >  1:5257965/1‑150 (MQ=255)
                                                                                                                                                      tctTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCGCCa  >  2:5273651/1‑150 (MQ=255)
                                                                                                                                                      tctTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCGCCa  >  1:7455091/1‑150 (MQ=255)
                                                                                                                                                      tctTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCGCCa  >  2:378134/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCGCCA  >  NC_002947/444977‑445275

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: