Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,202 C→T intergenic (+87/+27) 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,2020CT96.2% 167.6 / NA 53intergenic (+87/+27)yabI/thiQDedA family protein YabI/thiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/16);  minor base G (2/0);  total (37/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGG  >  NC_000913/71960‑72444
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
gcaaccgatcGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCAcgc                                                                                                                                                                                                                                            <  1:66717/242‑1 (MQ=255)
 caaccgttcGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCAcgcc                                                                                                                                                                                                                                           <  2:205196/243‑1 (MQ=255)
   accgatcGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGGGGCGCTGGTGGGGGTAATTCGCCACCCGTTGATGCCGGTGGATATCGATATTTTGCGGAAAGGGGGTGGGGGTTAGGGGCTTGTCTTGGGCGGGGTGGCTGATGGGACGCTGACCGCGGCTTTTCAGGCCTACAAAACGCACTGTCCGGAGGTCGGACAAGGAGTTCAcgccgc                                                                                                                                                                                                                                         >  2:32567/7‑251 (MQ=255)
               aGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGGGGCGCTGGTGGGGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGGTGGGGTTTAGGGGATTGTCTTTTGCGGGTTGCCTGATGCGGCGCTGACCGCGGCTTATTAGGCCTACAAAACGCACTGTCCGTAGGTCGGGTAAGGGGTTCACGCCGCATCTGACAgggc                                                                                                                                                                                                                             >  2:65469/1‑248 (MQ=255)
                  cATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCCGGTGGAGTTAATTCGCCCCCCGTTGAAGCCGGTGTATATCGATATTTTGCGTAAAGGGGGTGGGGATTATGGGATTGTCTTGTGCGGGTTGCCTGATGCGAAGCTGACCGCGTCTTATCAGGCCTCCAAAACGAACTGTCCGTAGGTCGGATAAGGTGTTCCCGCCGCATCTGACAGGGCAt                                                                                                                                                                                                                           >  1:251259/1‑250 (MQ=255)
                  cATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGGGGCGCTGGGGGGGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGATGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGGTAAGGTGTTCACGCCGCATCTGACAGTGCAt                                                                                                                                                                                                                           >  1:182143/1‑250 (MQ=255)
                                cggccGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCGGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATcc                                                                                                                                                                                                            <  2:136175/247‑1 (MQ=255)
                                   ccGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGAGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGGGTATATCGATATTTTGCGTCAAGGGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGGCTGATGCGAAGCTGAACGCGGCTTATCAGGCCTACAAAACGCACTGTCCGTAGGGCGGATAAGGTGTTTAAGCCGCATCTGACAGTGCATACTAACCCGTAATcccc                                                                                                                                                                                                          >  1:191055/2‑250 (MQ=255)
                                    cGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCaa                                                                                                                                                                                                        <  1:128202/251‑1 (MQ=255)
                                                    cGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTc                                                                                                                                                                                        >  1:27487/1‑251 (MQ=255)
                                                             tttCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGGTTAGGGGATTGTCTTGTGCGGGGTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGGCGGATAAGGGGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCCATAGTGCTGAAGCACTCGCCTTAccc                                                                                                                                                                               >  2:313110/1‑250 (MQ=255)
                                                                  gCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGGGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCt                                                                                                                                                                             <  2:91042/248‑1 (MQ=255)
                                                                        ggCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCAcgccg                                                                                                                                                                                                                                          <  2:309889/181‑1 (MQ=255)
                                                                        ggCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCAcgccg                                                                                                                                                                                                                                          >  1:309889/1‑181 (MQ=255)
                                                                                       gtggtgTTAATTCGCCGCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTg                                                                                                                                                     <  1:59790/251‑1 (MQ=255)
                                                                                                ggtCGCCCCCCGATGATGCTGGTGTATAGCGATATTTGGCGTAAAGGTGTTGGGTTTGAGTGGATTGTCTTGTGCGGGTGTCCTGAGGCGACGCTGCCCGCGTGTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGTTGTTCGCGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCGCAACAACGCATGGGTCTTACCCTGCCAGGCGAt                                                                                                                                            <  2:345378/249‑1 (MQ=255)
                                                                                                       aCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCg                                                                                                                                     <  1:63503/251‑1 (MQ=255)
                                                                                                                      tgtATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGaccaccacc                                                                                                                      >  1:162756/1‑244 (MQ=255)
                                                                                                                           atCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGCTTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAgcc                                                                                                                 >  2:163214/1‑250 (MQ=255)
                                                                                                                              gATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCg                                                                                                                                                                                                                  <  2:91066/151‑1 (MQ=255)
                                                                                                                              gATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCg                                                                                                                                                                                                                  >  1:91066/1‑151 (MQ=255)
                                                                                                                                     tGCGTAAAGTGGATGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATc                                                                                                       >  2:328277/1‑251 (MQ=255)
                                                                                                                                            aGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCgcgcgg                                                                                                >  1:370259/1‑249 (MQ=255)
                                                                                                                                              tGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCgg                                                                                                                                 >  1:175346/1‑216 (MQ=255)
                                                                                                                                              tGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCgg                                                                                                                                 <  2:175346/216‑1 (MQ=255)
                                                                                                                                              gggtgggggTTTAGGGGATTGTCGTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGCCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCGGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCa                                                                                              <  2:370910/246‑1 (MQ=255)
                                                                                                                                                       tttAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAc                                                                                      >  1:169701/1‑250 (MQ=255)
                                                                                                                                                             ggattgTCTTGTGCGGGTTGCCTGATGCGACTCTGACCGCGTCTTATCAGGACTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCATTACCGCTCAACAACTCATTGGTCTTACCCTGGCAGGCGATGCGCCCGTCGGCGACCCCCACCGAGCGCGTGGCGCTCCGCACGGCATCTTCAACACtgtgt                                                                               >  1:174063/6‑250 (MQ=255)
                                                                                                                                                                 tgTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGac                                                                            <  1:88981/249‑1 (MQ=255)
                                                                                                                                                                  gTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                          >  1:339558/1‑251 (MQ=255)
                                                                                                                                                                  gTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                          >  1:125592/1‑251 (MQ=255)
                                                                                                                                                                  gTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGaca                                                                           <  2:136190/250‑1 (MQ=255)
                                                                                                                                                                   tCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACAcc                                                                         >  1:284357/1‑251 (MQ=255)
                                                                                                                                                                           gggTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATa                                                                  >  1:164032/1‑250 (MQ=255)
                                                                                                                                                                                          acgctgaccGCGTCTTATCAGGCCTACAAAACGCACTGTCCGAAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTtgctg                                                  >  1:328718/8‑251 (MQ=255)
                                                                                                                                                                                                 ccGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGc                                             >  1:127592/1‑249 (MQ=255)
                                                                                                                                                                                                   gcggcTTATCAGGCCACCAAAACGCCCTGTCCGTAGGGCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCt                                         <  2:199131/247‑1 (MQ=255)
                                                                                                                                                                                                    cgTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTc                                        >  1:138545/1‑251 (MQ=255)
                                                                                                                                                                                                    cgTCTTATCAGGCCTACAAAACGAACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTc                                        >  2:324013/1‑251 (MQ=255)
                                                                                                                                                                                                    cgTCTTATCAGGCCTACAAAAAGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCAGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGACTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCATCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTc                                        >  1:28913/1‑251 (MQ=255)
                                                                                                                                                                                                                     aaaaCGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTc                         >  2:115187/1‑249 (MQ=255)
                                                                                                                                                                                                                     aaaaCGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTc                         >  2:137726/1‑249 (MQ=255)
                                                                                                                                                                                                                     aaaaCGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTc                         >  2:66521/1‑249 (MQ=255)
                                                                                                                                                                                                                          gcactgtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCAtct                  >  2:176318/8‑251 (MQ=255)
                                                                                                                                                                                                                          gcactgtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCAtct                  >  2:244891/8‑251 (MQ=255)
                                                                                                                                                                                                                             ctgtccGTAGGTCGGCTAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGTAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCCGCTCGTGCTCACCAGCGTCAGCATCTCCt               >  2:40563/5‑251 (MQ=255)
                                                                                                                                                                                                                              tgtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATCACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTg              >  2:71641/4‑251 (MQ=255)
                                                                                                                                                                                                                              tgtccGTAGGTAGGATAAGGTGTTCCCGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCCTCTTCCCCACTGTGTGCCACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTccgcatcttctg              >  2:23308/4‑240 (MQ=255)
                                                                                                                                                                                                                                 ccGTAGGACGGATAAGGTGTTCAAGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTCCCCTGCCAGGCGATGCGCCAGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCCTCAATCACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCCTCTCCTGACg           >  2:366955/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGCGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAAc       >  2:194193/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ggTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACg      >  2:368164/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cGGAGAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGGGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCg   <  2:328718/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ggATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCgg  >  2:72157/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
GCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGG  >  NC_000913/71960‑72444

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