Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,110,924 T→A 7.2% V317D (GTC→GAC)  PP_0970 → cyanate MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,110,9240TA 7.2% 149.5 / 3.1 56V317D (GTC→GAC) PP_0970cyanate MFS transporter
Reads supporting (aligned to +/- strand):  ref base T (27/25);  new base A (2/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-01

ATCGATCCTGATCGGCCGCGGCCTCACCCCTACCGAGGCAGGTTTGGTGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCTTTGCCCTGGCGCTGACCTTGATCGTGCTGCGCTCGA  >  NC_002947/1110781‑1111061
                                                                                                                                                   |                                                                                                                                         
aTCGATCCTGATCGGCCGCGGCCTCACCCCTACCGAGGCAGGTTTGGTGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGt                                                                                                                                      >  7:299165/1‑149 (MQ=255)
        tGATCGGCCGCGGCCTCACCCCTACCGAGGCAGGTTTGGTGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACg                                                                                                                              >  3:233766/1‑149 (MQ=255)
         gATCGGCCGCGGCCTCACCCCTACCGAGGCAGGTTTGGTGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGc                                                                                                                             >  7:137690/1‑149 (MQ=255)
                 gcgGCCTCACCCCTACCGAGGCAGGTTTGGTGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCgggg                                                                                                                     <  4:233766/149‑1 (MQ=255)
                 gcgGCCTCACCCCTACCGAGGCAGGTTTGGTGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCgggg                                                                                                                     <  3:325107/149‑1 (MQ=255)
                             cttccGAGGCAGGTTTGGTGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTTGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGc                                                                                                         <  8:137690/146‑1 (MQ=255)
                                ccGAGGCAGGTTTGGTGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCt                                                                                                       <  5:236740/148‑1 (MQ=255)
                                     gCAGGTTTGGTGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATgc                                                                                                 >  1:274607/1‑149 (MQ=255)
                                      cAGGTTTGGTGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATgcg                                                                                                >  6:279369/1‑149 (MQ=255)
                                       aGGTTTGGTGCTGTCAGGGTTGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATgcgc                                                                                               <  6:123351/149‑1 (MQ=255)
                                          tttGGTGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGt                                                                                                     >  1:136903/1‑140 (MQ=255)
                                          tttGGTGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGt                                                                                                     <  2:136903/140‑1 (MQ=255)
                                              gTGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCt                                                                                                       <  1:16338/134‑1 (MQ=255)
                                              gTGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCt                                                                                                       >  2:16338/1‑134 (MQ=255)
                                               tGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTcc                                                                                       >  5:209404/1‑149 (MQ=255)
                                                 cTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCgg                                                                                     <  8:328521/149‑1 (MQ=255)
                                                      aGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTtgt                                                                                >  6:136184/1‑149 (MQ=255)
                                                        ggTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCgggg                                                                                                                     >  1:242860/1‑110 (MQ=255)
                                                        ggTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCgggg                                                                                                                     <  2:242860/110‑1 (MQ=255)
                                                        ggTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGt                                                                                                                 <  7:295420/114‑1 (MQ=255)
                                                        ggTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGt                                                                                                                 >  8:295420/1‑114 (MQ=255)
                                                        ggTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGACATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTtgtgc                                                                              >  7:333770/1‑148 (MQ=255)
                                                        ggTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGACATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTgg                                                                              >  7:102173/1‑149 (MQ=255)
                                                             gTAATCGTGCAGTTGATCAGTTCGCTAAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTg                                                                         <  4:82826/149‑1 (MQ=255)
                                                               aaTCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGcc                                                                                                        <  4:103532/116‑1 (MQ=255)
                                                               aaTCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGcc                                                                                                        >  3:103532/1‑116 (MQ=255)
                                                                  atggCACTTGATCAGTTCGCTCAGCGCCCCCTCGCGGGCCACCCCCCCCAAGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTgg                                                                        <  1:298569/146‑1 (MQ=255)
                                                                    tGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCAtt                                                                  <  2:274607/149‑1 (MQ=255)
                                                                       aGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGACATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATtgtg                                                               <  8:102173/149‑1 (MQ=255)
                                                                       aGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGACATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGCGGGTGGGCCATtgtg                                                               <  8:333770/149‑1 (MQ=255)
                                                                                    gCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTgg                                                     <  3:200273/146‑1 (MQ=255)
                                                                                    gCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTgg                                                     >  4:200273/1‑146 (MQ=255)
                                                                                      tCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGCCCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGGTGTGGGGGTGGGCCATTGTGCTgg                                                           >  5:77864/1‑138 (MQ=255)
                                                                                      tCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCa                                                                    >  6:302752/1‑129 (MQ=255)
                                                                                      tCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCa                                                                    <  5:302752/129‑1 (MQ=255)
                                                                                      tCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTgg                                                           <  6:77864/138‑1 (MQ=255)
                                                                                      tCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAg                                                >  5:22654/1‑149 (MQ=255)
                                                                                         gcgcCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGc                                             >  1:232169/1‑149 (MQ=255)
                                                                                            cccccTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGt                                          <  6:22654/149‑1 (MQ=255)
                                                                                                ctggctggCCCCCCGC‑‑‑‑GGCAAGGACCCGCGGCTTGCAATCGTGGGGGTCATGTCGGGCCCGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCt                                      <  8:187076/149‑1 (MQ=255)
                                                                                                        ccACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCTTTGCCCTg                              >  2:45241/1‑149 (MQ=255)
                                                                                                        ccACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCTTTGCCCTg                              <  8:299165/149‑1 (MQ=255)
                                                                                                         cACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCGGGGGCAGGGCGGTACCTTTGCCCTgg                             <  5:136184/149‑1 (MQ=255)
                                                                                                                      cAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTGCCTTTGCCCTGGCGCTGACCt                    >  4:93971/1‑149 (MQ=255)
                                                                                                                         ggACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCTTTGCCCTGGCGCTGACCTTGa                 >  4:153412/1‑149 (MQ=255)
                                                                                                                          gACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCTTTGCCCTGGCGCTGACCTTGAt                >  3:150320/1‑149 (MQ=255)
                                                                                                                            ccAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCttt                                    <  4:98787/127‑1 (MQ=255)
                                                                                                                            ccAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCttt                                    >  3:98787/1‑127 (MQ=255)
                                                                                                                                cGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCTTTGCCCTGGCGCTGACCTTGATCGtgct          <  6:209404/149‑1 (MQ=255)
                                                                                                                                cGGCTGGCAATCGTGGGGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGGGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCTTTGCCCTGGCGCTGACCTTGATCGtgct          <  5:279369/149‑1 (MQ=255)
                                                                                                                                    tGGCAATCGTGGTGGTCATGTCGGGCACGCTGTCGGGGCTGATCTGTTGCCAGTATGCGCCACTTTCATGTTTGTGGTGGTGGTACATTGTGCTGTGCATGGGGCAGGGCGGTACCTTTCACCTGGCGCTGACCTAGATCGTGCTgcg       >  6:369618/1‑148 (MQ=255)
                                                                                                                                        aaTCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGTCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCTTTGCCCTGGAGCTGACCTTGATCGTGCTGCGCTCg   >  8:94699/1‑148 (MQ=255)
                                                                                                                                        aaTCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCTTTGccc                                <  3:174423/119‑1 (MQ=255)
                                                                                                                                        aaTCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCTTTGccc                                >  4:174423/1‑119 (MQ=255)
                                                                                                                                        aaTCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCTTTGCCCTGGCGCTGACCTTGATCGTGCTGCGCTCg   >  3:225494/1‑148 (MQ=255)
                                                                                                                                        aaTCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCTTTGCCCTGGCGCTGACCTTGATCGTGCTGCGCTCGa  >  8:328197/1‑149 (MQ=255)
                                                                                                                                                ttgtCATGTCGAGCACGCTGGCGGAGATGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTATGGGGGTGGGCCATTGTGATGGGCCTGGGGCAGGGCGGTACCTTTGCAATGGCGCTGACCTTGAACGTGCTGAGCTc    >  8:225142/3‑139 (MQ=255)
                                                                                                                                                tggtCATGTGGGGCACGCTGGCGGTGCTGTTCGGTTGCCTGTATGCGCAACTTTCCGGTTTGTGGTGGTGGGCTATTGTGTTGGGTCTGGGGCAGGGCGGTACCTTTGCCCTGGCGCTGACCTTGATCGTGCTGCGCTc    <  7:225142/139‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
ATCGATCCTGATCGGCCGCGGCCTCACCCCTACCGAGGCAGGTTTGGTGCTGTCAGGGTCGGTAATCGTGCAGTTGATCAGTTCGCTCAGCGCCCCCTGGCTGGCCACCCGC‑‑‑‑GGCAAGGACCAGCGGCTGGCAATCGTGGTGGTCATGTCGGGCACGCTGGCGGGGCTGTTCGGTTGCCTGTATGCGCCACTTTCCGGGTTGTGGGGGTGGGCCATTGTGCTGGGCCTGGGGCAGGGCGGTACCTTTGCCCTGGCGCTGACCTTGATCGTGCTGCGCTCGA  >  NC_002947/1110781‑1111061

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