Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,764,524 A→C 8.2% intergenic (+128/+28) trmL → / ← secB tRNA (cytosine(34)‑2'‑O)‑methyltransferase TrmL/protein‑export protein SecB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,764,5240AC 8.2% 702.1 / 40.2 215intergenic (+128/+28)trmL/secBtRNA (cytosine(34)‑2'‑O)‑methyltransferase TrmL/protein‑export protein SecB
Reads supporting (aligned to +/- strand):  ref base A (119/78);  new base C (15/3);  total (134/81)
Fisher's exact test for biased strand distribution p-value = 7.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 215 total aligned reads displayed.

TGATGTGTTGCCTGTGCGGGCCTCTTCGCGGGTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCGGGTCACCAGGC  >  NC_002947/5764394‑5764670
                                                                                                                                  |                                                                                                                                                  
tGATGTGTGGCGTGGGCGGGCCTCTCCGCGGGTAACCCCGCTCCTCCGCTGGCCGCGGGACCGGCGCAATCCCCGGACGGCGCGGCTTTACCCGCGAAGAGGCCCGCGCCGGTTTCAAGAAAGCGCCTCTCGGGGCGCTTTTTTgttgcc                                                                                                                                 <  2:4336179/150‑6 (MQ=12)
  aTGTGTTGCCTGTGCGGGCCTCTTCGCGGGTAAAACCCCTCCAACGCACGCCGAGGGAACGGCGCCAAACCCGGTAGGAGCGGGGTTAACCCCCAAGCGGCCCGCGCCCGCAACAAAAAAGCCCCCCGCGAGGCGCGGTCGTGAAAcccg                                                                                                                               >  2:6903938/1‑148 (MQ=255)
     tgtgGCGTGGGCGGGCCTCTTCGGGGGTAACCCCGCGCCTACCCCGGCCGAGGGAAGGGCGCACCCCCCTGTCGGAGGGGTTTTACCCGCTAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAg                                                                                                                            <  2:1114603/150‑1 (MQ=255)
         gCCTGTGCGGGCCTCTTCGCGGGTAAACCCGCTCCGAAGCAGGCCGCGGGAACGGCGCAATCCCCGGTAGGGGCGGGTTTAACCCCGAAGAGGCCCCCGCCGGCCACAAAAAAGCGCCCCCAGAGGCGCGTTCTTGGAACCGGCAGTAAt                                                                                                                        >  2:1652518/1‑150 (MQ=255)
          ccTGTGCGGGCCTCTTCGCGGTTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAAtt                                                                                                                       <  2:6714306/150‑1 (MQ=255)
                     CTCTTCGCGGGTAAAACCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCCGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAAACCGCCGTAATTACAGCACAGTa                                                                                                            <  2:2050612‑M1/150‑1 (MQ=255)
                      TCTTCGCGGGTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAg                                                                                                           <  2:5565779‑M1/150‑1 (MQ=255)
                      TCTTCGCGGGTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCGTTCTTGAAAACGGCAGTAATTACTGCACAGTAg                                                                                                           <  2:1995003‑M1/150‑1 (MQ=255)
                        ttcgcgggtaaacccgctcctacaggggattgcgccgttccctcggcctgcgtaggagcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGgc                                                                                                         >  2:6099526‑M2/84‑150 (MQ=255)
                         TCGCGGGTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACCGTAGgcg                                                                                                        <  2:1315631‑M1/150‑1 (MQ=255)
                           GCGGGTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCt                                                                                                      <  2:2487184‑M1/150‑1 (MQ=255)
                           GCGGGTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGGAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAATGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACCGGAGGCGCt                                                                                                      <  1:2933651‑M1/150‑1 (MQ=255)
                              GGTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCa                                                                                                   <  2:6000807‑M1/150‑1 (MQ=255)
                              GGTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCa                                                                                                   <  1:7358320‑M1/150‑1 (MQ=255)
                                 aatcccgctcctacaggggattgcgccgttccctcggcctgcgtaggagcgggtttacccgcgaagaggcccgcGCCGGCAACAAAATTGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCAccc                                                                                                >  2:4213892‑M2/75‑150 (MQ=255)
                                  AACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCg                                                                                               <  1:1065641‑M1/150‑1 (MQ=255)
                                    cccgctcctacaggggattgcgccgttccctcggcctgcgtaggagcgggtttacccgcgaagaggcccgcACAGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCt                                                                                             >  2:82433‑M2/72‑150 (MQ=255)
                                    CCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCt                                                                                             <  1:1186325‑M1/150‑1 (MQ=255)
                                     ccgctcctacaggggattgcgccgttccctcggcctgcgtaggagcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCtt                                                                                            >  1:625899‑M2/71‑150 (MQ=255)
                                     ccgctcctacaggggattgcgccgttccctcggcctgcgtaggagcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCtt                                                                                            >  1:3483859‑M2/71‑150 (MQ=255)
                                     ccgctcctacaggggattgcgccgttccctcggcctgcgtaggagcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCtt                                                                                            >  2:3399196‑M2/71‑150 (MQ=255)
                                     ccgctcctacaggggattgcgccgttccctcggcctgcgtaggagcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCtt                                                                                            >  1:4162707‑M2/71‑150 (MQ=255)
                                     ccgctcctacaggggattgcgccgttccctcggcctgcgtaggagcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCtt                                                                                            >  2:2862840‑M2/71‑150 (MQ=255)
                                     ccgctcctacaggggattgcgccgttccctcggcctgcgtaggagcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCtt                                                                                            >  1:2118829‑M2/71‑150 (MQ=255)
                                     ccgctcctacaggggattgcgccgttccctcggcctgcgtaggagcgggtttacccgcgaagaggcccgcACAGGCAACACAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCtt                                                                                            >  2:4424435‑M2/71‑150 (MQ=255)
                                     CCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCtt                                                                                            <  1:2325601‑M1/150‑1 (MQ=255)
                                      CGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTc                                                                                           <  1:1603772‑M1/150‑1 (MQ=255)
                                           ctccGCTGGCCGTGGGACGGGCGCAACCCCCTTTTGTTGGGGTTTAACCGGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGATACAGGCAGTAATTACTGCACAGTAGGCGCTTCCCCCGCTTCCTGCa                                                                                      <  2:3084375/147‑1 (MQ=255)
                                           CTACGCAGGCCGAGGGAACGGCGCAATCCCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTGTCTTGAAACCGGCAGTAATTACAGCACAGTAGGCGCTTCACCCGCTTCCAGCa                                                                                      <  1:2832372‑M1/150‑1 (MQ=255)
                                            TACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCAt                                                                                     <  2:5394669‑M1/150‑1 (MQ=255)
                                            TACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCAt                                                                                     <  2:5600134‑M1/150‑1 (MQ=255)
                                            TACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCAt                                                                                     <  1:4554879‑M1/150‑1 (MQ=255)
                                            TACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCAt                                                                                     <  2:425797‑M1/150‑1 (MQ=255)
                                                 gggattgcgccgttccctcggcctgcgtaggagcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTTATTACTTCACAGTAGGCGCTTCACCCGCTTCCTGCATCCGCt                                                                                >  2:1568935‑M2/59‑150 (MQ=255)
                                                       aGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCAttt                                                                          >  1:3595087/1‑150 (MQ=255)
                                                          gAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCtt                                                                       >  2:1604442/1‑150 (MQ=255)
                                                             cGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTgcg                                                                    >  1:2514370/1‑150 (MQ=255)
                                                              ggCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTgcgc                                                                   >  2:2533933/1‑150 (MQ=255)
                                                                 gcAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTa                                                                >  1:5566496/1‑150 (MQ=255)
                                                                  cAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTAc                                                               >  1:3331481/1‑150 (MQ=255)
                                                                   aaTCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACa                                                              >  2:3616855/1‑150 (MQ=255)
                                                                       cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGc                                                          >  2:3507705/1‑150 (MQ=255)
                                                                           ctcttcgcgggtaaacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCg                                                      >  1:1542333‑M2/33‑150 (MQ=255)
                                                                            taggagcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGa                                                     <  2:4555212‑M2/119‑1 (MQ=255)
                                                                               gagcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTCGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGt                                                  >  2:3542885‑M2/29‑150 (MQ=255)
                                                                               gagcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTc                                                       <  2:2984519‑M2/117‑1 (MQ=255)
                                                                               gagcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTc                                                       >  1:2984422‑M2/29‑145 (MQ=255)
                                                                               gagcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGt                                                  <  1:6099526‑M2/122‑1 (MQ=255)
                                                                                 gcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTgc                                                                     <  2:2916929‑M2/103‑1 (MQ=255)
                                                                                 gcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTg                                                >  1:3659706‑M2/27‑150 (MQ=255)
                                                                                 gcgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTg                                                >  1:2773775‑M2/27‑150 (MQ=255)
                                                                                  cgggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGa                                               >  1:1685477‑M2/26‑150 (MQ=255)
                                                                                   gggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGAc                                              <  2:2514370‑M2/126‑1 (MQ=255)
                                                                                   gggtttacccgcgaagaggcccgcACAGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGAc                                              <  2:2764997‑M2/126‑1 (MQ=255)
                                                                                    ggtttacccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGAcc                                             >  1:1024381‑M2/24‑150 (MQ=255)
                                                                                      tttacccgcgaataggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCgg                                           <  2:5700883‑M2/129‑1 (MQ=255)
                                                                                           ccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACa                                      >  2:4106477‑M2/17‑150 (MQ=255)
                                                                                           ccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACa                                      >  1:2398648‑M2/17‑150 (MQ=255)
                                                                                           ccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACa                                      >  1:7128951‑M2/17‑150 (MQ=255)
                                                                                           ccgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACa                                      >  1:3119937‑M2/17‑150 (MQ=255)
                                                                                            cgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAg                                     >  1:5375829‑M2/16‑150 (MQ=255)
                                                                                            cgcgaagaggcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAg                                     >  1:5319989‑M2/16‑150 (MQ=255)
                                                                                                     gcccgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAgcgc                            >  2:5843044‑M2/7‑150 (MQ=255)
                                                                                                       ccgcGCCGGCAACTAAAAAGCGCCCCGAGAGGCGCTTTCTTGAATCCGGCAGTAATTACTGCACAGTACGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCg                          <  1:6581222‑M2/146‑1 (MQ=255)
                                                                                                        cgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCgg                         <  1:6684621‑M2/147‑1 (MQ=255)
                                                                                                        cgcGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCgg                         >  2:1294446‑M2/4‑150 (MQ=255)
                                                                                                          cGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGa                       >  1:6389050‑M2/2‑150 (MQ=255)
                                                                                                            ccGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaac                     >  2:5887060/1‑150 (MQ=255)
                                                                                                             cGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaacg                    >  2:5423540/1‑150 (MQ=255)
                                                                                                             cGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaacg                    >  2:7171629/1‑150 (MQ=255)
                                                                                                               gtttcaaGAAAGCGCCTCTCGGGGCGCTTTTTTGTTGCCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaacgaa                  >  1:1044649/5‑150 (MQ=37)
                                                                                                               gtttcaaGAAAGCGCCTCTCGGGGCGCTTTTTTGTTGCCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaacgaa                  >  1:4772576/5‑150 (MQ=37)
                                                                                                               gCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAg                                     <  1:5644841/131‑1 (MQ=255)
                                                                                                                caacaaAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaacgaac                 >  2:5042044/1‑150 (MQ=255)
                                                                                                                 aacaaAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTATGCGCTACACACGCTTCCTGCAAAAGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaacgaaca                >  2:2323542/1‑149 (MQ=255)
                                                                                                                 aacaaAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAAcc                >  2:3380395/1‑150 (MQ=255)
                                                                                                                 aacaaAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCGTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCTACAGCATCAGCGCCGGGAACGAAcc                >  1:7256902/1‑150 (MQ=255)
                                                                                                                  acaaAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACcg               >  2:4620284/1‑150 (MQ=255)
                                                                                                                    aagaaaGCGCCTCTCGGGGCGCTTTTTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACcgcg             >  2:823257/4‑150 (MQ=255)
                                                                                                                    aagaaaGCGCCTCTCGGGGCGCTTTTTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCGTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACcgcg             >  1:792878/4‑150 (MQ=255)
                                                                                                                    aaaaaaGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCAc                                                                                                  >  2:6953882/1‑65 (MQ=255)
                                                                                                                    aaaaaaGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTTCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGACCGAACcgcg             >  1:1905816/1‑150 (MQ=255)
                                                                                                                    aaaaaaGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACcgcg             >  1:6782911/1‑150 (MQ=255)
                                                                                                                    aaaaaaGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACcgcg             >  2:4728069/1‑150 (MQ=255)
                                                                                                                    aaaaaaGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACcgcg             >  2:3098802/1‑150 (MQ=255)
                                                                                                                    aaaaaaGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACcgcg             >  2:1748296/1‑150 (MQ=255)
                                                                                                                    aaaaaaGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACcgcg             >  1:5355354/1‑150 (MQ=255)
                                                                                                                     agaaaGCGCCTCTCGGGGCGCTTTTTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCgg            >  1:5948287/3‑150 (MQ=255)
                                                                                                                     agaaaGCGCCTCTCGGGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCgg            >  1:444864/3‑150 (MQ=255)
                                                                                                                     aaaaaGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCgg            <  1:5893732/150‑1 (MQ=255)
                                                                                                                      gaaaGCGCCTCTCGGGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAgg                                                            <  1:4235844/100‑1 (MQ=255)
                                                                                                                       aaaGCGCCTCTCGGGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCGGGt          >  2:5344107/1‑150 (MQ=255)
                                                                                                                       aaaGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCGGGt          <  1:6923224/150‑1 (MQ=255)
                                                                                                                        aaGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAgg                                                            >  1:733924/1‑99 (MQ=255)
                                                                                                                        aaGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAgg                                                            <  2:733924/99‑1 (MQ=255)
                                                                                                                          gcgcCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTCGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGGAGACCGGCGACCGCATCAGCGCCGGGACCGAACCGCGGGTCAc       >  2:6651234/1‑150 (MQ=255)
                                                                                                                            gcCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCGGGTCACCa     >  1:1247502/1‑150 (MQ=255)
                                                                                                                              cccGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCGGGTCACCAgg   <  1:1652518/150‑1 (MQ=255)
                                                                                                                              cccGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCGGGTCACCAgg   >  2:3917657/1‑150 (MQ=255)
                                                                                                                               ccGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGACCGAACCGCGGGTCACCAGGc  >  2:2731407/1‑150 (MQ=255)
                                                                                                                                  |                                                                                                                                                  
TGATGTGTTGCCTGTGCGGGCCTCTTCGCGGGTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCGGGTCACCAGGC  >  NC_002947/5764394‑5764670

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: