Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 550,549 C→G 7.3% intergenic (+30/‑112) tufB → / → rpsJ elongation factor Tu‑B/30S ribosomal protein S10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947550,5490CG 7.3% 179.3 / 2.5 64intergenic (+30/‑112)tufB/rpsJelongation factor Tu‑B/30S ribosomal protein S10
Reads supporting (aligned to +/- strand):  ref base C (27/32);  new base G (3/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

ATTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGG  >  NC_002947/550403‑550694
                                                                                                                                                   |                                                                                                                                                  
aTTCAGATTACTGTTACCCTGATCAAGAACATCGAAATGCAAGACGGTCTGCGCTACGCAATCCGTGAAGGCGGTCGTACCCTCGGCGCCGGCCATGTGGCCAAAATCATCGAGTAATCACAAGACCGA‑TTGAAAAAACCCCCCCACAg                                                                                                                                                  >  6:42525/1‑149 (MQ=18)
    aGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCggg                                                                                                                                              <  5:42525/149‑1 (MQ=255)
    aGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCggg                                                                                                                                              <  3:105060/149‑1 (MQ=255)
           tGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑ttttt                                                                                                                                      >  5:186648/1‑149 (MQ=255)
           tGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑ttttt                                                                                                                                      >  3:301799/1‑149 (MQ=255)
           tGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑ttttt                                                                                                                                      >  3:124312/1‑149 (MQ=255)
           tGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑ttttt                                                                                                                                      >  2:175934/1‑149 (MQ=255)
             ttACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGACTTGAAAAAACCCCCGCTCAGCGGGGG‑tttttt                                                                                                                                     <  1:243727/149‑1 (MQ=255)
                cccTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTAt                                                                                                                                 >  6:101046/1‑149 (MQ=255)
                cccTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCCACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑ttttttttta                                                                                                                                 >  3:329768/1‑147 (MQ=255)
                 ccTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTGAGCGGGGG‑ttttttcaag                                                                                                                                 >  7:129664/1‑144 (MQ=255)
                 ccTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTGAGCGGGGG‑ttttttc                                                                                                                                    <  8:129664/145‑2 (MQ=255)
                 ccTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTAtt                                                                                                                                >  2:376572/1‑149 (MQ=255)
                     aTCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGTCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTggg                                                                                                                             <  2:265785/148‑1 (MQ=255)
                     aTCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTgg                                                                                                                              <  8:53908/147‑1 (MQ=255)
                        aaGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGa                                                                                                                         <  2:340543/149‑1 (MQ=255)
                        aaGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGa                                                                                                                         >  7:220719/1‑149 (MQ=255)
                               tCGCAATGGAAGACGTTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACAGTCGGCGCCGGCGTTGTGGCCAAAATCATAGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGATCAGAGGGGG‑TTTTTTTTATTGGTGTGAAaataat                                                                                                                  >  6:98947/1‑149 (MQ=25)
                               tCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAAt                                                                                                                  >  3:363029/1‑149 (MQ=255)
                                  cAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTgg                                                                                                               <  1:175934/149‑1 (MQ=255)
                                    aTGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTgggg                                                                                                             <  8:93488/149‑1 (MQ=255)
                                    aTGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTgggg                                                                                                             <  4:196161/149‑1 (MQ=255)
                                       gAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGc                                                                                                           <  5:41909/148‑1 (MQ=255)
                                       gAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCg                                                                                                          >  7:263288/1‑149 (MQ=255)
                                             ggTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCtata                                                                                                    <  7:336915/149‑1 (MQ=255)
                                             ggTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCtata                                                                                                    <  6:31123/149‑1 (MQ=255)
                                                   cgcTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTg                                                                                                                               >  8:123513/1‑116 (MQ=255)
                                                   cgcTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTg                                                                                                                               <  7:123513/116‑1 (MQ=255)
                                                           tATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTgggg                                                                                                             <  8:335357/126‑1 (MQ=255)
                                                           tATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTgggg                                                                                                             >  7:335357/1‑126 (MQ=255)
                                                           tATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGcctcct                                                                                      >  5:298555/1‑149 (MQ=255)
                                                             tCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGcctcct                                                                                      <  6:186648/147‑1 (MQ=255)
                                                                 tGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGc                                                                                >  4:48404/1‑149 (MQ=255)
                                                                 tGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGc                                                                                >  4:43385/1‑149 (MQ=255)
                                                                   aaGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCg                                                                               <  3:48404/148‑1 (MQ=255)
                                                                   aaGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCg                                                                               <  3:43385/148‑1 (MQ=255)
                                                                    aGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCggg                                                                             <  2:360799/149‑1 (MQ=255)
                                                                          tCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGt                                                                       <  5:335019/149‑1 (MQ=255)
                                                                          tCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGt                                                                       <  2:95748/149‑1 (MQ=255)
                                                                            gTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGc                                                                     <  4:301799/149‑1 (MQ=255)
                                                                            gTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGc                                                                     <  4:124312/149‑1 (MQ=255)
                                                                                     gcgcCGGCGTTGTGGCCAAAATCATCTAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTGGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTgg                                                            <  5:98947/149‑1 (MQ=255)
                                                                                      cgcCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTggg                                                           <  6:283679/149‑1 (MQ=255)
                                                                                                    ccAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAAACCACGCTCAGGGGGGG‑TTTTTTTTTTTGGGTTGACACTTATTGGGGGGGGCGATAGAATAACGCCACCATTTAGCCGTCGTAATGAGTGCCTTGGGAACACCGTGGAGtg                                             >  1:64065/1‑148 (MQ=255)
                                                                                                     cAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCt                                            >  6:260202/1‑149 (MQ=255)
                                                                                                            gtcGAGTGATTACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATc                                      >  7:96176/2‑148 (MQ=255)
                                                                                                               gAGTAATCCCTCGACCGT‑TTGAAAAAACCCCCGCTGAGCGGGGG‑TTTTTTTTAATGGGTTGACACTCATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCaa                                   <  8:295845/148‑1 (MQ=255)
                                                                                                                 gTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAAt                                  <  4:329768/148‑1 (MQ=255)
                                                                                                                 gTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGc                                <  1:376572/149‑1 (MQ=255)
                                                                                                                   aaTCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCa                               <  8:263288/148‑1 (MQ=255)
                                                                                                                     tCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCaaaa                            >  1:264100/1‑149 (MQ=255)
                                                                                                                     tCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCaaa                             >  4:341791/1‑148 (MQ=255)
                                                                                                                         tCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGGACGTTGGGAACAGACTGGAGTCTGAAATCCAATGAAatat                           >  8:253874/1‑146 (MQ=255)
                                                                                                                           gACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAg                        <  4:363029/147‑1 (MQ=255)
                                                                                                                                 ataaaaaaaaCCCCCGCTGAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATcc                 >  6:390782/5‑149 (MQ=255)
                                                                                                                                  ttGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATTCAAAATCAGCAAATCCg                >  8:12576/1‑149 (MQ=255)
                                                                                                                                     aaaaaaCCCCCGCTCAGCGGGGG‑‑TTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTAt             >  1:314316/1‑148 (MQ=255)
                                                                                                                                      aaaaaCCCCCGCTGAGCGGGGG‑TTTTTTTTATTGGGTTCACACTAATTGGGGGCGTCTATAGAATCACGCATCCTTTTAGCGGGAGTAGTGCGTCCGTTGGGAACAGCATGGAGTCTGAAAACCAATGCAAAATCAGCAAATCCGTAt             >  7:134035/1‑148 (MQ=255)
                                                                                                                                      aaaaaCCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTAt             >  7:188355/1‑148 (MQ=255)
                                                                                                                                      aaaaaCCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATc            <  5:101046/149‑1 (MQ=255)
                                                                                                                                           cccccGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGATTCTATATAATCACGCATCCTTTTAGCGGTCGTAGTGCGTCCGTTGTGAACATCCTTGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGtt       <  2:64065/149‑1 (MQ=255)
                                                                                                                                            ccccGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAAt                           >  5:172584/1‑128 (MQ=255)
                                                                                                                                            ccccGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAAt                           <  6:172584/128‑1 (MQ=255)
                                                                                                                                                gCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAgg  <  5:260202/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGA‑TTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGG  >  NC_002947/550403‑550694

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑