Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,757,531 G→A 14.0% S99F (TCT→TTT)  yiaY ← L‑threonine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,757,5310GA14.0% 148.2 / 2.4 60S99F (TCT→TTT) yiaYL‑threonine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (37/14);  new base A (3/6);  total (40/20)
Fisher's exact test for biased strand distribution p-value = 4.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCGCTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCAT  >  NC_000913/3757298‑3757760
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                       
tgTCATAGACAGAAAGGAGGGGGTTGACATTTTAAAAAAAACGCGACATTTGTATGGGACGCGTTGCGTAAGTGATGATGCAAAAGGGGTTCTTTTCAGGGGCCGTCCCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAacc                                                                                                                                                                                                                                     <  2:1145685/237‑1 (MQ=255)
agTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCGCTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCATGCCCATACCTTTTGGCCAGCCGTGGGGGGAacc                                                                                                                                                                                                                                     >  2:1185605/1‑238 (MQ=255)
agTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCGCTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCCCCCTTGGCTGCCACCAGCGCAATACCTTTTGCCCAGTCGTGTGGAGAacc                                                                                                                                                                                                                                     >  2:486526/1‑238 (MQ=255)
agTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCGCTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTCCGGCTTTGACAGGCGGTCAACGCCTTCGTAATCCGGAATATCACCGCCATTGGCTGCGCCCATCCCAAAACCTTTTGCACAGCAGTGTGGAAAaca                                                                                                                                                                                                                                     >  1:138554/1‑237 (MQ=255)
  tCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCGCTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCCCCACCGCAATACCTTTTGCGCAGTCGTGTGGAGAaccac                                                                                                                                                                                                                                   >  2:600087/1‑238 (MQ=255)
  tCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCGCTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGAGGTTTTGCAGAGCGGTCAACGCCCTCGTAATCGCGAATATCGCCGCCCTTTCGTGCCCCCAGCGCAATACCTTTTGCGCAAGCGCGTGGAGGaccac                                                                                                                                                                                                                                   >  2:932679/1‑238 (MQ=255)
         aaaaaagcggggggtgaaCATTTTTACCAACAATCGCCTTTTGGATAGGCCGCGCTTCGGCAGTGATGATGCAAAAGGGGGTCATTCCAGGGCCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAg                                                                                                                                                                                                                          <  2:455745/222‑1 (MQ=255)
          aGAAAGCAGCGAAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCGCTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGACGTGTGGGGAACCACCGCCCAAgg                                                                                                                                                                                                                         >  2:906180/1‑239 (MQ=255)
              agcagcGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCGCTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGGGAACCACCGCCTAAGgaga                                                                                                                                                                                                                      >  2:395147/1‑238 (MQ=255)
                   cGGAGTGACATGTTTATCAACAATCGACATTTTGATATGACGTGCTTCGTCAGTGCTGATGCAGAAACGGGTCATTTCAGAGGCCGTACCAGCCGTGGTATTGCTGGCGATAATCGTCAGATTCGGTTTTGTCGTGCGGTGAACGGCTTCCAAAGCCGGGATGTCGCCGCCATCTGCAGCCCCAAAGTGCCTACCCTTTGCGCCTTCGTGTGgagacataagcctaaaagggaaacgga                                                                                                                                                                                                                 >  1:875051/1‑216 (MQ=255)
                    ggAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCGCTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCCCCATTGGCTGCCACCAGCCCAATACCTTTTTCGCAGTCGTGTGGGAGACCACCGCCTAAGGAGATCACGCt                                                                                                                                                                                                               >  2:582421/1‑239 (MQ=255)
                      aGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCGCTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCCATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCtt                                                                                                                                                                                                              >  2:1106262/1‑237 (MQ=255)
                       gTGACTTTTTTAACAACAACCGCTGTTGTGATAGGCCGCGCTTCGTCAGTGAGGGTGGAAAAGCGGTTTTTTACAGGGGCCTACCCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTAt                                                                                                                                                                                                             <  2:516877/238‑1 (MQ=255)
                               tttacaaccatccgcttttttGATATGCGCGCCTTGGTCAGTATGGTTGCAGACCCGGGTTATTCAAGGGCCCTTCCCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAattat                                                                                                                                                                                                   <  2:1242715/222‑1 (MQ=255)
                                  aTCAACAATCGCCATTTTGATATGACGCGCTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGGGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAACGCAATACCTTTTGCGCAGTCGTGTGGCAAACCACCACCTACggcaatacggcattccgaatttttt                                                                                                                                                                                                  >  2:469717/1‑215 (MQ=255)
                                          tccgctttttGATAGGCGGCGCTTCGTCAGGTAGTATGAAGAAACGGTTAATTTCGGGGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGt                                                                                                                                                                                        <  2:1277402/233‑1 (MQ=255)
                                          tCGCCATTTTGATATGACGCGCTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCCCCAGCGCAATACCTTTTGCGCAGTCGTGGGGAGAACCACCGCCTAAGGAGATCACGCTATCCCAATTATTCTCTTTAAGt                                                                                                                                                                                        >  2:1170260/1‑239 (MQ=255)
                                                        tGACGCGCTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCCCGCTATCGCACTTATTCTCTTTAAGTAATTTCAAACCTgc                                                                                                                                                                          >  2:1275056/1‑239 (MQ=255)
                                                           cgcgCTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTgcggc                                                                                                                                                                       >  2:1200775/1‑239 (MQ=255)
                                                                cTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATTCCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTGTTGCGCAGTCGTGTGGAGAACCCCCGCCAAAGGGGATCACACCATCGCACATTTTCCTTTTTAAGAATTTAAAACCTagggcggag                                                                                                                                                                   >  1:1080883/1‑236 (MQ=255)
                                                                cTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTTTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCCATTATTCTCTTTAAGTAATTTCAAACCTGCGGGGACg                                                                                                                                                                   >  2:18479/1‑239 (MQ=255)
                                                                   cGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCCCCGCCTAAGGAGATCACCACATCCCAATTATTCTCTTTAAGGAATTTCAAACCTGCGGGGACGttt                                                                                                                                                                >  2:1114431/1‑239 (MQ=255)
                                                                     gccatgtagggtgaaaaaCGGGGCATTTGAGGGGCCCTGCCCGGCCGTGGTATGGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGtttt                                                                                                                                                               <  2:122763/224‑1 (MQ=255)
                                                                           atgatgCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCgtggt                                                                                                                                                        >  2:702732/1‑239 (MQ=255)
                                                                              atgCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCCCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTggg                                                                                                                                                     >  2:784136/1‑239 (MQ=255)
                                                                                    aaggggtgtttttaggtggcctgaCCCGCGGGGGTTATGGTGGCGATCACCGGCAGCTGCGTTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAgg                                                                                                                                               <  2:957044/216‑1 (MQ=255)
                                                                                      aCGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGtt                                                                                                                                             >  2:858077/1‑239 (MQ=255)
                                                                                           tatttcgagggccgcgaCCCGCGTTGGTATGGATGGCGATCATCGGACGCTGCGGTTTTGCAGGGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTg                                                                                                                                        <  2:573819/223‑1 (MQ=255)
                                                                                           cttttcgtgggggcttccccGCCGTGGTTTGTATGGCGATCGCCGGGATCTGCGGTTTTGAGAGGCGGTAAACGCCGTCTTAATCCGGAATTCCGCCGCCATTGGGTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTg                                                                                                                                        <  1:18479/222‑1 (MQ=255)
                                                                                            ctttccgggggcgttACCCGCCGTGGTATTGATGGCGATCATCGGAAGCTGCGGTTTTGCAGGGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTg                                                                                                                                        <  2:1053227/224‑1 (MQ=255)
                                                                                                tCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAg                                                                                                                                                                                                                          >  2:809598/1‑153 (MQ=255)
                                                                                                            acccggtggtggtttGATGGGGATCATCGGCAGGTGCGTTTTGGCAGAGCGGTCAACGCTTTCGTAACCGCGAATAGCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAATCCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGc                                                                                                                        <  2:533244/226‑1 (MQ=255)
                                                                                                                    tGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAata                                                                                                                >  2:1171860/1‑239 (MQ=255)
                                                                                                                              ggCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGttct                                                                                                       >  2:959910/1‑238 (MQ=255)
                                                                                                                                             cTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTTCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTTAAAAATTTTCGTTCTTCCAGTGCTTTTTGc                                                                                        >  2:1298481/1‑238 (MQ=255)
                                                                                                                                               gCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGGGGAGAACCCCCGCCTAAGGAGATCACGCTATCGCAATGATTCTCTTTAAGTAATTTCAAACCTGCGGCCACGTTTTCCGTGTGGGGGTTAGGTTGGGGGCCACCATAAAAAACCACAAAAAAATTTGGTTTTTcccagcgttttttaact                                                                                     >  1:721156/1‑224 (MQ=255)
                                                                                                                                                gggTTTTGCAGGGGGGGCAACGTCTTGATAACGGGAAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGAGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGcaca                                                                                     <  2:12342/237‑1 (MQ=255)
                                                                                                                                                               cGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATAcc                                                                       >  2:1047308/1‑239 (MQ=255)
                                                                                                                                                               cGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGCAATTTCAAACCTGCGGCGACGTTTTCCGTGGCGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCCCATAGCCCGCCATAcc                                                                       >  2:328529/1‑239 (MQ=255)
                                                                                                                                                                         cTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCCGTGCCTTTTGCACATCGCCCGCCATACCCAATTTct                                                               >  2:90947/1‑236 (MQ=255)
                                                                                                                                                                         cTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGGGCCATCCTAAATAACGCTAAAAATATTGCTTTCTTCCAGGGCTTTTTGCACATCTCCCGCCATAACTAATTTCGtt                                                             >  2:851607/1‑239 (MQ=255)
                                                                                                                                                                           ctggaacccacaatttcccctataTTGGCTCCAACCAGCGCAACTCCTTTGGCGCAGTGGTGTGGAGACCCCCCGCCTAAGGAGATCACGAGTTCGAAATTATTCTCTTTAAGTAATTTCAAACCGGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGTTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTaa                                                           <  2:34479/216‑1 (MQ=255)
                                                                                                                                                                                aTCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACatat                                                      >  2:241500/1‑239 (MQ=255)
                                                                                                                                                                                 tCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCCTTTGCGCAGTCGTGTGAAGAACCACCGCCTAAGGAGATCACGCTATTGCAATTATTCTCTTTAAGTAATTTCAAACCTTCGGCGACGTTTTCCTTGGTGGGTCTCGGTAGGGTGCCATCATAAATAACGTTAAAAAAATTGCGTTCTTCTAGTGCTTTTTGTACTTAGCCGCACATACCTAATTTTGTaaaattata                                                     >  1:1263247/1‑233 (MQ=255)
                                                                                                                                                                                              gcccATTGGTTGCCCCCAGCGCAATCCCTTTTGCGCATTCGTGTGGAAACCCCCCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCGGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGTTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAAtt                                        <  2:721156/237‑1 (MQ=255)
                                                                                                                                                                                                   gtGGCGGCCACCCGCGCAATCCCTTGTGCGCTGTTGTGTGGAAACCCCCCGCCAAAGGAGATCACGTTATCGAAATTATTCCTTTTAAGTATTGTCAACCCGGCGGCGACGTTTTCCGGGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGt                                   <  2:672703/238‑1 (MQ=255)
                                                                                                                                                                                                           caccaGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTaa                           >  2:592419/1‑239 (MQ=255)
                                                                                                                                                                                                                aGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAACATATTGCGTTCTTCCAGGGCTTTTTTCACATCGCCCGCCAAAACTAATTTCGTTAACCAATGTTCAGTGACAATAAAGGTACGGGAAAAAcc                       >  1:535082/1‑239 (MQ=255)
                                                                                                                                                                                                                      aTACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAAt                  >  2:652667/1‑238 (MQ=255)
                                                                                                                                                                                                                        cccTTTTGCGAAGTCGTGTGGAAACCCCCCCCCTAAGGAGTTCACCCTATCGCAATTATTCTTTTTAAGTAATTCCAAACCTGCGGCGACGTTTCCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATTTTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTg               <  2:197255/238‑1 (MQ=255)
                                                                                                                                                                                                                         cttttttccgtgTCGTGGGGGAAACCCCCCCCTTAAGGGTATCCCGTTATCGCATTTATTCTCTTTAAGTAATGACAACCCGGGGGCACCGTTTTCGGTGGGGGGTTTAGGTTGGGGTCATTCATAAATAACGATAAAATTTTTGCTTTCTTCCAGTGCTTTTTGCCCATCTCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTg               <  1:959910/228‑1 (MQ=255)
                                                                                                                                                                                                                            tttGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAACAATATTGCGTTCTTCCAGTGCTTTTTGCCCATCGCCCGCCATACCTAATTTCGTTACCATATTGTCAGGGACAATTAAGGTACGGGTaaaacaaaaaaatcccaa           >  1:197255/1‑224 (MQ=255)
                                                                                                                                                                                                                            tttGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCcat           >  2:206684/1‑239 (MQ=255)
                                                                                                                                                                                                                                 gcAGTCGGTGGGAGAACCCCCTCCAAAGGAGATCACGCTATCGCATTTATTCTCTTTAAGTAATTTCAACCCTGCGGCGACGTTTTCCGTGGTGGGTTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCata       <  2:753528/238‑1 (MQ=255)
                                                                                                                                                                                                                                   aGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAg                                                >  2:409985/1‑195 (MQ=255)
                                                                                                                                                                                                                                    ggcggggggAAAACCACCGCTGAGGGAGATCCCGTTATCGCAATTATTCTCTTAAAGTAATTCCAACCCGGCGGCGACGTTTCCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCa   <  2:249176/232‑1 (MQ=255)
                                                                                                                                                                                                                                    ggcGTGGGGAAACCCCCCGACTAAGGAGATCAGGTTTTCGAAATTTTTCTCTTAAAGTAATTTAAAACCTGGGGCGACGTTTTCCGTGGTGGGGTTAGTTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCa   <  2:138554/237‑1 (MQ=255)
                                                                                                                                                                                                                                    gTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTCACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCAAAATTTGCcaatacattaa   >  1:254644/1‑230 (MQ=255)
                                                                                                                                                                                                                                    gTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGAGCTTTTTGCACATCGCCCTCCAACCCTAATTTCGTTAACATATTTTCATCGGCAATTAAGGTACGGTGAAAACCATAAACTGAcataacataaa   >  1:677542/1‑231 (MQ=255)
                                                                                                                                                                                                                                     tCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATTTTGTCAGTGACAATTCAGGAACCGGTAAACCCAAAATATCCCATCATAATCAt  >  1:310368/1‑239 (MQ=255)
                                                                                                                                                                                                                                     tCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATATGCCATCATATTCAt  >  2:40324/1‑239 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                       
AGTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCGCTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCCGCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTGCAGAGCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCAT  >  NC_000913/3757298‑3757760

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