Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,139,288 C→T D141N (GAC→AAC)  fsaB ← fructose‑6‑phosphate aldolase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,139,2880CT98.4% 211.6 / NA 64D141N (GAC→AAC) fsaBfructose‑6‑phosphate aldolase 2
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (37/25);  minor base G (0/1);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGAGATGAGTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTC  >  NC_000913/4139047‑4139532
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
tagagtagAGTTGTGCCAAATGCGGCATTCCAGTCGTGGTCGGACTTCTTTAGAGCTGACTCTCCCGCAGGGTTTTTGGGGCTTTGTTGCGCAAATTCTAAGGGCAGGGTGATGGATTCGCATCCTGCCAGTAAACAGTCCATCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTGACAGCGTTTGCATCTCCGGAACCGTACGAATGCCGTTTCCGCCCt                                                                                                                                                                                                                                              <  2:16755/243‑1 (MQ=255)
tAGAGATGAGTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCt                                                                                                                                                                                                                                              <  1:190875/250‑1 (MQ=255)
tAGAGATGAGTAGTGCCAAATGCGGCATTCCAGTCGAGTTCGAACTTCTCTAGAGCTGACTCTACAGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGCTTCACATCCTACCAGCAACCAGTCCAGCGCCTGACGCGGCGTTTTAAAACTGGCTGCCAGCAACATGCTTTCTGGCGCGTGCCATTCTACCAGCGTTTGCAGCTCCTGAACAGTTCGAATGCCGTCTCCGCCCtt                                                                                                                                                                                                                                             >  2:125596/1‑250 (MQ=255)
           agtGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTAcgc                                                                                                                                                                                                                                  >  1:463590/1‑251 (MQ=255)
              gCCAAATGCGGCCTTCCAGTCTTGTTCGATCTTCTCTATAGCTGACTCTCCCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCCGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGt                                                                                                                                                                                                                               <  1:379653/251‑1 (MQ=255)
                 acaTGCGGCCTTCCCGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCTCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACCGCGTTTGCAGCCCCTGAACCGTACGACTGCCGTTTCCGCCCTTGGCATCTACGCGGTTa                                                                                                                                                                                                                             >  2:227649/3‑250 (MQ=255)
                        gCATTACAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCAAGCGCCTGACGCGGCGTTTTAAAGATGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCAGCCCTGGGCATCTACGCGGTTAACATACgg                                                                                                                                                                                                                     >  1:417014/1‑251 (MQ=255)
                          aTTCCAGTCGTGTTCGAACGTCTCTATAGCTGACTCTACCGCAGGGTTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAg                                                                                                                                                                                                                   <  1:332711/251‑1 (MQ=255)
                              cAGTCGTGTTCGAACGTCTCTATAGCGGCCTCGACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAAc                                                                                                                                                                                                               <  2:233141/251‑1 (MQ=255)
                               aGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACg                                                                                                                                                                                                              >  1:434046/1‑251 (MQ=255)
                                gTCGTGTTCGACCTTCTCTATAGCTGACTCAACCGCCGGGGTGTTGAGCATTTGTTGCTCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACgg                                                                                                                                                                                                             >  2:221650/1‑250 (MQ=255)
                                gTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACgg                                                                                                                                                                                                             >  2:496272/1‑250 (MQ=255)
                                gTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGt                                                                                                                                                                                                             >  1:182975/1‑251 (MQ=255)
                                  cGTGTTCGATCTTCTCTATAGCTTACTCTACCGCAGGGTTGTTGAGCATTTGTTGCGCTACATCGAAGGGCAGGGTGATGGATTCGCATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTAt                                                                                                                                                                                                           <  1:249727/251‑1 (MQ=255)
                                     gtTCGACCTTCACTATAGCTGACTCTACCGCAGGGGGGTTGAGCATTTGTTGCGCTACATCTAAAGGCATGGAGATGGATTCACATCCTGCCAGTAAACAGTCCAGAGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCGGGCGCGTGCATTTCTAACAGCAATTGCAGCTCCTGAACCGTACGCATGCCGTTTCCGCCCTGGCCCTNAACGCGGTTGACATACgacgacccgaatttt                                                                                                                                                                                                        >  2:252494/1‑237 (MQ=255)
                                     gtTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTAtttt                                                                                                                                                                                                        <  1:390654/251‑1 (MQ=255)
                                              tctcTATAGCTGACTCTACCGAAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGACACGTATTTTCCCCCTGcc                                                                                                                                                                                               >  2:141505/1‑251 (MQ=255)
                                                   aacgCTGACTCTCCGTCAGGGTTGTTGCGCATTGTTTGCGCTCCATCAAAGGGCAGGGGGATGGTTTCACACCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGATTGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGAAACGTATTTTGCCCCTGCCAGt                                                                                                                                                                                            <  2:463590/246‑1 (MQ=255)
                                                        tGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCt                                                                                                                                                                                      >  1:332517/1‑250 (MQ=255)
                                                             ctcccGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTCAGGGCAGGGTGATGTATGCACATCCGGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCTTGCATTTCTAACATCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCtaata                                                                                                                                                                                  <  2:325920/246‑1 (MQ=255)
                                                                ccGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCt                                                                                                                                                                             <  1:265179/251‑1 (MQ=255)
                                                                    cAGGGGGGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGAGATGGATTCACATCCTGCCCGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTTCAGCTCCTGACACGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCATGCCATGGCGGCTAATAACCCTtgt                                                                                                                                                                          >  2:111439/1‑251 (MQ=255)
                                                                     aGGGGTGTTGAGCATTTGTTGCGCTACATATAAGGGCAGGGTGATGGATTCACATCCTGCCAGAAAACAGTCCAGCGCCTGACGCGGAGTTTTAACGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGCACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTACCATACGGAGCAACGTATTTTGCCCCTGCCAGGGCGGTTAATAACCCTtgtg                                                                                                                                                                         >  2:340360/1‑251 (MQ=255)
                                                                          ttttGAGCATTTGTTGAGCTACAACTAAGGGCAGGGTTATGGATTCCCCCCCGTCCAGTAAACATGCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACTAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATATCCCTTGTgcgga                                                                                                                                                                     <  2:31519/248‑2 (MQ=255)
                                                                           gtTGAGCTTTTCGTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGAAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGAATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGgcgc                                                                                                                                                                   >  2:334038/1‑251 (MQ=255)
                                                                              gAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCta                                                                                                                                                                 <  1:20536/250‑1 (MQ=255)
                                                                                     tgttgCGCTCCATCTAAGGGCAGGGTGAGGGATTCACATCCTGCCAGTAAACAGTCAAGCGCCTGACGCGGCGTTTTAAAGCTGGCTTCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACTGCGTTTGCAGCTCCTGAACAGTACGAAAGCCGTTTCCGCTCATGGCATCAACGCGGTTACCATCCGGATCAACGAATTTTGCCCCTGCCAGTGAGGCTAATCAACATTGTGCGGCGCACTATAcagc                                                                                                                                                         >  2:110942/1‑251 (MQ=255)
                                                                                          cgcTCCATCTAGGGGCAGGGTGAGGGATTCACCTCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGc                                                                                                                                                    <  2:182975/251‑1 (MQ=255)
                                                                                          cgcTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGc                                                                                                                                                    >  1:520420/1‑251 (MQ=255)
                                                                                             tACATCTAAGAGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTCCGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGCATTTTGCCCCTGCCAGCGCGGCTAATAAACCTTGTGAGGCGCTCGATACAGAAGTGCCa                                                                                                                                                  >  2:63565/1‑250 (MQ=255)
                                                                                                aTCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATCACCCTTGGGCGGCGCTATTTACCGCCGTGCCAAg                                                                                                                                                >  2:543250/1‑249 (MQ=255)
                                                                                                  cgaaGGGCGGGGTGTTGTATCCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCGTGTGCGGCGCTATATACAGCAGTGCCAAgtgt                                                                                                                                             <  2:565059/248‑1 (MQ=255)
                                                                                                    aaGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAgtgt                                                                                                                                             <  1:536565/248‑1 (MQ=255)
                                                                                                     aGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGGGCGGCGCTATATACAGCAGTGCCAAGTgtagt                                                                                                                                          >  2:370292/1‑250 (MQ=255)
                                                                                                      gggCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTgtagta                                                                                                                                         >  1:98375/1‑250 (MQ=255)
                                                                                                      gggCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTaa                                                                                                                                        <  1:258408/251‑1 (MQ=255)
                                                                                                       agCAGGGGGATGAATTCACATCCTGCCAGTAAACCGTACAACGCCTGACGCGGCGTTTTCAAGCCGGCTGCCAGCACCCTGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCATCTCCTGACCCGTACGAATGCCGTTTCCGCACTGGGCATCTTCGCGGTTAACATACGGCACAACGGATTTTACCACGGCCAGCGCGCCTTATACCCCTTGTGGGNTGCTAGATACAGCAGTGCCAAGTgtagaaa                                                                                                                                        >  2:274078/2‑247 (MQ=255)
                                                                                                          aGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATAcc                                                                                                                                    >  1:171018/1‑251 (MQ=255)
                                                                                                               gATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCTAGTgta                                                                                                                                            >  2:428984/1‑238 (MQ=255)
                                                                                                               gATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTgta                                                                                                                                            <  1:428984/238‑1 (MQ=255)
                                                                                                                 gggATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCGTTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCtttt                                                                                                                             <  1:375585/250‑1 (MQ=255)
                                                                                                                  ggATTCACATCCGTCCAGTAAACAGTCCAGCGCCTGAGGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTTTAGTAATACCCTCttttt                                                                                                                            <  2:492676/251‑1 (MQ=255)
                                                                                                                   gATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCtttttt                                                                                                                           >  1:272174/1‑251 (MQ=255)
                                                                                                                     ttCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGAGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGCAGTAATACCCTCTTTTTTc                                                                                                                          >  1:499121/1‑250 (MQ=255)
                                                                                                                        acaTCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGATTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACAGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGAATTTTGCCCCTGACAGAGCGGCTAATAACCCTTGTGCGGCGCTAAATACAGCAGTGCCAAGTGTAGTAAAACCCTCTTTTTTCAGt                                                                                                                       >  1:566552/1‑250 (MQ=255)
                                                                                                                                         tgTCCAGCGCCTGACACGGCGTTTTAAAGCTAGCTGCCAGCACCATGCTTTCTGGCGCCTGCATTTCTAACAGAGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCCACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGAAGCAATCCCCTCTTTTTTAAGTATTTTACTTGTTGccggc                                                                                                     >  2:228173/2‑248 (MQ=255)
                                                                                                                                            ccAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCTCGTGCATTACTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTGTCCGCCCTGGGCATCTATGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGATAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTTCCAGACCt                                                                                                  <  2:141714/251‑1 (MQ=255)
                                                                                                                                                cgcCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTGTAATTGCTGCCAGACCTTCgg                                                                                              <  1:354131/251‑1 (MQ=255)
                                                                                                                                                         cggcgTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCg                                                                                     >  2:95672/1‑251 (MQ=255)
                                                                                                                                                              ttttAAAGCTGGCTGCCATCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAACCCCTGAACCGTACGAATGCCGTTTCCGCCCTCTGCAACTACCCGCATAACATACGAAGCACCGTATTTTGCCCCCGCCAGTGCGGCTAATACCCCTTGTCCCGCGCTATATACAGCAGTGCCCAGGGTCGTAATAACCTCTTATTTCAGTATTTTACTTTCTTCCAGTCACTCGGAAGTCACCggactt                                                                                >  2:474675/1‑247 (MQ=255)
                                                                                                                                                                 tAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTc                                                                             >  1:20500/1‑251 (MQ=255)
                                                                                                                                                                      ctggctgCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGGGCGGCTAATAACCTTTGTGCGGCGCTATATACGGCAGTGCCAAGTGTAGTAATACCCTCTTGTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAATTCCCCGGGTTTTTcaccac                                                                        >  2:512490/1‑251 (MQ=255)
                                                                                                                                                                         ggtgCCAGCACCATGGTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCTTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGATCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCGATATACAGCAGTGCCAAGTGTAGTAAGGCCCTCTTTTTTCAGTAGTTTAATTGCGGCCGGACCTTCGGAAGTCACCGGGATTGTCACCACAAt                                                                     <  2:157107/249‑1 (MQ=255)
                                                                                                                                                                                       gCTTTCGGGCGCGTGAAATTCTAACGGCGGTAGCGGCTCCTGAACCGTACGAGTTGCGTTTCCGCCCTGGGTATCTACGCGTGTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATAACCGCTTTTTTCAGTATTTTACTTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCTCCACAATACCCGGAATAGCGt                                                       <  2:500395/251‑1 (MQ=255)
                                                                                                                                                                                        cTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGCTAACATACGAAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGTTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACCATACCCTGACTAGCGTc                                                      >  2:504266/1‑251 (MQ=255)
                                                                                                                                                                                                     gCATTTCAAACAGCGGTTGCAGCTCCTGAACCGAACGAATGCCGTTTACGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGAGCGGCTAATACCCCTTGTGCGGCGCTATATACACCAGTGCCAAGTGTAGTCAAACCCTCTTTTTCCAGTATTTTAGTTTCTGCAAGACCTCCGGAGTTCACCGGGGTTTTCCCCACAATACCCGGAATAgagtgccgcacgagcttc                                         >  2:413526/1‑234 (MQ=255)
                                                                                                                                                                                                     cgattacgcaCAGAGTTTGCAGCACCTGAACCGTACGAATGCCGTTTCCGCCCGGGGCTTCTACGCGGTTAACATTCGGAGCAACTTATTTTGCCCCTGCCAGTGCGGCTAATATCCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGGAATACCCTCTTTTTTCAGTATGTTAATTGCTGCCAGACCGTCGGAAGTCACCGGGATGTTCACCACAATTCCCGGAATAGGGTGGCGTATGcgcttc                                         <  2:33441/242‑1 (MQ=255)
                                                                                                                                                                                                           cTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGcttct                                   >  1:127512/1‑251 (MQ=255)
                                                                                                                                                                                                                   gTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGCTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCACCGCTATATACAGCAGAGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTCCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATATCGTCGCGCGGGCGCTTCGCTTCTTCCACCAt                           >  1:355031/1‑251 (MQ=255)
                                                                                                                                                                                                                         aGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCACCGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATccctgc                     >  2:413042/1‑248 (MQ=255)
                                                                                                                                                                                                                           cTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTgcgc                  >  1:408314/1‑251 (MQ=255)
                                                                                                                                                                                                                             ccTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTgcgc                  <  1:221650/249‑1 (MQ=255)
                                                                                                                                                                                                                                        cGAATGCCGTTTCCGCCCTGGGCATCTACGCGGGTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCCCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATg     >  2:389166/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aTGCCGTTTCCGCCCGGGGCATTTACGCGGTTAACAAACGGAGCAACGTAGTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACATCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATATCGTCGCGCAGGCGCTGCGCTTCTTCCATCATCCCTGCGCGTCGCGGCTCATGGTc  <  2:350998/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
TAGAGATGAGTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTC  >  NC_000913/4139047‑4139532

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