Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,786,444 C→T I395I (ATC→ATT)  gpmM → 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,786,4440CT95.6% 152.3 / NA 45I395I (ATC→ATT) gpmM2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/21);  minor base G (0/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGTTAACACCTTCGGCGAGTGGATGGCGAAAAACGACAAAACTCAGTTGCGTATTTCCGAAACCGAAAAATATGCCCACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGC‑GAAAGCGGTTG‑AGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGG‑TA‑ACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGAT  >  NC_000913/3786198‑3786684
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
ggTTAACACCTTCGGCGAGTGGATGGCGAAAAACGACAAAACTCAGTTGCGTATTTCCGAAACCGAAAAATATGCCCACGTTACTTTCTTCTGCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCAGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTTACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCt                                                                                                                                                                                                                                                   <  1:205816/251‑1 (MQ=255)
   tAACACCGTCGGCGAGTGGATGGCGAAAAACGACAAAACTCAGTTGCGTATTTCCGAAACCGAAAAATATGCCCACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGt                                                                                                                                                                                                                                                 <  2:56504/250‑1 (MQ=255)
           gcGGCGAATGGATGGCTAAAAACGACAAAACTCAGTTGCGTATTTCCGATTTCGAAAAATGAGCCCACGATACTTTATCCGGCAACGGTGGCGTAGAAGAGTCGTTCAAAGCCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTTCCTACGATCTGCAACCTGAAATTAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCGATAGCGACAACTACGACACCATTATCTGTTTCTATcc                                                                                                                                                                                                                                         <  2:101398/249‑1 (MQ=255)
                     gATGGCGAAAAACGACAAAACTCAGTTGCGTATTTCCGAAACCGAAAAATATGCCCACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACa                                                                                                                                                                                                                              <  1:10273/251‑1 (MQ=255)
                       tGGCGCAAAATGACACAACTCAGTTGCGTATTTCTGAACCCGAAAAATATGCCCCCGTTACTTTCTTCTTCAACGGTGGCGTAGAAGCGTAGTTCAAAGGCGAAGATCGCATTCTGATCCACTCACCGAAAGTGGCTACCTCCGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATAAGACACCATTATCTGTAACTATCCGAACGGTGACATg                                                                                                                                                                                                                            >  2:286466/1‑251 (MQ=255)
                          cGAAAAACGACAAAACTCAGTTGCGTATTTCCGAAACCGAAAAATATGCCCACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGt                                                                                                                                                                                                                          >  1:84200/1‑250 (MQ=255)
                               aaCGACAAAACTCAGTTGCGTATTTCCGAAACCGAAAAATATGCCCACGTTACTTTCTTCTTCAAAGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGGACGGGGACACGGTAGGTc                                                                                                                                                                                                                     >  2:18440/1‑250 (MQ=255)
                                                 cGTATTTCCGAAACCGAAAAATATGCCCACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAgc                                                                                                                                                                                                  <  2:16879/251‑1 (MQ=255)
                                                        ccgaaaccaaaaAATATTCCCACGTTACTTTCTTCTTCAACTGTGGCGTAGAAGAGTCGTGCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTTCGCAGAGCTGACCGAAAAACGGGGTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGGAACTATCCGAACGTTGACATGGTAGGTCACACCGGTGTGATGGAAGCGGCGGt                                                                                                                                                                                            <  2:358070/250‑1 (MQ=255)
                                                          gaaaccgaaaAATATGCCCACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTa                                                                                                                                                                                          >  2:424103/1‑250 (MQ=255)
                                                                  aaaaTATGCCCACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTg                                                                                                                                                                                >  1:302189/1‑251 (MQ=255)
                                                                       atGCCCACGTTACATTCTTCTTCAACGGTGGCGGAGAAGAGTCGTGCAAAGTCGAAGATCGCGTTCTGATCAACCCACCGAAAGTGGCTACCTACGAGCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCAGCAAGAGCGGCAAATACGACACCATGATCTGTAACTACCCGAACGTTGACATGGTAGGTCACCCCGGGGTGATGGGCCCGGCGGTTAAAG‑CGGTTTAAgca                                                                                                                                                                           <  2:248748/251‑2 (MQ=255)
                                                                       aggtccACGTTACTTCCGTCTTCAACGGTTGCGTAGAAGAGTCTTTCAAAGGCGAAGATCGCATTCTGAGCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGTCAAATACGACACCATTATCTGTAACTATCCGAACGGTGAGATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTGAAAG‑CGGTTGacgga                                                                                                                                                                           <  2:424434/247‑5 (MQ=255)
                                                                           ccACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTaaaaagggtgaagaactgg                                                                                                                                                                        >  1:255655/1‑235 (MQ=255)
                                                                             aCGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTGAAGCACTGGAt                                                                                                                                                                     <  1:383157/251‑1 (MQ=255)
                                                                                  aCTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATGATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTGAAGCACTGGATCACTg                                                                                                                                                                <  1:286466/251‑1 (MQ=255)
                                                                                   cTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTGAAGCACTg                                                                                                                                                                        <  2:39692/242‑1 (MQ=255)
                                                                                   cTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGGTGAAGCACTgagatcggg                                                                                                                                                                >  1:39692/1‑242 (MQ=255)
                                                                                    tttCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAA‑CGGGTGAAGCACTGGAACACcgcg                                                                                                                                                              >  1:249538/1‑251 (MQ=255)
                                                                                        ttcttcAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTCTCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGGTAAAG‑AGGTTggaacaacgggttactgcggggg                                                                                                                                                          >  1:238199/1‑229 (MQ=255)
                                                                                        ttcttcAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGGCACACCGGGGTGCTGGAAGCGGCGGTTAAAG‑CGGGTGaaagactggagaactgcgtggg                                                                                                                                                          >  1:103653/1‑231 (MQ=255)
                                                                                         tcttcAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTGAAGCACTGGATCACTGCGTg                                                                                                                                                            <  1:413721/248‑1 (MQ=255)
                                                                                          cttcAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGGTAAAG‑AGGGTGAAGCGCttgaactatgcgtggg                                                                                                                                                          >  2:187612/1‑234 (MQ=255)
                                                                                                     gCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGGTGGAGGACTTGGTCACCGCGTGGAAGAAGTCGG‑Gaa                                                                                                                                             >  1:83782/1‑250 (MQ=255)
                                                                                                       gTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTGAAGCACTGGATCACTGCGGTGAAGAAGTCGC‑GAAAGc                                                                                                                                          <  1:132391/251‑1 (MQ=255)
                                                                                                                        aGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑AGGGTGGAGCACTGGATCACTTCGTGGGAAAAGTCGC‑GAAAGAGGATG‑AGTCCgtggggg                                                                                                                        >  1:431999/1‑249 (MQ=255)
                                                                                                                                             cAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCGCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACAACGGGGTGATGGAAGCGGCGGTTAAAG‑CGTTTGAAGCACTGGGTCACTGCGTGGAAGAATTCGC‑GAAAGCGGTTG‑AGTCCGTGGGTGGACAACTGCTGATCACCGCt                                                                                                    <  2:430556/250‑1 (MQ=255)
                                                                                                                                               aCTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGGTAAAG‑CGGTTGGAACACTGGGTCACTGCGGGGAAGAAGACGC‑GAAAGCGGGTG‑GGTCCGTGGGGGGACAAATGCCGATCACCCCTgcc                                                                                                 >  1:80870/1‑249 (MQ=255)
                                                                                                                                                     cGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTGAAGCACTGGATCACTGCGTGGAAGAAGTCGC‑GAAAGCGGTTG‑AGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGG‑t                                                                                          <  2:171238/251‑1 (MQ=255)
                                                                                                                                                      gAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTGAAGCACTGGATCACTGCGTGGAAGAAGTCGC‑GAAAGCGGTTG‑AGTCCGTGGGTGGACAACTGCTGATc                                                                                                          <  1:433248/235‑1 (MQ=255)
                                                                                                                                                       aaaGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTGAAGCACTGGATCACTGCGTGGAAGAAGTCGC‑GAAAGCGGTTG‑AGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGG‑Ta‑a                                                                                       <  1:319441/251‑1 (MQ=255)
                                                                                                                                                              ctacctacGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGGTGAAGCACCGGGGCACTGCGTGGGAGAAGGCGC‑GGAAGGGGGTG‑AGGCCGGGGGGGGGCAAACGGTGATCACCGCTGGCCACGG‑GA‑ACGCTgg                                                                                 >  2:75628/1‑249 (MQ=255)
                                                                                                                                                               tacctacGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAA‑CGGTTGTAGCACTCGGTCACTCAGTGGAGAAAGTGGC‑GAAAGAGGGTGTAGTC‑GGGGGGGGGCAACAGCTGATCACCGCCGCACACGG‑TA‑ACGCTGaa                                                                                >  2:14192/1‑249 (MQ=255)
                                                                                                                                                                acctacGATCTGCAACCGGAAATGAGCCCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTGAAGCACTGGATCACTGCGTGGAAGAAGTCGC‑GAAAGCGGTTG‑AGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGG‑TA‑ACGCTGaac                                                                               <  1:382749/250‑3 (MQ=255)
                                                                                                                                                                 cctacGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTGAAGCACTGGATCACTGCGTGGAAGAAGTCGC‑GAAAGCGGTTG‑AGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGG‑TA‑ACGCTGAACAg                                                                             <  1:187612/251‑1 (MQ=255)
                                                                                                                                                                  ctacGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAA‑CGGGTGGAGCACTGGATCACTGCGTGGACGAAGTCGC‑GCAAGCGGGTG‑AGGGCGGGGGTGGACAAATGCCGGTCCCCGCCCACCCAGG‑TA‑ACGCTGaaacg                                                                             >  2:118169/1‑246 (MQ=255)
                                                                                                                                                                             aaCCGGAAATGAGCTCGGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAAACCGACACCTTTAGCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTGGAAGCACTGGATCACTGCGTGGAAGAAGTCGC‑GATAGCGGTTG‑AGTCCGTGGGTGGACATCTGCTGATCACCGCTGTCCACGG‑TA‑ACGCTGATCAGATGCGCGATCCg                                                                 <  2:434434/251‑1 (MQ=255)
                                                                                                                                                                                                          gAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTGAAGCACTGGATCACTGCGTGGAAGAAGTCGC‑GAAAGCGGTTG‑AGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGG‑TA‑ACGCTGAACAGATGCGCGATCCGGCAACGGGTCAggcacacacggcacacac                                    <  2:83782/251‑1 (MQ=255)
                                                                                                                                                                                                                         gCCATCAAGAGCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTGAAGCACTGGATCACTGCGTGGAAGAAGTCGC‑GAAAGCGGTTG‑AGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGG‑TA‑ACGCTGAACAGAGGCGCGATCCGGCAACGGGTCAGGCACACACGTCACACACCAACCTGCCAGTTcc                     <  2:84200/251‑1 (MQ=255)
                                                                                                                                                                                                                           cATCAAGAGCGGCAAATACGACACCATTATCTGAAACTATCCGAACCGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGAGGGTTGAA‑CACTTGATAACTGCGGGGGAACAGTGGG‑GACAGGGGTTG‑AGTCCGTGGGGGGAAAAATGATGGTTCCCGCTGAACACGG‑GA‑ACGCCGGACAGAAGCGCGCTCAGGCAAAGGGGCAGGGACACACGGCACACACCACCCCGCCAGTTCCGc                   >  1:276338/1‑251 (MQ=255)
                                                                                                                                                                                                                                agagCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTTAAG‑CGGTTGAAGCACTGGATCACTGCGTGGGAGACGTCGC‑GCACAAGGTTG‑GGTCCGTGGGGGGACAACCGCTGATCACAGCTGAACACGG‑GACAC‑CTGCAAAGCAGCGCGGTCCCGCAACGGGGCAGGCACACACGGCACACACCACCCtacccggtcccctggtt              >  2:391933/1‑235 (MQ=255)
                                                                                                                                                                                                                                agagCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTGAAGCACTGGATCACTGCGTGGAAGAAGTCGC‑GAAAGCGGTTG‑AGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGG‑TA‑ACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGAtt              <  2:302189/251‑1 (MQ=255)
                                                                                                                                                                                                                                agagCGGCAAATACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGGTGTAGCACTGGATCACTTCGTGGGAGAAGTCGC‑GCAAGCGGGTG‑AGTCCGTGGGTGGGGAACTGCTGGTCACCGCCGACCACGG‑GA‑ACGCTGGACAGATGGGCGGTCCGGGCACGGGGCAggcacacaccgcacacacacacccgccccttcccccgatt              >  2:221403/1‑229 (MQ=255)
                                                                                                                                                                                                                                         aaTACGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTTAAG‑CGGGTGAAGCACTGGATCACTGCGGGGGAGAAGTCGC‑GGAAGAGGTTG‑AGTCCGTGGGGGGGCAACTGCTGATCACCGCTGACCACGG‑TA‑ACGCTGGAAAGATGCGCGCTCCGGCAACGGGGCAGGCACACACGGCACCCACCAACCCGCCAGTTCCGCTGATTTTCGtttgt     >  1:102222/1‑248 (MQ=255)
                                                                                                                                                                                                                                             cGACACCATTATCTGTAACTATCCGAACGGTGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGGTAAAG‑CGGATGTAGCCCTGGATCCCTGCGTGGGAGCAGGCGCAGAAA‑AGGGTG‑AGTACGCGGGTGGGAAAATGCTGATACCCGCTGACCACGGGTA‑A‑GCTGCACAGGTGCGCGGTCCGGCCACGGGGCCAGCACACACCGCAACCACCACCATGCCACCTCCCCTGATTTACGtgtgtggt  >  2:224347/1‑243 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
GGTTAACACCTTCGGCGAGTGGATGGCGAAAAACGACAAAACTCAGTTGCGTATTTCCGAAACCGAAAAATATGCCCACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAG‑CGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGC‑GAAAGCGGTTG‑AGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGG‑TA‑ACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGAT  >  NC_000913/3786198‑3786684

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