Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,240,008 T→C 7.3% N114S (AAC→AGC)  mglB ← methyl‑galactoside transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,240,0080TC 7.3% 160.7 / 3.4 56N114S (AAC→AGC) mglBmethyl‑galactoside transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (28/23);  new base C (2/2);  total (31/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 = 8.33e-01

GACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTGGTCATTCT  >  NC_000913/2239862‑2240155
                                                                                                                                                  |                                                                                                                                                   
gACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTttgtt                                                                                                                                                    >  8:350693/1‑148 (MQ=255)
gACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTttgtt                                                                                                                                                    <  7:350693/148‑1 (MQ=255)
  ccGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTgaa                                                                                                                                                 <  1:212062/149‑1 (MQ=255)
     tCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTgaaga                                                                                                                                               >  5:20649/1‑148 (MQ=255)
         tGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAac                                                                                                                                           <  6:20649/148‑1 (MQ=255)
             cAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCg                                                                                                                                      <  5:256461/148‑1 (MQ=255)
                               cgccgcCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCAcgcg                                                                                                                    >  8:214652/1‑149 (MQ=255)
                                 ccgccCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCg                                                                                                                                      >  3:347264/1‑129 (MQ=255)
                                 ccgccCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCg                                                                                                                                      <  4:347264/129‑1 (MQ=255)
                                    cccAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTtc                                                                                                               <  7:214652/149‑1 (MQ=255)
                                     ccAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGCTGAAGAAAac                                                                                                                                           <  7:24966/120‑1 (MQ=255)
                                     ccAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGCTGAAGAAAac                                                                                                                                           <  7:30648/120‑1 (MQ=255)
                                     ccAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGCTGAAGAAAac                                                                                                                                           >  8:30648/1‑120 (MQ=255)
                                     ccAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGCTGAAGAAAac                                                                                                                                           >  8:24966/1‑120 (MQ=255)
                                      cAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTtct                                                                                                              >  5:153743/1‑148 (MQ=255)
                                                 aaTCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTAc                                                                                                  >  3:46108/1‑149 (MQ=255)
                                                     aaaTCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCg                                                                                               >  8:67572/1‑148 (MQ=255)
                                                            cTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGcc                                                                                       >  1:6033/1‑149 (MQ=255)
                                                            cTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCAATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGcc                                                                                       >  2:108219/1‑149 (MQ=255)
                                                                aataatGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGt                                                                                   <  2:6033/149‑1 (MQ=255)
                                                                aataatGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGt                                                                                   >  7:189623/1‑149 (MQ=255)
                                                                 ataatGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTc                                                                                  <  4:46108/149‑1 (MQ=255)
                                                                             tctTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCTCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGAt                                                                      >  2:254402/1‑149 (MQ=255)
                                                                               tttGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGAtt                                                                     >  5:134284/1‑148 (MQ=255)
                                                                               tttGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTg                                                                    >  2:313601/1‑149 (MQ=255)
                                                                                 tGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGAtt                                                                     <  8:189623/146‑1 (MQ=255)
                                                                                 tGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGAtt                                                                     <  6:134284/146‑1 (MQ=255)
                                                                                         gTGCCAACGTAGTAGGCTTTCTCGTAGCAATCCAGCGCCTTACTAGACTGTTATTTGATGAATAAAACCACCGGCACATTTTTCACAAACGCTATCTCAATCACCGTAACCGCAGCTGCCGGGTCAACAAGGTTGATTGCCATTGCCt                                                           >  8:375707/1‑148 (MQ=255)
                                                                                            ccAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCa                                                        <  1:313601/148‑1 (MQ=255)
                                                                                             cAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCAcc                                                      >  2:297774/1‑149 (MQ=255)
                                                                                              aaCGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCTGCACGTTTTGCCAACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTTCAAGTGCCTTCAcca                                                     >  6:258447/1‑148 (MQ=255)
                                                                                              aaCGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCAccc                                                     >  4:284241/1‑149 (MQ=255)
                                                                                              aaCGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTGGATTGCCAGTGCCTTCAccc                                                     >  5:8119/1‑149 (MQ=255)
                                                                                                     tagGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACTGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTGAACCCCTTTCGc                                              >  3:287835/1‑149 (MQ=255)
                                                                                                     tagGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCa                                                                 <  8:257442/130‑1 (MQ=255)
                                                                                                     tagGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCa                                                                 >  7:257442/1‑130 (MQ=255)
                                                                                                       ggCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCTCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCa                                            <  1:254402/149‑1 (MQ=255)
                                                                                                        gCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAg                                           <  7:375707/149‑1 (MQ=255)
                                                                                                          tttGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGc                                          >  6:147661/1‑148 (MQ=255)
                                                                                                             gTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGAtt                                                                     >  8:363808/1‑118 (MQ=255)
                                                                                                             gTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGAtt                                                                     <  7:363808/118‑1 (MQ=255)
                                                                                                              tCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATAc                                     >  8:336253/1‑149 (MQ=255)
                                                                                                                  aGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCg                                                                                      <  5:319250/96‑1 (MQ=255)
                                                                                                                  aGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCg                                                                                      >  6:319250/1‑96 (MQ=255)
                                                                                                                     tATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCg                                                                                                     <  3:99672/78‑1 (MQ=255)
                                                                                                                     tATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCg                                                                                                     >  4:99672/1‑78 (MQ=255)
                                                                                                                          aGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATc                         <  6:8119/149‑1 (MQ=255)
                                                                                                                               cTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGa                           <  1:14029/142‑1 (MQ=255)
                                                                                                                               cTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGa                           >  2:14029/1‑142 (MQ=255)
                                                                                                                                ttACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAg                                           >  5:45459/1‑125 (MQ=255)
                                                                                                                                ttACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAg                                           <  6:45459/125‑1 (MQ=255)
                                                                                                                                   cGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCtt                <  7:336253/149‑1 (MQ=255)
                                                                                                                                    gagaCGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTg               >  2:378884/1‑149 (MQ=255)
                                                                                                                                       aCGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGAc            <  7:351716/149‑1 (MQ=255)
                                                                                                                                        cGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACt           >  3:222192/1‑149 (MQ=255)
                                                                                                                                           ttCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTg          >  5:208288/1‑147 (MQ=255)
                                                                                                                                           ttCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTg          <  6:208288/147‑1 (MQ=255)
                                                                                                                                           ttCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTTGt        <  1:297774/149‑1 (MQ=255)
                                                                                                                                              gttgttgAAGAAACCCACCGGCACGGTTTGCCCACGCGCTTTCTCATTCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGCTTGCCAGTGCCTTCACGCCTTTCGCCAGCAATCCGGCGATCTGCTCGTGCTGCTGGGACTGGGCAt     >  8:207693/2‑149 (MQ=255)
                                                                                                                                                tgttgAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTGGTCAtt    <  5:147661/148‑1 (MQ=255)
                                                                                                                                                tgttgAAGAAAACCACCGGCAAGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTGGTCATTc   >  7:51236/1‑149 (MQ=255)
                                                                                                                                                  ttgAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTg          >  5:19838/1‑140 (MQ=255)
                                                                                                                                                  ttgAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTg          <  6:19838/140‑1 (MQ=255)
                                                                                                                                                  ttgAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTGGTCATTCt  <  1:108219/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTGGTCATTCT  >  NC_000913/2239862‑2240155

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