Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,151 2 bp→GA intergenic (+36/+77) yabI → / ← thiQ DedA family protein YabI/thiamine ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,1510TG93.1% 151.8 / 1.5 59intergenic (+36/+78)yabI/thiQDedA family protein YabI/thiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base G (29/25);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01
*NC_00091372,1520GA91.4% 156.4 / 3.4 59intergenic (+37/+77)yabI/thiQDedA family protein YabI/thiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (28/25);  total (32/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 = 7.50e-01

GTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTG‑TGCGGGTTGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTG  >  NC_000913/71922‑72393
                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                                  
gtggCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTAc                                                                                                                                                                                                                               <  1:604295/251‑1 (MQ=255)
               ggCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACtgtcc                                                                                                                                                                                                                 <  1:400523/250‑5 (MQ=255)
                 cGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGGGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACtgtccgt                                                                                                                                                                                                               >  1:543732/1‑244 (MQ=255)
                   gTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACtgtccttagg                                                                                                                                                                                                            <  1:514121/251‑10 (MQ=255)
                      aTGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGGGGCGCTGGTGGGGTTAATTCGCCACCCGTTGATGCCGGGGTATATCGAAATTTTGCGTAAAGTGGGTGGGGTTTAGGGGATTTTCTTT‑TGCGGGTTGCCTGATGCGACG‑TTGACCCCGTCTTATCAAGCCTACAAAACGCACTGTCCGTAGGTCg                                                                                                                                                                                                         >  1:700534/1‑251 (MQ=255)
                        ggCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGGTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCCGGCCTACAAAACGCACTGTCCGTAGGTCggg                                                                                                                                                                                                       >  1:183566/1‑250 (MQ=255)
                                   aaagcaaccgatcGTTTAAGTCATTATTTTTCGCGCGGACGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTTGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTc                                                                                                                                                                                             <  2:177059/238‑1 (MQ=255)
                                   aaagcaaccgatcGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTc                                                                                                                                                                                             <  1:306171/238‑1 (MQ=255)
                                     aGCGACTGACCGTTTTAGTCATTATTTGTCCCGCGGTCGTTTGTGGCGGCTGACGTCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCGCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTG‑TGCGGGTTGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAATGCGTTCAcg                                                                                                                                                                                          <  2:287141/251‑1 (MQ=255)
                                                  ttaaGTCATGATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACa                                                                                                                                                                              <  1:53349/247‑1 (MQ=255)
                                                    aagctattattTGTCGCGTGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCTTCGGCGGGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATAGCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCAAGCCGCATCTGACAGt                                                                                                                                                                            <  1:684587/245‑1 (MQ=255)
                                                       tCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGGGGTGTTAATTCGCCACCCGTTGGTGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGGCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGGTAAGGTGTTccagccgcctctgaccgagcat                                                                                                                                                                        >  1:538745/1‑230 (MQ=255)
                                                                gtcgcgcggcggTTTGTTGTGGGTGACGCCGTTGATTTCGGCCATCGGCGTGGTGGCGCGGGTGGTGTTAATTCGCCACCCGTTGATGCCGGGGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAAccc                                                                                                                                                                <  2:170254/239‑1 (MQ=255)
                                                                  cgcgcgGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGGTAAATCGCCACCCGTTTATGCCGGTGTATATCGATATTTTGCGGCAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGGGGCGCTGT‑CCGCGTCTTATACGGGCTCCAAAACGCACTGGCCGTAAGGAGGGTAAGGGGTTCACCCAGCATCCGACAGTGCAACCTAAccc                                                                                                                                                                >  2:485996/3‑248 (MQ=255)
                                                                     gcgGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCGTCTGACAGTGCATACTAACCCGTAATc                                                                                                                                                          <  1:536028/251‑1 (MQ=255)
                                                                       ggccGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGGTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATccc                                                                                                                                                        >  1:288289/4‑251 (MQ=255)
                                                                       ggccGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGGGGTAATTCGCCACCCGTTTATGCCGGTGAATATCGATATTTTGCGTAAAGGGGTTGGGGTTTAGGGGGTTGTTTTG‑TGCGGGGTGCCTGGTTCGACGCTGA‑CCGCGTCTTATCCAGCCTACAAAACCCCCTGTACGGAGGTGGGATAAGGTGTTAACGCCGCCTCCTAAAGTGCATCCTAACCCGTAAAccc                                                                                                                                                        >  2:496612/4‑251 (MQ=255)
                                                                              tgttgtGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCCCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAg                                                                                                                                                  <  1:308789/250‑1 (MQ=255)
                                                                                     gCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCGGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAg                                                                                                                                          <  1:297154/251‑1 (MQ=255)
                                                                                            ccGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGGGTTAATTCGCCACCCGGTGATGCCGGTGTATATCGGTATTTTGCGTAAAGTGGTTGGGGGTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTATTATCAGGCCTACAAAACGCACTGTCCGGAGGGTGGCAAAGGTGTTCAGGCCGCCTCTGACCGGGCATACAAACACGTAATCCCCCATAGTGCTGaaacccccgc                                                                                                                                   >  2:637950/1‑243 (MQ=255)
                                                                                                tGATTTCTGCCAGCGGCGTGGTGGCGCTGGTGGTGTTAATGCGCCGCCCGTGGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGCCGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCtt                                                                                                                                <  2:223891/250‑1 (MQ=255)
                                                                                                 gATTGCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTGGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTa                                                                                                                               <  1:708057/250‑1 (MQ=255)
                                                                                                   tttCTGCCTGCGGCTTGGTGGCGCTGTTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATTGTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTAcc                                                                                                                             <  2:188197/250‑1 (MQ=255)
                                                                                                        gCCATCGGCGTGGTTGCGCGGTTGGTTTTAAGTCGCAACCCGTTGATGCCGGTGTATATCGATATTTTGTGTAAAGGGGTGGGGGTTGAGGGGATTGTCTTG‑TGCGGGGTGCCTGATGCGCCGCTGA‑CCGCGTCTTATCAGGCCTACAAACCGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCt                                                                                                                          <  2:38484/248‑1 (MQ=255)
                                                                                                            tCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCAt                                                                                                                                                                                    <  2:636460/186‑1 (MQ=255)
                                                                                                            tCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCAt                                                                                                                                                                                    >  1:636460/1‑186 (MQ=255)
                                                                                                                       gtgttggtggGGTTAATTCGCCACCCGGTGATGCCGGGGTTTATCGATATTTTGGGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTctt                                                                                                        <  2:183566/247‑3 (MQ=255)
                                                                                                                          cTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATa                                                                                                                                                   <  1:494034/205‑1 (MQ=255)
                                                                                                                          cTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATa                                                                                                                                                   >  2:494034/1‑205 (MQ=255)
                                                                                                                           tggtggTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTAccc                                                                                                    >  1:566792/1‑251 (MQ=255)
                                                                                                                                     aaTTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGAGGTTGGGGTTTAGGGGATTGTCTTG‑TGAGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCCTACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCg                                                                                           >  2:82102/1‑250 (MQ=255)
                                                                                                                                             aCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCACTAGTGCTGAAGCACTCTCCTTACCGCTCACCAACTCATTCGTCTTACCCTGCCAGGCGCTGCTcccc                                                                                  >  2:419050/1‑250 (MQ=255)
                                                                                                                                                cGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCg                                                                               <  1:112324/251‑1 (MQ=255)
                                                                                                                                                        ccGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCgcccac                                                                       >  1:215203/1‑246 (MQ=255)
                                                                                                                                                        ccGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGaccac                                                                       >  2:146419/1‑248 (MQ=255)
                                                                                                                                                               atatCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGATCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGACCGTCGGCGACCACCACCg                                                                  >  1:420059/1‑249 (MQ=255)
                                                                                                                                                                    gATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGt                                                           >  1:246578/1‑251 (MQ=255)
                                                                                                                                                                      tatTTTGCTGAAAGTGTTTGGGGTTTAGGGTATTGTCTTG‑GCCGGGTTGCTGGATCCGCCGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTg                                                          <  2:12912/250‑1 (MQ=255)
                                                                                                                                                                       atTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTg                                                          >  1:12912/1‑249 (MQ=255)
                                                                                                                                                                        ttttGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGCCCACCACCGAGCCCGTGGCg                                                       >  2:649887/1‑251 (MQ=255)
                                                                                                                                                                          ttGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGAt                                                     >  1:364894/1‑251 (MQ=255)
                                                                                                                                                                                aaaGGGGTTGGGGGTTAGGGGATTGTCTGGGTG‑GGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGGTAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGGCTTTCCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCgcgc                                               <  2:441971/251‑1 (MQ=255)
                                                                                                                                                                                        tGGGGTTTAGGAGATAGTCTTG‑TGCGGGTTGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCtt                                       >  1:659108/1‑251 (MQ=255)
                                                                                                                                                                                         ggggTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCCTTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGATCGCGTGGCGCGCCGCGAGGCATCTTc                                      >  2:139350/1‑251 (MQ=255)
                                                                                                                                                                                          gggTTTAGGGGATTGTCTTG‑TGGGGTTTGCTTGATGCGCCG‑TGGACCGCGTCTTATCAGGCCTTAAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTcc                                     <  2:543732/251‑1 (MQ=255)
                                                                                                                                                                                          gggTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTcc                                     >  1:374851/1‑251 (MQ=255)
                                                                                                                                                                                              ttAGGTGATTGTCTTG‑TGCGGGTTGCCTGATGGGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGCCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAca                                  <  1:663953/250‑2 (MQ=255)
                                                                                                                                                                                                aggggattgTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGGCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGCCCACCACCGAGCGCGTGGCGGTCCGCGCGGCATCTTCCACACtg                               >  2:614622/9‑251 (MQ=255)
                                                                                                                                                                                                 ggggattgTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACtgt                              >  1:207465/8‑251 (MQ=255)
                                                                                                                                                                                                    gattgTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCCCACtgtg                             >  2:348239/5‑249 (MQ=255)
                                                                                                                                                                                                        gTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGCGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                       >  1:130820/1‑251 (MQ=255)
                                                                                                                                                                                                        gTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACCCTGTGTGacac                       >  1:348505/1‑251 (MQ=255)
                                                                                                                                                                                                            tg‑tgCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACAcc                      <  1:496612/248‑1 (MQ=255)
                                                                                                                                                                                                                   ggTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTg                                                                                                             >  2:117255/1‑155 (MQ=255)
                                                                                                                                                                                                                   ggTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTg                                                                                                             <  1:117255/155‑1 (MQ=255)
                                                                                                                                                                                                                    gTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGACGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGCCACCATCAATaacg            >  2:422700/1‑248 (MQ=255)
                                                                                                                                                                                                                      tGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCACTAACGTc          >  1:402160/1‑251 (MQ=255)
                                                                                                                                                                                                                            aTGCGACGCTTG‑CCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCAttttt    >  1:23102/1‑251 (MQ=255)
                                                                                                                                                                                                                              gCGACGCTCA‑CCGCGCCTTATCCGGCCTACAAAACGCACTGCCCGTAGGTCGGATAAGTTTTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACATCCCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATTCGCGCGGCATCTTCCACACTGTGTGACACCAGCAATAACGTCATTTTTtg  <  1:207806/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                                  
GTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTG‑TGCGGGTTGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTG  >  NC_000913/71922‑72393

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