Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,183 2 bp→GT intergenic (+68/+45) 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,1830AG92.5% 187.8 / 3.6 67intergenic (+68/+46)yabI/thiQDedA family protein YabI/thiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base G (34/28);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01
*NC_00091372,1840CT90.9% 179.3 / 6.1 67intergenic (+69/+45)yabI/thiQDedA family protein YabI/thiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base T (33/27);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

ATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTG‑TGCGGGTTGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTG  >  NC_000913/71944‑72432
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                        
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                                                                                                                                                    <  1:117255/155‑1 (MQ=255)
                                                                                                                                                                                             ggTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTg                                                                                                                                                    >  2:117255/1‑155 (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)
                                                                                                                                                                                                             cgctgaccGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTtgctgc                                     <  1:82102/245‑1 (MQ=255)
                                                                                                                                                                                                                    cgcgTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGc                              <  1:391518/251‑1 (MQ=255)
                                                                                                                                                                                                                           tATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGc                        <  1:649887/250‑1 (MQ=255)
                                                                                                                                                                                                                            aTCAGGCCTACAAAACGCACTACCCGTAGGTCGGCTAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCACCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCCTTTTTTGCTGCTGGCAGCTCGGGCTc                      >  2:325396/1‑251 (MQ=255)
                                                                                                                                                                                                                               aGGCCTACAACACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTACCCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGACGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGCGTGACACCATAAATAACGTCATTTTTTGCTGCTGGCACCTCGGGCTAAcc                   >  2:533882/1‑251 (MQ=255)
                                                                                                                                                                                                                               aGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCAcc                   >  1:253871/1‑251 (MQ=255)
                                                                                                                                                                                                                               aGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGGCATTTTTTGCTGCTGGCAGCTCGTGCTCAcc                   >  2:172565/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCg               >  1:552869/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aCAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTc             <  1:419050/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCa            <  1:538629/251‑1 (MQ=255)
                                                                                                                                                                                                                                       caaaCGCACTGTCCGTAGGCCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGCCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTc             >  2:87633/2‑249 (MQ=255)
                                                                                                                                                                                                                                                tACCCGTAGGTCAGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTg  >  1:55711/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                        
ATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTG‑TGCGGGTTGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTG  >  NC_000913/71944‑72432

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