Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,336,411 2 bp→AT coding (83‑84/1260 nt) murA ← UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,336,4110GA90.0% 165.5 / 8.9 63I28N (ATC→AATmurAUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base A (27/27);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01
*NC_0009133,336,4120AT90.3% 173.2 / 8.8 63I28N (ATC→AAT) murAUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base T (30/26);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

GCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGAC  >  NC_000913/3336173‑3336652
                                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                                
gCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGagc                                                                                                                                                                                                                                        <  1:219536/250‑1 (MQ=255)
 ccccTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGCTAGAAGCACGCATGGGTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGGCGTCTTTCAGTTTCGGGACGGTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCCAAAAGATTCGGCAGagccg                                                                                                                                                                                                                                      >  2:559295/1‑249 (MQ=255)
 ccccTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGagcag                                                                                                                                                                                                                                      <  1:250605/251‑1 (MQ=255)
       aCCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCAtttt                                                                                                                                                                                                                                >  2:236844/1‑251 (MQ=255)
         cAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTa                                                                                                                                                                                                                              <  1:161897/251‑1 (MQ=255)
              cgcgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCTCAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGAGCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGcc                                                                                                                                                                                                                         <  2:93827/251‑1 (MQ=255)
                     ccgggggCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATGCATTAACGTCGCGGGCATCAATATTCACATAACCATAACGTTCTACTTTCGCCCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAAtg                                                                                                                                                                                                                  <  2:577183/246‑2 (MQ=255)
                          ggCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAACGATTAGGCAGAGCAGCATTTTTCGCGCCGGAAATTGTGAc                                                                                                                                                                                                             >  2:366898/1‑251 (MQ=255)
                           gCCCCAGCGCCCAAATAGAAGCACTCATTGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTACCGTCGCGGTCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACt                                                                                                                                                                                                            <  2:513204/251‑1 (MQ=255)
                                  cgcCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCt                                                                                                                                                                                                     >  2:327059/1‑251 (MQ=255)
                                  cgcCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCt                                                                                                                                                                                                     <  1:366898/251‑1 (MQ=255)
                                  cgcCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCt                                                                                                                                                                                                     <  1:495428/251‑1 (MQ=255)
                                             aaGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACACAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGtt                                                                                                                                                                                          <  1:412195/251‑1 (MQ=255)
                                                cACGCATGGTTTTAACCAGATCGTAAGGGGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCGGGAGCTTCGTTGGc                                                                                                                                                                                       >  2:39226/1‑251 (MQ=255)
                                                    aaTGGTTTTAACCAGAGCGTAAGGTGCGCAGAATACATAAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCGGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCt                                                                                                                                                                                   <  2:408231/250‑1 (MQ=255)
                                                           ttAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAacac                                                                                                                                                                             <  1:327059/250‑1 (MQ=255)
                                                                   aTCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATtg                                                                                                                                                                                                                 >  1:42369/1‑206 (MQ=255)
                                                                   aTCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTCCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATtg                                                                                                                                                                                                                 <  2:42369/206‑1 (MQ=255)
                                                                     cGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATCTGCACAGAACCATTACTTTCTCCTTGCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGTTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGTCTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATc                                                                                                                                                                   <  2:588246/250‑1 (MQ=255)
                                                                      gTAAGGTGCGCAGAATACCTTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCACACCCAGCTGGCTTAGCAGCTTAATTGCTGTATCGACGTCTTTCAGTTTCAGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGACATTGTGCCTTCGACCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATca                                                                                                                                                                  >  2:648365/1‑249 (MQ=255)
                                                                              cgcAGAATTCAGTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTa                                                                                                                                                                                                                              <  2:648047/182‑1 (MQ=255)
                                                                              cgcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTa                                                                                                                                                                                                                              >  1:648047/1‑182 (MQ=255)
                                                                              cgcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGttt                                                                                                                                                         <  1:388768/251‑1 (MQ=255)
                                                                              cgcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTATACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGAAGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGCAATTTATCCATTTCGttt                                                                                                                                                         >  1:417481/1‑251 (MQ=255)
                                                                                            aCGTCGCGGGCAGCAATCTGCACAGAACCATTACGTTCTACTTTCGCACTCAGCTGGGTTAGCAGCTTCATTGATGTATCGACTTCTTGCAGTTGCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCATTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAAGTGTGACTTCGCCCTGTAGCTGCGTTGGCCCCTGAACACGAAATTTATCCATTTAGGTTGTTCTCAGTTAACa                                                                                                                                           <  1:540353/251‑1 (MQ=255)
                                                                                               tCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGCCGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGCTTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGACCACGAATTTTATCCATTTAGTTTGTTCTCAGTTAACAAtt                                                                                                                                        >  2:181387/1‑251 (MQ=255)
                                                                                               gcgcgGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGGTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAAtt                                                                                                                                        <  1:335010/250‑1 (MQ=255)
                                                                                                cgcgGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAAtt                                                                                                                                        <  1:296970/250‑1 (MQ=255)
                                                                                                                 acaGAACCATTACGTTCTACTTTCGCTCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCCGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCCTTTTTAGCGCCGGACATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGATATTTATCCATTTAGTTTGTTCTCAGTTACCAATTCATATCCGCTACCGGCg                                                                                                                       >  2:299003/1‑250 (MQ=255)
                                                                                                                 acaGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGa                                                                                                                      <  1:371142/251‑1 (MQ=255)
                                                                                                                     aaCCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTtcgca                                                                                                                                                  <  1:506331/219‑4 (MQ=255)
                                                                                                                     aaCCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCa                                                                                                                                                  >  2:506331/1‑219 (MQ=255)
                                                                                                                      aCCATTACTTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGGATCGACGTCTTGCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAacac                                                                                                                                                                             <  2:509411/191‑1 (MQ=255)
                                                                                                                      aCCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAacac                                                                                                                                                                             >  1:509411/1‑191 (MQ=255)
                                                                                                                            aCGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCAAATCGCCCAt                                                                                                             >  2:615885/1‑249 (MQ=255)
                                                                                                                                   aCTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCaaa                                                                                                     >  1:337574/1‑250 (MQ=255)
                                                                                                                                        cGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGt                                                                                               >  1:130067/1‑251 (MQ=255)
                                                                                                                                            cccAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCGTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAg                                                                                           <  1:353454/251‑1 (MQ=255)
                                                                                                                                              cAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTCGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGAGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGACATTTATCCATTTAGTGTGTTCTCAGTTAACATTTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGtt                                                                                         >  2:279016/1‑251 (MQ=255)
                                                                                                                                                   ggCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGAGAGGCAGAGCAGTATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTGCTCAGTTAACAATTAATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTCCAGTTTGc                                                                                      <  1:588611/249‑1 (MQ=255)
                                                                                                                                                     cTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTAACGGCGAAACGCCCATAGCTCAAAAGCCGTTCAGTTTGCGacc                                                                                  >  2:538952/1‑249 (MQ=255)
                                                                                                                                                      ttAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATcg                                                                                 >  1:116103/1‑251 (MQ=255)
                                                                                                                                                        agcagcTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTg                                                                                                                                                        >  2:282370/1‑178 (MQ=255)
                                                                                                                                                        agcagcTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTg                                                                                                                                                        <  1:282370/178‑1 (MQ=255)
                                                                                                                                                              ttCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGccc                                                                          >  1:255054/1‑250 (MQ=255)
                                                                                                                                                                    gATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGc                                                                   >  1:348209/1‑251 (MQ=255)
                                                                                                                                                                     aTGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCg                                                                                                                          >  1:435655/1‑195 (MQ=255)
                                                                                                                                                                     aTGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCg                                                                                                                          <  2:435655/195‑1 (MQ=255)
                                                                                                                                                                         aTCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAgggg                                                              >  1:645227/1‑251 (MQ=255)
                                                                                                                                                                            gACGTCTTTCATTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGtat                                                           <  2:255054/251‑1 (MQ=255)
                                                                                                                                                                             aCGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGtata                                                          <  2:116103/251‑1 (MQ=255)
                                                                                                                                                                                                 gTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTCCCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCAt                                      >  2:692073/1‑251 (MQ=255)
                                                                                                                                                                                                    cTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCACCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGa                                    >  1:441043/1‑250 (MQ=255)
                                                                                                                                                                                                      ggATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCTGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAAt                                  <  1:279016/250‑1 (MQ=255)
                                                                                                                                                                                                         tctcTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCgg                              >  2:291842/1‑251 (MQ=255)
                                                                                                                                                                                                         cctctACCGTGTCGTCCTCCAGTAGTGCGGCAAAAAGAGTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTGCGCTCAGGAGGTTCGTTTGCCCCGTACCACGAAATTGATCCATTTATTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCgg                              <  2:355248/250‑1 (MQ=255)
                                                                                                                                                                                                              aCCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTa                           <  1:615885/249‑1 (MQ=255)
                                                                                                                                                                                                              aCCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGACAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATc                         >  1:601460/1‑251 (MQ=255)
                                                                                                                                                                                                                            cAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGACGGAAATTGTGACTTAGCCCTGGAGCTTCGTTGGCCCCTGCACACGCAATTTATCCATTTAGTTTGTTCTCAGTTAACACTTCATATCCGCTACCGGCGAATCGCACATAGCACAAAAGCCGTTCATTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAAACCGGTTACCCGCATTATATTcc            >  2:709740/1‑250 (MQ=255)
                                                                                                                                                                                                                             agtagtGCGGCAAAAAGCTTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATcgcgcg                                                                             >  2:300850/1‑184 (MQ=255)
                                                                                                                                                                                                                             agtagtGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATcgcgcg                                                                             <  1:300850/184‑1 (MQ=255)
                                                                                                                                                                                                                                agtGCGGCACAAAGCTTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCATGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTATCAGTTCACAATTCATATACGCTACCGGCGAATCGCCCATAGCTCAAACGCCGTTCAGTTTGCGATCGCGCGCCCACTCAGCAAGGGTATACGATTTGATCGACACAGCATGAATGCGGTTCTCCGCAATATATTCCATc         >  2:635269/1‑249 (MQ=255)
                                                                                                                                                                                                                                     ggCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACACTTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCAGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGTGGTATCCGCTTTGATCAAACCAGCATGAATGCGCTTATCCGCAATCTATTCCATCAGCGGAc  >  2:230813/1‑251 (MQ=255)
                                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                                
GCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGAC  >  NC_000913/3336173‑3336652

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