Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,583,529 T→A 12.7% K735* (AAG→TAG) ‡ parC ← DNA topoisomerase IV subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,583,5290TA12.7% 155.8 / 7.7 60K735* (AAG→TAG) ‡parCDNA topoisomerase IV subunit A
Reads supporting (aligned to +/- strand):  ref base T (26/26);  new base A (4/4);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCATATCATCCGCGAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGATCGGTAACAAATTCTTCACGGCTGG  >  NC_002947/5583387‑5583669
                                                                                                                                              |                                                                                                                                            
ggCATATCATCCGCGAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGcc                                                                                                                                        >  3:123622/1‑149 (MQ=255)
 gCATATCATCCGCGAAATAACCGTCCGACTCCAGCATTTCTGCCCGAACCGGCGTTGAAATTGCCGTGCCAGACATTACAACCGGTTATGTCGGGACTTCCGCCGGCAACCCCTCGACGCGCGGGAACCCGCGTGGCAGCTTGCTGCcg                                                                                                                                       >  6:189129/1‑149 (MQ=255)
     atcatcCGCGAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGcc                                                                                                                                   <  4:18192/149‑1 (MQ=255)
         tcCGCGAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGAcc                                                                                                                               >  5:187759/1‑149 (MQ=255)
           cgcgAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGTGCGGCAGCAAGCTGCCacgcgggttcc                                                                                                                             <  4:273707/149‑8 (MQ=255)
           cgcgAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTtgctg                                                                                                                                          <  8:12068/136‑1 (MQ=255)
           cgcgAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTtgctg                                                                                                                                          >  7:12068/1‑136 (MQ=255)
            gagAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTTACGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGAAGGAACCCGCGTGGCAGCTTGCTGCCCCGACGCCcgcg                                                                                                                            >  4:57497/3‑149 (MQ=255)
              gAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCt                                                                                                                          >  5:323671/1‑149 (MQ=255)
                        tcctaatCCAGGATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTa                                                                                                                <  3:57497/143‑1 (MQ=255)
                        tCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTg                                                                                                                  >  8:105246/1‑147 (MQ=255)
                        tCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTg                                                                                                                  >  5:99336/1‑147 (MQ=255)
                            aCTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGc                                                                                                            <  1:130931/149‑1 (MQ=255)
                            aCTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGc                                                                                                            <  4:123622/149‑1 (MQ=255)
                               ccAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTc                                                                                                          >  5:164740/1‑148 (MQ=255)
                                cAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTc                                                                                                          <  8:61676/147‑1 (MQ=255)
                                 aGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGCGGCAGCAAGCTGCCacgcg                                                                                                                                   <  5:188557/121‑2 (MQ=255)
                                 aGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGCGGCAGCAAGCTGCCacgcg                                                                                                                                   >  6:188557/1‑120 (MQ=255)
                                 aGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGCGGCAGCAAGCTGCCacgcg                                                                                                                                   >  7:137772/1‑120 (MQ=255)
                                 aGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGCGGCAGCAAGCTGCCacgcg                                                                                                                                   <  8:137772/121‑2 (MQ=255)
                                 aGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGCGGCAGCAAGCTGCCacgcg                                                                                                                                   <  1:177682/121‑2 (MQ=255)
                                 aGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGCGGCAGCAAGCTGCCacgcg                                                                                                                                   >  2:177682/1‑120 (MQ=255)
                                 aGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACATGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGATTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAg                                                                                                       >  7:173212/1‑149 (MQ=255)
                                   cATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGgt                                                                                                     >  4:56408/1‑149 (MQ=255)
                                       tCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGgtcgt                                                                                                  <  3:56408/148‑1 (MQ=255)
                                           cccGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGcc                                                                                             <  6:99336/149‑1 (MQ=255)
                                           cccGAAACGGCGTTGAAATTGCAGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTcc                                                                                                <  2:86931/146‑1 (MQ=255)
                                            ccGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCt                                                                                            >  6:96887/1‑149 (MQ=255)
                                            ccGAAACGGCGTTGAAATTGCCGTTCCACACATTACAACCCGTTATGTCGGCACTTCCACCTCCAAACCATCGACGCGCAGGAACCCGCGTGGCAGCTAGCTGCCGCGCCGACCGCGCTCTCCCTTGTAGTGCTCCAGGTCGTCCCCCt                                                                                            >  3:185953/1‑149 (MQ=255)
                                             cGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCtt                                                                                           <  5:96887/149‑1 (MQ=255)
                                              gAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTc                                                                                          <  7:219671/149‑1 (MQ=255)
                                               aaaCGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGCTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCa                                                                                         <  5:326414/149‑1 (MQ=255)
                                                 aCGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCaaa                                                                                       <  6:164740/149‑1 (MQ=255)
                                                      gTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTc                                                                                          <  6:309768/141‑1 (MQ=255)
                                                      gTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTc                                                                                          >  5:309768/1‑141 (MQ=255)
                                                                   ttCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGAAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCt                                                                                                           >  3:32335/1‑111 (MQ=255)
                                                                       aGACATTACAACCCGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTCCGCTTGCCGGCCt                                                                 >  8:253214/1‑149 (MQ=255)
                                                                        gACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTg                                                                >  7:239912/1‑149 (MQ=255)
                                                                        cgacctACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTg                                                                <  3:270315/144‑1 (MQ=255)
                                                                           aTTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCaa                                                             >  5:209816/1‑149 (MQ=255)
                                                                           aTTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCaa                                                             <  8:239912/149‑1 (MQ=255)
                                                                               cAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTg                                                                >  4:270315/1‑142 (MQ=255)
                                                                                aaCCGGTTATGTCGGCACTTCCACCTGCACCCCCTCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACaa                                                        <  4:44608/149‑1 (MQ=255)
                                                                                 aCCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAg                                                       <  2:182060/149‑1 (MQ=255)
                                                                                  cgggTTATGTCGGGACTTCCAGCTGCTACCCATCGACGCGCTGGACCCCCCGTTGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTTCGCTTGCCGGCCTGCAATACAAgg                                                      <  4:185953/147‑1 (MQ=255)
                                                                                     gTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTc                                                   >  6:200581/1‑149 (MQ=255)
                                                                                     gTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTc                                                   <  6:187759/149‑1 (MQ=255)
                                                                                           tCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACTAGGGTCGCGCCt                                             >  5:32506/1‑149 (MQ=255)
                                                                                           tCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCt                                             <  6:323671/149‑1 (MQ=255)
                                                                                            cGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCtt                                            <  7:105246/149‑1 (MQ=255)
                                                                                              gCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCg                                          <  6:209816/149‑1 (MQ=255)
                                                                                                      aCCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCAc                                  <  5:200581/149‑1 (MQ=255)
                                                                                                          gCAACCCATCGACGCGCTCGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCAAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGTGTCGCGCCTTCGGGAATCACGGc                               >  6:317560/1‑148 (MQ=255)
                                                                                                             aCCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTg                                                                <  2:36009/112‑1 (MQ=255)
                                                                                                             aCCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTg                                                                >  1:36009/1‑112 (MQ=255)
                                                                                                              cccATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGAt                          <  8:173212/149‑1 (MQ=255)
                                                                                                                cATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGc                               >  3:53089/1‑142 (MQ=255)
                                                                                                                              aaCCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGATCGGTAACAAAttctt           >  1:80285/1‑148 (MQ=255)
                                                                                                                                 ccGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGATCGGTAACAAATTCTTCACg       >  2:282690/1‑149 (MQ=255)
                                                                                                                                  cgcgTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGATCGGTAACAAATTCTTCACgg      >  3:91790/1‑149 (MQ=255)
                                                                                                                                   gcgTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGATCGGTAACAAATTCTTCACGGc     >  7:170924/1‑149 (MQ=255)
                                                                                                                                      tGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGATCGGTAACAAATTCTTCACGGCTgg  <  4:53089/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                            
GGCATATCATCCGCGAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGATCGGTAACAAATTCTTCACGGCTGG  >  NC_002947/5583387‑5583669

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: ‑