Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,939,297 G→T G137V (GGG→GTG)  fucU → L‑fucose mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,939,2970GT94.6% 104.2 / ‑6.5 38G137V (GGG→GTG) fucUL‑fucose mutarotase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base T (14/21);  minor base A (2/0);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGATTATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGCGCAATCATGGACGCGCTAAATATATCC  >  NC_000913/2939052‑2939542
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ggCGATTATCCCGTTATTTGAAATGGACAGTTATGCACCGCAGCTGGTGATGATGGCGGCGGTAGAAGATGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGAGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGGACGGGAACATTCTTTTAAAAAAAGTGGTa                                                                                                                                                                                                                                                   >  2:345135/1‑250 (MQ=255)
  cGATTATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAac                                                                                                                                                                                                                                                 >  2:129787/1‑250 (MQ=255)
          ccgttattgtaaattgACAGTTACGCTCCGACGCTGGTGATGATGGCGGCGGTAGAAGGTGATACTCTCGATCCGGAAGTAGAACGACGTTATCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGACGTACGGGAATATGCTTTTAAAAAAAGTGTTAACACCTTAAtc                                                                                                                                                                                                                                        <  1:408587/236‑1 (MQ=255)
           cGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGGCATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAAtct                                                                                                                                                                                                                                       >  1:312939/1‑251 (MQ=255)
                  tGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCg                                                                                                                                                                                                                                >  1:23847/1‑251 (MQ=255)
                  tGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGGGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCg                                                                                                                                                                                                                                >  1:354985/1‑251 (MQ=255)
                         gACAGTTATGCCCCGCCGCTGGTGAGGATGGCGGCGGTAGAAGGTGACACTCTCGATCCGGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTGTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGcc                                                                                                                                                                                                                         <  1:163502/251‑1 (MQ=255)
                             gTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCa                                                                                                                                                                                                                     <  2:23847/251‑1 (MQ=255)
                                           cTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGtt                                                                                                                                                                                                       >  1:107073/1‑251 (MQ=255)
                                           cTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGtt                                                                                                                                                                                                       >  1:331559/1‑251 (MQ=255)
                                              gggatgatGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGATATc                                                                                                                                                                                                    <  2:405502/249‑1 (MQ=255)
                                                   gatgGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACAGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAat                                                                                                                                                                                                <  1:218445/250‑1 (MQ=255)
                                                                        cacTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTAACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAg                                                                                                                                                                          <  2:472316/251‑1 (MQ=255)
                                                                                                acgTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCg                                                                                                                                                  <  1:84542/251‑1 (MQ=255)
                                                                                                     atcGTAATGCGCTTTCACGACAAGCCCCGTGTCCTGACATCATCCGCATCAATCTTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACct                                                                                                                                              <  1:394692/248‑1 (MQ=255)
                                                                                                           aTGCGCTTTCACTACAAGCCAAGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGAGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGCCGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCGAGCGATAGTCGACctgctggtg                                                                                                                                       >  2:366613/1‑248 (MQ=255)
                                                                                                                 tttCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACg                                                                                                                                                                                                                           >  1:143591/1‑161 (MQ=255)
                                                                                                                 tttCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACg                                                                                                                                                                                                                           <  2:143591/161‑1 (MQ=255)
                                                                                                                     ccTACAAGCCCCGGGTCCCGACCTCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTCCGGTAATATTCTTTTAAAAAAAGTGGTAACTCCGTAAGCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAATACATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAGCCATTCCa                                                                                                                             <  2:187486/250‑1 (MQ=255)
                                                                                                                         cAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGccc                                                                                                                         >  2:394028/1‑250 (MQ=255)
                                                                                                                                 gtgtCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACgg                                                                                                                 <  1:259609/251‑1 (MQ=255)
                                                                                                                                                      cAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTa                                                                                             <  1:129787/250‑1 (MQ=255)
                                                                                                                                                         tCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGAAAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTaacgg                                                                                         >  1:124293/1‑248 (MQ=255)
                                                                                                                                                              ttGCGATTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGAACGGGAATATTCTTTTAAAAAAAGAGGTAACACCGTAATCTCATACCGGTACGCCCGAATGACGCGGGCGGTTATCGCATGATGGGGTGAACAATATGAAAGCGGCCCGCCAGCAAGCGATAGTCGACCTGCTGCTACAACAAACCAGCCTGACCACGGAAACTCTCTCTGAACAGCTAAAGGTc                                                                                       >  2:482711/1‑248 (MQ=255)
                                                                                                                                                                    tttATGAACGGGCGCAAAAAGCCTTTGCTATCGTTATCACAGGCGAACGAGCGAAGTACGGTAATATTCTTTTAAAAAAAGTGGTAACGCCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGTGGTGAAAAATATGAAAGCTGCACGCCAGCAATCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGa                                                                               <  1:258561/250‑1 (MQ=255)
                                                                                                                                                                          aaCGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCAtt                                                                        <  2:124293/251‑1 (MQ=255)
                                                                                                                                                                          aaCGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGAGGTAACACCGTAATCTCATACCGGTCCGACCGCATGACGAGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGGACCATAACAGACTGACCACGGACGCTCTCTCTGATCAGCTAAAGGTCAGTAAAGAAACCAtt                                                                        >  1:503349/1‑251 (MQ=255)
                                                                                                                                                                            cGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATtcg                                                                      >  2:313326/1‑251 (MQ=255)
                                                                                                                                                                            cGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTAACGAATGATGGGGTGAAAAATATGAAATCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACTGAAGCTCTCTCTGATCAGCTAAAGGTCAGTAAAGAAACCATtcg                                                                      <  2:482481/251‑1 (MQ=255)
                                                                                                                                                                              ggCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATtcgtc                                                                    <  1:169726/251‑1 (MQ=255)
                                                                                                                                                                                                    gTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACCCCGTAATCTCATACCGATACGCCCTCATGACGCGGGCTGTTATCGATGGAGGTGGTGAAAAAGATGAAAGCGGCACGCCATCAAGCGTTAGTCGTCCTGCTGCTGAACCATACCAGCCTTACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCATTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGAc                                              <  1:457440/251‑1 (MQ=255)
                                                                                                                                                                                                               cGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAAACATACCAGCCTGACCACGGAAGCTCTCTCTGAACAACTAAAGGCCAGTAAAGAAACCATTCGTAGCGATTTCAATGAATTAAAGACGCAGGGaaaa                                    >  1:10734/1‑250 (MQ=255)
                                                                                                                                                                                                                 aaCGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAAtt                                 <  1:349312/251‑1 (MQ=255)
                                                                                                                                                                                                                  aCGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAGTGAATTACAGACGCAGGGTAAAATTc                                <  1:89217/251‑1 (MQ=255)
                                                                                                                                                                                                                       cGAAGTACGGGCATATTCCTTTACAAAAAGTGGTCACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAAAATGAACGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACAACGGAAACTCTCTCTGAAAAGCTAAAGGTCAGTAAAGAAACCAATCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTgcgc                           >  2:79663/1‑251 (MQ=255)
                                                                                                                                                                                                                                ggAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATTCCAGCCTGACCACGGAAGCTCTCTCTGAACATCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGCGCAATCATGGa                  <  1:310620/251‑1 (MQ=255)
                                                                                                                                                                                                                                          tttAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGCGCAATCATGGACGCGCTaa          <  2:279266/249‑1 (MQ=255)
                                                                                                                                                                                                                                                aaaaGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGCGCAATCATGGACGCGCTAAATATATcc  <  2:107073/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
GGCGATTATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGCGCAATCATGGACGCGCTAAATATATCC  >  NC_000913/2939052‑2939542

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