Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,127 C→A 13.4% intergenic (+97/‑191) 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,1270CA13.4% 440.5 / 55.7 173intergenic (+97/‑191)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (93/55);  new base A (10/13);  total (103/70)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

Only 100 of 184 total aligned reads displayed.

CGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCGCCA  >  NC_002947/444979‑445275
                                                                                                                                                      |                                                                                                                                                    
 gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCtt                                                                                                                                                    <  2:7671043/150‑1 (MQ=255)
 gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCtt                                                                                                                                                    >  2:8354663/1‑150 (MQ=255)
    aaGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTTCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgc                                                                                                                                                 <  1:7924046/150‑1 (MQ=255)
      ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                               <  2:7180481/150‑1 (MQ=255)
      ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                               <  2:3904066/150‑1 (MQ=255)
      ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                               <  2:3488369/150‑1 (MQ=18)
      ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                               <  1:4546502/150‑1 (MQ=18)
      ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                               <  1:2370291/150‑1 (MQ=18)
      ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                               <  1:7031647/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)
                                                   atgcggggCGTCACCCCCTATGCCCGGCCCGGCGTTGTCGGCGGGGGGTGCCCCGCGACGGCGCCCTCTCACACAACACTTGTGTGTGGCGGGGCCGCGGGGGGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                                  <  2:7900112/143‑1 (MQ=255)
                                                             tCATCCCCTATTCCCGCGCCTGCCTTGTCGGAGCGGGCTTGCCCGCGAATGCGTCCGCTCAGAAAACAATTTTGTGTGGTCGGGCCTCTGCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTc                                                                                        <  2:7551977/150‑1 (MQ=255)
                                                               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)
                                                                   ccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCATTCCCGCGAAAGCCTACCATAAGACAACAATTATGTCTGGTCGGGCCGCTTCGCCGGCAAGACCCCACCGACAAGGGTTTTTGGCAAGGCATGACTGATTTTCCAATTTTCGTGGCg                                                                                  >  2:7524415/1‑150 (MQ=255)
                                                                        ttCCGCGCCGGCCTTGTCGGCGGGGGCTGCCCCGCCACTGCGCCCGAACACACAAAATTTTGTTCGGGTCGGGCCTCTCCGGGGCCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAAt                                                                             <  2:3296616/150‑1 (MQ=255)
                                                                                 tgcttTGTCGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGAAAACATTTTGGTCTGGTCGGGCCTCTTTGCGGGAAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATa                                                                    <  1:2474849/150‑1 (MQ=255)
                                                                                 gcctttgttgGCGCGGGCGGCCCCGCGACTGCGTCCCAACCAACACCATTTTGGGGTGGTGGGGCCTCCCCGCGGGCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATa                                                                    <  1:7537879/146‑1 (MQ=255)
                                                                                         tggAGCGGGGTTGCCCGCTAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCATGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATcc                                                            <  2:6740798/149‑1 (MQ=255)
                                                                                         tggAGCGGGCTTGCGCTCTAAGGCTTCCTCTCCTCAAACATTTTGGTCGGGGGGGGCCTCTTGGCGGGAATGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATcc                                                            <  2:1968839/149‑1 (MQ=255)
                                                                                           gcgggggCTGGCCCGCCACTGCGCCCCACCACACAACAATTGTGTGGGGGCGGCCCTCTCCGCGGGCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                                          <  1:6742231/146‑1 (MQ=255)
                                                                                           gAGCGGGCTTGCCCGCTAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTCCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                                          <  2:4564600/150‑1 (MQ=255)
                                                                                            gcgggggTTGCCCGGGTATTGCGTCCTATCACACAAAATTTTTGTGGGGCGGGGCCTCTGCGCGGCCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                                          <  1:2986595/147‑1 (MQ=255)
                                                                                             gCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGGTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCCTAGTCCATCCTTgc                                                        >  2:5801179/1‑149 (MQ=255)
                                                                                                 gCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                                                    >  2:694331/1‑150 (MQ=255)
                                                                                                   ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                                  >  1:7281445/1‑150 (MQ=255)
                                                                                                      cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                               >  2:7230839/1‑150 (MQ=255)
                                                                                                      cccGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                               >  1:5512363/1‑150 (MQ=18)
                                                                                                       ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                                             <  1:7185886/71‑1 (MQ=255)
                                                                                                       ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                              >  1:444634/1‑150 (MQ=255)
                                                                                                       ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                              >  1:6503236/1‑150 (MQ=255)
                                                                                                        cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                             >  2:1368129/1‑150 (MQ=255)
                                                                                                           gAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCt                                          >  2:1437650/1‑150 (MQ=255)
                                                                                                              agaggcccGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTaa                                        >  1:9042438/7‑150 (MQ=21)
                                                                                                                  gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTaa                                        >  2:1601660/1‑146 (MQ=255)
                                                                                                                  gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt                                    >  1:1923021/1‑150 (MQ=255)
                                                                                                                   tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa                                   >  1:7493975/1‑150 (MQ=255)
                                                                                                                    ccGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                             >  2:1263837/1‑139 (MQ=255)
                                                                                                                    ccGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                             <  1:1263837/139‑1 (MQ=255)
                                                                                                                    ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                                  <  1:7055818/86‑1 (MQ=18)
                                                                                                                    ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                                  >  2:7055673/1‑86 (MQ=18)
                                                                                                                     cGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAg                                 >  1:7225067/1‑150 (MQ=38)
                                                                                                                      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)
                                                                                                                         cAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc                             >  2:5219736/1‑150 (MQ=255)
                                                                                                                          aGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCa                            >  2:2677624/1‑150 (MQ=255)
                                                                                                                           gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                           >  2:979310/1‑150 (MQ=255)
                                                                                                                           gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                           >  2:6478689/1‑150 (MQ=255)
                                                                                                                           gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                           >  1:1677932/1‑150 (MQ=255)
                                                                                                                           gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                           >  1:6678504/1‑150 (MQ=255)
                                                                                                                              aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcat                                                                                                                >  1:4220457/1‑62 (MQ=255)
                                                                                                                              aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAtaca                        >  1:3635340/1‑150 (MQ=255)
                                                                                                                              aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                        >  2:7732837/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)
                                                                                                                              aaaattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                        >  1:1189921/8‑150 (MQ=37)
                                                                                                                              aaaattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                        >  1:4642264/8‑150 (MQ=37)
                                                                                                                               acatTTTGGTGTGGTGGGGCCTCTGCGCGGGAACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa                       <  1:8139219/150‑1 (MQ=255)
                                                                                                                               acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa                       >  1:808252/1‑150 (MQ=255)
                                                                                                                                aattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac                      >  2:4183467/6‑150 (MQ=38)
                                                                                                                                aattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac                      >  2:8556896/6‑150 (MQ=38)
                                                                                                                                 aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGATGAGCTAAACGTAAGTTTCAAACAACt                     >  1:7083362/1‑150 (MQ=255)
                                                                                                                                 aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCAACCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt                     >  2:6075372/1‑150 (MQ=255)
                                                                                                                                  aTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTg                    >  2:9091468/1‑150 (MQ=255)
                                                                                                                                   ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACGGt                   >  2:7853607/1‑150 (MQ=255)
                                                                                                                                    tttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGtt                  >  1:1407600/1‑150 (MQ=255)
                                                                                                                                    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             >  2:345749/1‑150 (MQ=255)
                                                                                                                                          cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATc             >  1:8031584/1‑150 (MQ=255)
                                                                                                                                          cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTATCAAACAGTTGTTTGaaa             <  2:6706467/150‑3 (MQ=255)
                                                                                                                                             gTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGc          >  1:3131552/1‑150 (MQ=255)
                                                                                                                                             gTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAAATGTTTGATCGGc          >  1:8459563/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)
                                                                                                                                                gggCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTc       >  2:2561397/1‑150 (MQ=255)
                                                                                                                                                  accTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCg     >  2:2281443/2‑150 (MQ=255)
                                                                                                                                                     tctTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCGcgc  >  2:5273713/1‑148 (MQ=255)
                                                                                                                                                     tctTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCGCCa  >  1:4009359/1‑150 (MQ=255)
                                                                                                                                                     tctTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCGCCa  >  2:378134/1‑150 (MQ=255)
                                                                                                                                                     tctTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCGCCa  >  1:5167985/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTCCGCCA  >  NC_002947/444979‑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: