Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,836,443 T→C 7.0% intergenic (‑154/‑110) yohC ← / → ttcA inner membrane protein/tRNA 2‑thiocytidine(32) synthetase TtcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,836,4430TC 7.0% 165.1 / 2.3 57intergenic (‑154/‑110)yohC/ttcAinner membrane protein/tRNA 2‑thiocytidine(32) synthetase TtcA
Reads supporting (aligned to +/- strand):  ref base T (30/23);  new base C (2/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

TCACTCCTCTCAATGTTGCCAAACGCCCCCCAATATGGAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATGGGCACCCTCTCGGTCAACCAGAACAAACT  >  NC_002947/1836295‑1836584
                                                                                                                                                     |                                                                                                                                              
tCACTCCTCTCAATGTTGCCAAACGCCCCCCAATATGGAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCt                                                                                                                                                <  1:136484/149‑1 (MQ=255)
    tcctcTCAATGTTGCCAAACGCCCCCCAATATGGAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGt                                                                                                                                              >  8:177054/1‑147 (MQ=255)
       tctcAATGTTGCCAAACGCCCCCCAATATGGAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCt                                                                                                                                                <  4:20474/142‑1 (MQ=255)
       tctcAATGTTGCCAAACGCCCCCCAATATGGAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCt                                                                                                                                                >  3:20474/1‑142 (MQ=255)
           aaTGTTGCCAAACGCCCCCCAATATGGAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGc                                                                                                                                    >  6:223472/1‑149 (MQ=255)
           aaTGTTGCCAAACGCCCCCCAATATGGAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGc                                                                                                                                    <  7:220922/149‑1 (MQ=255)
              gTTGCCAAACGCCCCCCAATATGGAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGt                                                                                                                                 >  8:160717/1‑149 (MQ=255)
               ttGCCAAACGCCCCCCAATATGGAGCGTGAAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGt                                                                                                                                 >  1:22056/1‑149 (MQ=255)
                  ccAAACGCCCCCCAATATGGAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGAt                                                                                                                             <  8:117496/149‑1 (MQ=255)
                           ccccAATATGGAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGaa                                                                                                                    <  7:160717/149‑1 (MQ=255)
                            cccAATATGGAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCTAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGAt                                                                                                                             >  5:216342/1‑139 (MQ=255)
                            cccAATATGGAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGAt                                                                                                                             <  6:216342/139‑1 (MQ=255)
                            cccAATATGGAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAg                                                                                                                   >  8:21282/1‑149 (MQ=255)
                             ccAATATGGAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGt                                                                                                                  >  6:198795/1‑149 (MQ=255)
                                     gAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCg                                                                                                          <  5:223472/149‑1 (MQ=255)
                                        cGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGaa                                                                                                       >  2:237170/1‑149 (MQ=255)
                                             aGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGc                                                                                                    <  5:198795/148‑1 (MQ=255)
                                                 gCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTa                                                                                                 >  7:350252/1‑147 (MQ=255)
                                                    cGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAgcg                                                                                            <  7:177054/149‑1 (MQ=255)
                                                       ggCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGaa                                                                                                                    >  1:137123/1‑122 (MQ=255)
                                                       ggCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGaa                                                                                                                    <  2:137123/122‑1 (MQ=255)
                                                          ccTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTcg                                                                                      >  3:187820/1‑149 (MQ=255)
                                                          ccTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTcg                                                                                      >  1:134917/1‑149 (MQ=255)
                                                                        aaGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTaa                                                                         >  3:191648/1‑148 (MQ=255)
                                                                         aGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGc                                                                             >  3:357245/1‑143 (MQ=255)
                                                                         aGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGc                                                                             <  4:357245/143‑1 (MQ=255)
                                                                         aGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCCGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTaaaa                                                                       >  4:355150/1‑149 (MQ=255)
                                                                         aGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCCGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTaaaa                                                                       >  4:360733/1‑149 (MQ=255)
                                                                             gTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGc                                                                   <  3:57637/149‑1 (MQ=255)
                                                                              tAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGcc                                                                  <  7:21282/149‑1 (MQ=255)
                                                                              tAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCAGTGCTTTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGAGTCGCGTTTGTGCGTAAAATAGcc                                                                  >  5:290079/1‑149 (MQ=255)
                                                                                  gtgGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCt                                                              >  2:155572/1‑149 (MQ=255)
                                                                                             tCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACAcc                                                   >  1:133138/1‑149 (MQ=255)
                                                                                                 aGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTcgc                                               <  1:268554/149‑1 (MQ=255)
                                                                                                   tataGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGGTGAATACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCCTAAAATAGCCGGCTTTTGTCACACCtagag                                              >  5:291279/1‑144 (MQ=255)
                                                                                                       ggAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATcc                                         <  1:155572/149‑1 (MQ=255)
                                                                                                            ccGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCAAACCTCGCGGATCCAAgcg                                    >  4:64490/1‑149 (MQ=255)
                                                                                                                 aCCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTAt                                <  4:191648/148‑1 (MQ=255)
                                                                                                                 aCCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCCGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTAt                                <  3:355150/148‑1 (MQ=255)
                                                                                                                 aCCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCCGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTAt                                <  3:360733/148‑1 (MQ=255)
                                                                                                                  ccTACATAAACCGCTCTTTTTAGAGCAAATTTCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAATTAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATg                               <  6:291279/148‑1 (MQ=255)
                                                                                                                  ccTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATgg                              <  3:64490/149‑1 (MQ=255)
                                                                                                                   cTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATggg                             >  5:101751/1‑149 (MQ=255)
                                                                                                                   cTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATgg                              <  2:113844/148‑1 (MQ=255)
                                                                                                                   cTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATgg                              >  1:113844/1‑148 (MQ=255)
                                                                                                                     aCATAAACCGCGCTTTTTAGAGCAAATTGCGCTGTGAAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATggg                             >  1:54908/1‑148 (MQ=255)
                                                                                                                     aCATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATGGTCa                           >  5:276124/1‑149 (MQ=255)
                                                                                                                      gaTAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGTCGTGCTGTTGATTACATTGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATGGGCa                           <  6:290079/147‑1 (MQ=255)
                                                                                                                      cATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATGGGCAc                          >  5:158188/1‑149 (MQ=255)
                                                                                                                      cATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGATATGGGCAc                          >  8:343784/1‑149 (MQ=255)
                                                                                                                       aTAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATGGGCAcc                         <  1:237170/149‑1 (MQ=255)
                                                                                                                        tAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATGGGCAccc                        >  7:196911/1‑149 (MQ=255)
                                                                                                                              gcgcTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATGGGCACCCTCTCgg                  >  8:326995/1‑149 (MQ=255)
                                                                                                                                    ctcacagCAAATTGCGCTGT‑CAGCCGCGCTGTTGATTACAATGAAGTCTTTCTCGGACTGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATGGGCACCCTCTCGGTCaa              <  8:350252/142‑1 (MQ=255)
                                                                                                                                        gagCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATGGGCACCCTCTCGGTCAACCAGaa        <  2:133138/149‑1 (MQ=255)
                                                                                                                                         agCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATGGGCACCCTCTCGGTCAAc             <  6:351265/143‑1 (MQ=255)
                                                                                                                                         agCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATGGGCACCCTCTCGGTCAAc             >  5:351265/1‑143 (MQ=255)
                                                                                                                                         agCAAATTGCGCTGT‑AAGCCGTGCTGTTGAATACAATGAAGTCTATATCGGAATGCCTACAGCGGCGTCGCCTTTGTGCGTAAAATAGACGGCTTTTGACACATCTCGCGGATCCAAGCGCTATGGGCACCCTCACGGTCAAACAGaa        >  3:265931/1‑148 (MQ=255)
                                                                                                                                           cAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATGGGCACCCTCTCGGTCAACCAGAACa      >  5:124362/1‑148 (MQ=255)
                                                                                                                                              aTTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATGGGCACCCTCTCGGTCAACCAGAACAAACt  <  2:134917/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
TCACTCCTCTCAATGTTGCCAAACGCCCCCCAATATGGAGCGT‑AAGCGGCCCGAGGCCCTAACACCGATGCAAGCGGTAATGTGGCCTTATGTCGCAGTATAGGAAGCCGTTACCTACATAAACCGCGCTTTTTAGAGCAAATTGCGCTGT‑AAGCCGTGCTGTTGATTACAATGAAGTCTTTATCGGAATGCCTACAGCGGCGTCGCGTTTGTGCGTAAAATAGCCGGCTTTTGTCACACCTCGCGGATCCAAGCGCTATGGGCACCCTCTCGGTCAACCAGAACAAACT  >  NC_002947/1836295‑1836584

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