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

AGTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCG‑CTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCAT  >  NC_000913/3757298‑3757760
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                       
tgTCATAGACAGAAAGGAGGGGGTTGACATTTTAAAAAAAACGCGACATTTGTATGGGACGCG‑TTGCGTAAGTGATGATGCAAAAGGGGTTCTTTTCAGGGGCCGTCCCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAacc                                                                                                                                                                                                                                     <  2:1145685/237‑1 (MQ=255)
agTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCG‑CTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CATGCCCATACCTTTTGGCCAGCCGTGGGGGGAacc                                                                                                                                                                                                                                     >  2:1185605/1‑238 (MQ=255)
agTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCG‑CTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCCCCCTTGGCTGCCAC‑CAGCGCAATACCTTTTGCCCAGTCGTGTGGAGAacc                                                                                                                                                                                                                                     >  2:486526/1‑238 (MQ=255)
agTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCG‑CTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTCCGGCTTTGACAG‑‑GCGGTCAACGCCTTCGTAATCCGGAATATCACCGCCATTGGCTGCGCC‑CATCCCAAAACCTTTTGCACAGCAGTGTGGAAAaca                                                                                                                                                                                                                                     >  1:138554/1‑237 (MQ=255)
  tCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCG‑CTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCCC‑CACCGCAATACCTTTTGCGCAGTCGTGTGGAGAaccac                                                                                                                                                                                                                                   >  2:600087/1‑238 (MQ=255)
  tCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCG‑CTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGAGGTTTTG‑CAGA‑GCGGTCAACGCCCTCGTAATCGCGAATATCGCCGCCCTTTCGTGCCCC‑CAGCGCAATACCTTTTGCGCAAGCGCGTGGAGGaccac                                                                                                                                                                                                                                   >  2:932679/1‑238 (MQ=255)
         aaaaaagcggggggtgaaCATTTTTACCAACAATCGCCTTTTGGATAGGCCGCG‑CTTCGGCAGTGATGATGCAAAAGGGGGTCATTCCAGGGCCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAg                                                                                                                                                                                                                          <  2:455745/222‑1 (MQ=255)
          aGAAAGCAGCGAAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCG‑CTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGACGTGTGGGGAACCACCGCC‑CAAgg                                                                                                                                                                                                                         >  2:906180/1‑239 (MQ=255)
              agcagcGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCG‑CTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGGGAACCACCGCC‑TAAGgaga                                                                                                                                                                                                                      >  2:395147/1‑238 (MQ=255)
                   cGGAGTGACATGTTTATCAACAATCGACATTTTGATATGACGTG‑CTTCGTCAGTGCTGATGCAGAAACGGGTCATTTCAGAGGCCGTACCA‑GCCGTGGTATTGCTGGCGATAATCGTCAGATTCGGTTTTG‑TCGT‑GCGGTGAACGGCTTCCAAAGCCGGGATGTCGCCGCCATCTGCAGCCCCAAAGTGC‑CTACCCTTTGCGCCTTCGTGTGgagacataagcctaaaagggaaacgga                                                                                                                                                                                                                 >  1:875051/1‑216 (MQ=255)
                    ggAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCG‑CTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCCCCATTGGCTGCCAC‑CAGCCCAATACCTTTTTCGCAGTCGTGTGGGAGACCACCGCC‑TAAGGAGATCACGCt                                                                                                                                                                                                               >  2:582421/1‑239 (MQ=255)
                      aGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCG‑CTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCCATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCtt                                                                                                                                                                                                              >  2:1106262/1‑237 (MQ=255)
                       gTGACTTTTTTAACAACAACCGCTGTTGTGATAGGCCGCG‑CTTCGTCAGTGAGGGTGGAAAAGCGGTTTTTTACAGGGGCCTACCCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTAt                                                                                                                                                                                                             <  2:516877/238‑1 (MQ=255)
                               tttacaaccatccgcttttttGATATG‑CGCGCCTTGGTCAGTATGGTTGCAGACCCGGGTTATTCAAGGGCCCTTCCCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑attat                                                                                                                                                                                                   <  2:1242715/222‑1 (MQ=255)
                                  aTCAACAATCGCCATTTTGATATGACGCG‑CTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGG‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAACGCAATACCTTTTGCGCAGTCGTGTGGCAAACCACCACC‑TACggcaatacggcattccgaatttttt                                                                                                                                                                                                  >  2:469717/1‑215 (MQ=255)
                                          tccgctttttGATAGGCGGCG‑CTTCGTCAGGTAGTATGAAGAAACGGTTAATTTCGGGGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGt                                                                                                                                                                                        <  2:1277402/233‑1 (MQ=255)
                                          tCGCCATTTTGATATGACGCG‑CTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCCC‑CAGCGCAATACCTTTTGCGCAGTCGTGGGGAGAACCACCGCC‑TAAGGAGATCACGCTATCCCA‑ATTATTCTCTTTAAGt                                                                                                                                                                                        >  2:1170260/1‑239 (MQ=255)
                                                        tGACGCG‑CTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCCCGCTATCGCA‑CTTATTCTCTTTAAGTAATTTCAAACCTgc                                                                                                                                                                          >  2:1275056/1‑239 (MQ=255)
                                                           cgcg‑CTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTgcggc                                                                                                                                                                       >  2:1200775/1‑239 (MQ=255)
                                                                cTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATTCCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTGTTGCGCAGTCGTGTGGAGAACCCCCGCC‑AAAGGGGATCACACCATCGCACATTTTCCTTTTTAAG‑AATTTAAAACCTagggcggag                                                                                                                                                                   >  1:1080883/1‑236 (MQ=255)
                                                                cTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTTTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCC‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGGGACg                                                                                                                                                                   >  2:18479/1‑239 (MQ=255)
                                                                   cGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCCCCGCC‑TAAGGAGATCACCACATCCCA‑ATTATTCTCTTTAAGGAATTTCAAACCTGCGGGGACGttt                                                                                                                                                                >  2:1114431/1‑239 (MQ=255)
                                                                     gccatgtagggtgaaaaaCGGGGCATTTGAGGGGCCCTGCCCGGCCGTGGTATGGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGtttt                                                                                                                                                               <  2:122763/224‑1 (MQ=255)
                                                                           atgatgCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCgtggt                                                                                                                                                        >  2:702732/1‑239 (MQ=255)
                                                                              atgCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCCCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTggg                                                                                                                                                     >  2:784136/1‑239 (MQ=255)
                                                                                    aaggggtgtttttaggtggcctgaCCC‑GCGGGGGTTATGGTGGCGATCACCGGCAGCTGCGTTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAgg                                                                                                                                               <  2:957044/216‑1 (MQ=255)
                                                                                      aCGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGtt                                                                                                                                             >  2:858077/1‑239 (MQ=255)
                                                                                           tatttcgagggccgcgaCCC‑GCGTTGGTATGGATGGCGATCATCGGACGCTGCGGTTTTG‑CAGG‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTg                                                                                                                                        <  2:573819/223‑1 (MQ=255)
                                                                                           cttttcgtgggggcttcccc‑GCCGTGGTTTGTATGGCGATCGCCGGGATCTGCGGTTTTG‑‑AGAGGCGGTAAACGCCGTCTTAATCCGGAATTCCGCCGCCATTGGGTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTg                                                                                                                                        <  1:18479/222‑1 (MQ=255)
                                                                                            ctttccgggggcgttACCC‑GCCGTGGTATTGATGGCGATCATCGGAAGCTGCGGTTTTG‑CAGG‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTg                                                                                                                                        <  2:1053227/224‑1 (MQ=255)
                                                                                                tCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAg                                                                                                                                                                                                                          >  2:809598/1‑153 (MQ=255)
                                                                                                            acccggtggtggtttGATGGGGATCATCGGCAGGTGCGTTTTGG‑CAGA‑GCGGTCAACGCTTTCGTAACCGCGAATAGCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAATCCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGc                                                                                                                        <  2:533244/226‑1 (MQ=255)
                                                                                                                    tGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAata                                                                                                                >  2:1171860/1‑239 (MQ=255)
                                                                                                                              ggCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGttct                                                                                                       >  2:959910/1‑238 (MQ=255)
                                                                                                                                             cTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTTCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTTAAAAATTTTCGTTCTTCCAGTGCTTTTTGc                                                                                        >  2:1298481/1‑238 (MQ=255)
                                                                                                                                               gCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGGGGAGAACCCCCGCC‑TAAGGAGATCACGCTATCGCA‑ATGATTCTCTTTAAGTAATTTCAAACCTGCGGCCACGTTTTCCGTGTGGGGGTTAGGTTGGGGGCCACCATAAAAAACCACAAAAAAATTTGGTTTTTcccagcgttttttaact                                                                                     >  1:721156/1‑224 (MQ=255)
                                                                                                                                                gggTTTTG‑CAGG‑GGGGGCAACGTCTTGATAACGGGAAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGAGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGcaca                                                                                     <  2:12342/237‑1 (MQ=255)
                                                                                                                                                               cGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATAcc                                                                       >  2:1047308/1‑239 (MQ=255)
                                                                                                                                                               cGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGCAATTTCAAACCTGCGGCGACGTTTTCCGTGGCGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCCCATAGCCCGCCATAcc                                                                       >  2:328529/1‑239 (MQ=255)
                                                                                                                                                                         cTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCCGTGCCTTTTGCACATCGCCCGCCATACCCAATTTct                                                               >  2:90947/1‑236 (MQ=255)
                                                                                                                                                                         cTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGGGCCATCCTAAATAACGCTAAAAATATTGCTTTCTTCCAGGGCTTTTTGCACATCTCCCGCCATAACTAATTTCGtt                                                             >  2:851607/1‑239 (MQ=255)
                                                                                                                                                                           ctggaacccacaatttcccctataTTGGCTCCAAC‑CAGCGCAACTCCTTTGGCGCAGTGGTGTGGAGACCCCCCGCC‑TAAGGAGATCACGAGTTCGAA‑ATTATTCTCTTTAAGTAATTTCAAACCGGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGTTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTaa                                                           <  2:34479/216‑1 (MQ=255)
                                                                                                                                                                                aTCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACatat                                                      >  2:241500/1‑239 (MQ=255)
                                                                                                                                                                                 tCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCCTTTGCGCAGTCGTGTGAAGAACCACCGCC‑TAAGGAGATCACGCTATTGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTTCGGCGACGTTTTCCTTGGTGGGTCTCGGTAGGGTGCCATCATAAATAACGTTAAAAAAATTGCGTTCTTCTAGTGCTTTTTGTACTTAGCCGCACATACCTAATTTTGTaaaattata                                                     >  1:1263247/1‑233 (MQ=255)
                                                                                                                                                                                              gcccATTGGTTGCCCC‑CAGCGCAATCCCTTTTGCGCATTCGTGTGGAAACCCCCCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCGGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGTTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAAtt                                        <  2:721156/237‑1 (MQ=255)
                                                                                                                                                                                                   gtGGCGGCCAC‑CCGCGCAATCCCTTGTGCGCTGTTGTGTGGAAACCCCCCGCC‑AAAGGAGATCACGTTATCGAA‑ATTATTCCTTTTAAGTATTGTCAACCCGGCGGCGACGTTTTCCGGGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGt                                   <  2:672703/238‑1 (MQ=255)
                                                                                                                                                                                                           cac‑caGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTaa                           >  2:592419/1‑239 (MQ=255)
                                                                                                                                                                                                                aGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAACATATTGCGTTCTTCCAGGGCTTTTTTCACATCGCCCGCCAAAACTAATTTCGTTAACCAATGTTCAGTGACAATAAAGGTACGGGAAAAAcc                       >  1:535082/1‑239 (MQ=255)
                                                                                                                                                                                                                      aTACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAAt                  >  2:652667/1‑238 (MQ=255)
                                                                                                                                                                                                                        cccTTTTGCGAAGTCGTGTGGAAACCCCCCCCC‑TAAGGAGTTCACCCTATCGCA‑ATTATTCTTTTTAAGTAATTCCAAACCTGCGGCGACGTTTCCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATTTTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTg               <  2:197255/238‑1 (MQ=255)
                                                                                                                                                                                                                         cttttttccgtgTCGTGGGGGAAACCCCCCCCTTAAGGGTATCCCGTTATCGCA‑TTTATTCTCTTTAAGTAATGACAACCCGGGGGCACCGTTTTCGGTGGGGGGTTTAGGTTGGGGTCATTCATAAATAACGATAAAATTTTTGCTTTCTTCCAGTGCTTTTTGCCCATCTCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTg               <  1:959910/228‑1 (MQ=255)
                                                                                                                                                                                                                            tttGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAACAATATTGCGTTCTTCCAGTGCTTTTTGCCCATCGCCCGCCATACCTAATTTCGTTACCATATTGTCAGGGACAATTAAGGTACGGGTaaaacaaaaaaatcccaa           >  1:197255/1‑224 (MQ=255)
                                                                                                                                                                                                                            tttGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCcat           >  2:206684/1‑239 (MQ=255)
                                                                                                                                                                                                                                 gcAGTCGGTGGGAGAACCCCCTCC‑AAAGGAGATCACGCTATCGCA‑TTTATTCTCTTTAAGTAATTTCAACCCTGCGGCGACGTTTTCCGTGGTGGGTTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCata       <  2:753528/238‑1 (MQ=255)
                                                                                                                                                                                                                                   aGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAg                                                >  2:409985/1‑195 (MQ=255)
                                                                                                                                                                                                                                    ggcggggggAAAACCACCGCT‑GAGGGAGATCCCGTTATCGCA‑ATTATTCTCTTAAAGTAATTCCAACCCGGCGGCGACGTTTCCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCa   <  2:249176/232‑1 (MQ=255)
                                                                                                                                                                                                                                    ggcGTGGGGAAACCCCCCGAC‑TAAGGAGATCAGGTTTTCGAA‑ATTTTTCTCTTAAAGTAATTTAAAACCTGGGGCGACGTTTTCCGTGGTGGGGTTAGTTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCa   <  2:138554/237‑1 (MQ=255)
                                                                                                                                                                                                                                    gTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTCACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCAAAATTTGCcaatacattaa   >  1:254644/1‑230 (MQ=255)
                                                                                                                                                                                                                                    gTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGAGCTTTTTGCACATCGCCCTCCAACCCTAATTTCGTTAACATATTTTCATCGGCAATTAAGGTACGGTGAAAACCATAAACTGAcataacataaa   >  1:677542/1‑231 (MQ=255)
                                                                                                                                                                                                                                     tCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATTTTGTCAGTGACAATTCAGGAACCGGTAAACCCAAAATATCCCATCATAATCAt  >  1:310368/1‑239 (MQ=255)
                                                                                                                                                                                                                                     tCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATATGCCATCATATTCAt  >  2:40324/1‑239 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                       
AGTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACGCG‑CTTCGTCAGTGATGATGCAGAAACGGGTCATTTCAGAGGCCGTACCC‑GCCGTGGTATTGATGGCGATCATCGGCAGCTGCGGTTTTG‑CAGA‑GCGGTCAACGCCTTCGTAATCGCGAATATCGCCGCCATTGGCTGCCAC‑CAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCC‑TAAGGAGATCACGCTATCGCA‑ATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCAT  >  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: ‑