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

GTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTG  >  NC_000913/71922‑72393
                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                                  
gtggCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTAc                                                                                                                                                                                                                               <  1:604295/251‑1 (MQ=255)
               ggCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACtgtcc                                                                                                                                                                                                                 <  1:400523/250‑5 (MQ=255)
                 cGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGGGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACtgtccgt                                                                                                                                                                                                               >  1:543732/1‑244 (MQ=255)
                   gTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACtgtccttagg                                                                                                                                                                                                            <  1:514121/251‑10 (MQ=255)
                      aTGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGGGGCGCTGGTGGGGTTAATTCGCCACCCGTTGATGCCGGGGTATATCGAAATTTTGCGTAAAGTGGGTGGGGTTTAGGGGATTTTCTTTTGCGGGTTGCCTGATGCGACGTTGACCCCGTCTTATCAAGCCTACAAAACGCACTGTCCGTAGGTCg                                                                                                                                                                                                         >  1:700534/1‑251 (MQ=255)
                        ggCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGGTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCCGGCCTACAAAACGCACTGTCCGTAGGTCggg                                                                                                                                                                                                       >  1:183566/1‑250 (MQ=255)
                                   aaagcaaccgatcGTTTAAGTCATTATTTTTCGCGCGGACGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTTGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTc                                                                                                                                                                                             <  2:177059/238‑1 (MQ=255)
                                   aaagcaaccgatcGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTc                                                                                                                                                                                             <  1:306171/238‑1 (MQ=255)
                                     aGCGACTGACCGTTTTAGTCATTATTTGTCCCGCGGTCGTTTGTGGCGGCTGACGTCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCGCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAATGCGTTCAcg                                                                                                                                                                                          <  2:287141/251‑1 (MQ=255)
                                                  ttaaGTCATGATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACa                                                                                                                                                                              <  1:53349/247‑1 (MQ=255)
                                                    aagctattattTGTCGCGTGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCTTCGGCGGGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATAGCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCAAGCCGCATCTGACAGt                                                                                                                                                                            <  1:684587/245‑1 (MQ=255)
                                                       tCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGGGGTGTTAATTCGCCACCCGTTGGTGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGGCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGGTAAGGTGTTccagccgcctctgaccgagcat                                                                                                                                                                        >  1:538745/1‑230 (MQ=255)
                                                                gtcgcgcggcggTTTGTTGTGGGTGACGCCGTTGATTTCGGCCATCGGCGTGGTGGCGCGGGTGGTGTTAATTCGCCACCCGTTGATGCCGGGGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAAccc                                                                                                                                                                <  2:170254/239‑1 (MQ=255)
                                                                  cgcgcgGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGGTAAATCGCCACCCGTTTATGCCGGTGTATATCGATATTTTGCGGCAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGGGGCGCTGTCCGCGTCTTATACGGGCTCCAAAACGCACTGGCCGTAAGGAGGGTAAGGGGTTCACCCAGCATCCGACAGTGCAACCTAAccc                                                                                                                                                                >  2:485996/3‑248 (MQ=255)
                                                                     gcgGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCGTCTGACAGTGCATACTAACCCGTAATc                                                                                                                                                          <  1:536028/251‑1 (MQ=255)
                                                                       ggccGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGGTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATccc                                                                                                                                                        >  1:288289/4‑251 (MQ=255)
                                                                       ggccGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGGGGTAATTCGCCACCCGTTTATGCCGGTGAATATCGATATTTTGCGTAAAGGGGTTGGGGTTTAGGGGGTTGTTTTGTGCGGGGTGCCTGGTTCGACGCTGACCGCGTCTTATCCAGCCTACAAAACCCCCTGTACGGAGGTGGGATAAGGTGTTAACGCCGCCTCCTAAAGTGCATCCTAACCCGTAAAccc                                                                                                                                                        >  2:496612/4‑251 (MQ=255)
                                                                              tgttgtGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCCCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAg                                                                                                                                                  <  1:308789/250‑1 (MQ=255)
                                                                                     gCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCGGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAg                                                                                                                                          <  1:297154/251‑1 (MQ=255)
                                                                                            ccGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGGGTTAATTCGCCACCCGGTGATGCCGGTGTATATCGGTATTTTGCGTAAAGTGGTTGGGGGTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTATTATCAGGCCTACAAAACGCACTGTCCGGAGGGTGGCAAAGGTGTTCAGGCCGCCTCTGACCGGGCATACAAACACGTAATCCCCCATAGTGCTGaaacccccgc                                                                                                                                   >  2:637950/1‑243 (MQ=255)
                                                                                                tGATTTCTGCCAGCGGCGTGGTGGCGCTGGTGGTGTTAATGCGCCGCCCGTGGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGCCGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCtt                                                                                                                                <  2:223891/250‑1 (MQ=255)
                                                                                                 gATTGCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTGGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTa                                                                                                                               <  1:708057/250‑1 (MQ=255)
                                                                                                   tttCTGCCTGCGGCTTGGTGGCGCTGTTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATTGTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTAcc                                                                                                                             <  2:188197/250‑1 (MQ=255)
                                                                                                        gCCATCGGCGTGGTTGCGCGGTTGGTTTTAAGTCGCAACCCGTTGATGCCGGTGTATATCGATATTTTGTGTAAAGGGGTGGGGGTTGAGGGGATTGTCTTGTGCGGGGTGCCTGATGCGCCGCTGACCGCGTCTTATCAGGCCTACAAACCGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCt                                                                                                                          <  2:38484/248‑1 (MQ=255)
                                                                                                            tCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCAt                                                                                                                                                                                    <  2:636460/186‑1 (MQ=255)
                                                                                                            tCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCAt                                                                                                                                                                                    >  1:636460/1‑186 (MQ=255)
                                                                                                                       gtgttggtggGGTTAATTCGCCACCCGGTGATGCCGGGGTTTATCGATATTTTGGGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTctt                                                                                                        <  2:183566/247‑3 (MQ=255)
                                                                                                                          cTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATa                                                                                                                                                   <  1:494034/205‑1 (MQ=255)
                                                                                                                          cTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATa                                                                                                                                                   >  2:494034/1‑205 (MQ=255)
                                                                                                                           tggtggTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTAccc                                                                                                    >  1:566792/1‑251 (MQ=255)
                                                                                                                                     aaTTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGAGGTTGGGGTTTAGGGGATTGTCTTGTGAGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCCTACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCg                                                                                           >  2:82102/1‑250 (MQ=255)
                                                                                                                                             aCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCACTAGTGCTGAAGCACTCTCCTTACCGCTCACCAACTCATTCGTCTTACCCTGCCAGGCGCTGCTcccc                                                                                  >  2:419050/1‑250 (MQ=255)
                                                                                                                                                cGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCg                                                                               <  1:112324/251‑1 (MQ=255)
                                                                                                                                                        ccGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCgcccac                                                                       >  1:215203/1‑246 (MQ=255)
                                                                                                                                                        ccGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGaccac                                                                       >  2:146419/1‑248 (MQ=255)
                                                                                                                                                               atatCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGATCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGACCGTCGGCGACCACCACCg                                                                  >  1:420059/1‑249 (MQ=255)
                                                                                                                                                                    gATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGt                                                           >  1:246578/1‑251 (MQ=255)
                                                                                                                                                                      tatTTTGCTGAAAGTGTTTGGGGTTTAGGGTATTGTCTTGGCCGGGTTGCTGGATCCGCCGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTg                                                          <  2:12912/250‑1 (MQ=255)
                                                                                                                                                                       atTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTg                                                          >  1:12912/1‑249 (MQ=255)
                                                                                                                                                                        ttttGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGCCCACCACCGAGCCCGTGGCg                                                       >  2:649887/1‑251 (MQ=255)
                                                                                                                                                                          ttGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGAt                                                     >  1:364894/1‑251 (MQ=255)
                                                                                                                                                                                aaaGGGGTTGGGGGTTAGGGGATTGTCTGGGTGGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGGTAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGGCTTTCCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCgcgc                                               <  2:441971/251‑1 (MQ=255)
                                                                                                                                                                                        tGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCtt                                       >  1:659108/1‑251 (MQ=255)
                                                                                                                                                                                         ggggTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCCTTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGATCGCGTGGCGCGCCGCGAGGCATCTTc                                      >  2:139350/1‑251 (MQ=255)
                                                                                                                                                                                          gggTTTAGGGGATTGTCTTGTGGGGTTTGCTTGATGCGCCGTGGACCGCGTCTTATCAGGCCTTAAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTcc                                     <  2:543732/251‑1 (MQ=255)
                                                                                                                                                                                          gggTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTcc                                     >  1:374851/1‑251 (MQ=255)
                                                                                                                                                                                              ttAGGTGATTGTCTTGTGCGGGTTGCCTGATGGGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGCCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAca                                  <  1:663953/250‑2 (MQ=255)
                                                                                                                                                                                                aggggattgTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGGCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGCCCACCACCGAGCGCGTGGCGGTCCGCGCGGCATCTTCCACACtg                               >  2:614622/9‑251 (MQ=255)
                                                                                                                                                                                                 ggggattgTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACtgt                              >  1:207465/8‑251 (MQ=255)
                                                                                                                                                                                                    gattgTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCCCACtgtg                             >  2:348239/5‑249 (MQ=255)
                                                                                                                                                                                                        gTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGCGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                       >  1:130820/1‑251 (MQ=255)
                                                                                                                                                                                                        gTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACCCTGTGTGacac                       >  1:348505/1‑251 (MQ=255)
                                                                                                                                                                                                            tgtgCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACAcc                      <  1:496612/248‑1 (MQ=255)
                                                                                                                                                                                                                   ggTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTg                                                                                                             >  2:117255/1‑155 (MQ=255)
                                                                                                                                                                                                                   ggTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTg                                                                                                             <  1:117255/155‑1 (MQ=255)
                                                                                                                                                                                                                    gTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGACGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGCCACCATCAATaacg            >  2:422700/1‑248 (MQ=255)
                                                                                                                                                                                                                      tGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCACTAACGTc          >  1:402160/1‑251 (MQ=255)
                                                                                                                                                                                                                            aTGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCAttttt    >  1:23102/1‑251 (MQ=255)
                                                                                                                                                                                                                              gCGACGCTCACCGCGCCTTATCCGGCCTACAAAACGCACTGCCCGTAGGTCGGATAAGTTTTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACATCCCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATTCGCGCGGCATCTTCCACACTGTGTGACACCAGCAATAACGTCATTTTTtg  <  1:207806/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                                  
GTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTG  >  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: