Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 364,170 2 bp→AA coding (2135‑2136/3075 nt) lacZ ← beta‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913364,1700GA97.1% 208.8 / ‑5.9 70A712V (GCC→GTTlacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (29/39);  minor base C (2/0);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913364,1710GA98.6% 208.6 / ‑6.3 70A712V (GCC→GTT) lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/38);  minor base C (0/1);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCGGCTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCC  >  NC_000913/363935‑364406
                                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                            
ccAGCGGTGCACGGGTGAACTGATCGCTCAGCGGCGTCAGCAGTGGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCACCTTCTGACCATGc                                                                                                                                                                                                                                 <  1:456791/250‑1 (MQ=255)
        gCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAAGCCAATTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCTGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCTACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCACGCGGTCGCGtt                                                                                                                                                                                                                       <  2:456875/251‑1 (MQ=255)
         cACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGTCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTc                                                                                                                                                                                                                      <  1:573214/251‑1 (MQ=255)
                  aCTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTGCTGACCATGCGGGCGCGTTCGGTTGCACt                                                                                                                                                                                                             <  1:501768/251‑1 (MQ=255)
                     gATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCGCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTAc                                                                                                                                                                                                           <  1:56677/250‑1 (MQ=255)
                           gcAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACt                                                                                                                                                                                                    <  1:481828/251‑1 (MQ=255)
                             aGCGGCGACTGCAGTTGGTTTTTATCGCCAATCCACATCTGTTAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATGTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACGTCTGACCATGCGTTCTCGTTCGGTTGCACTACGCGTACtgt                                                                                                                                                                                                  <  2:207214/251‑1 (MQ=255)
                              gcgccgTCAGCAGTTGTTGTGTTTCTCCATTCCACAGCTGTGAAAGTAATCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGAGGCAAAAATCCAGTTCGCTGGTGGGCAGATGCGGGATGGCGTGTGACGCGGCGGGGTGCGTCGCACTGAGGTGTTCCGGCTGTCGCCACTGCTGCCAGGCGCTGATGTTCGTAACTGCTGACCATGCGTTCGCGTTCGTTTGCACTACGCTTACtgtg                                                                                                                                                                                                 <  2:424673/247‑1 (MQ=255)
                               cggcgTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGTCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGGGGTCAGATGCGTGATGGCGTGGGACGCTGCGGGGAGCGTCACACTGAGTTTTGCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGGCCATGCGGGCGCGGTCTGTTTCACTACGCGTGCTGTGa                                                                                                                                                                                                <  2:301661/251‑1 (MQ=255)
                                ggcgTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGTTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCg                                                                                                                                                                                                                                <  1:418413/219‑1 (MQ=255)
                                             gttTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCCTCTCGATGCAAAAAACCATTTCGCTGGTGGTTCGATGCGGGGTGGCGTGGGGCGCGGCGGGGGGGGGCACCCTTGGGGTTTTTCCCCGACGCCCCCGCTGCCCGGCGCGGGTGTGCCCCACTTCTGTACATACGGTCGCGTTCGGGTGCCCCACCCCTCCTGTGGGCCCGAGTTGcccc                                                                                                                                                                                  >  1:429841/1‑250 (MQ=255)
                                                    aTCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCCAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGctct                                                                                                                                                                           <  1:109969/251‑1 (MQ=255)
                                                      cGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTcc                                                                                                                                                                         <  1:94974/251‑1 (MQ=255)
                                                             ccACATCTGTGAAAGAAAGCCTGACTTGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGCGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTTCCCGGCGCTCTCCGGCTGCg                                                                                                                                                                  <  1:514449/251‑1 (MQ=255)
                                                                     gtgAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCCTTTTCCTGGTGGTCAGATGAGGGATGGCGTGGGACGCGGCGGGGGGCGGCACACTTGGGGTTTCCCCCCGACGGCCCCTCTGCCCGGCGCCGCATGGGCCAAATTCTGACCATGCGGGTGGGGTTGGGTGCACCACCCGTACTGTGAGGCAGAGTTGGCCGGCGCGCTCCGGcgtcgggaggtt                                                                                                                                                           >  1:82226/1‑239 (MQ=255)
                                                                     gtgAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGt                                                                                                                                                             <  2:502737/248‑1 (MQ=255)
                                                                     gggaaagaaaGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTc                                                                                                                                                           <  2:337221/248‑1 (MQ=255)
                                                                      ggaaagaaaGCCTGGCTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAg                                                                                                                                                         <  1:597727/250‑1 (MQ=255)
                                                                           gaaagCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGt                                                                                                                                                    <  2:388077/251‑1 (MQ=255)
                                                                             aagCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCATGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTtc                                                                                                                                                  <  1:207102/251‑1 (MQ=255)
                                                                                ccTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGtt                                                                                                                                                                                       <  1:111368/211‑1 (MQ=255)
                                                                                         cGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCTGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGttt                                                                                                                                      <  2:478162/251‑1 (MQ=255)
                                                                                           gTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGCTGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGCGGTTTTTCGCCAGACGCCACTGGTGCCAGGCGCTGATGTGCCCCACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGCACTGTGAGCCAGAGTTGGCCGGCGCTCTCCGGCTGCGGTAGTTAAGGCAGTTCAATCAACTGttttc                                                                                                                                    >  2:620487/1‑249 (MQ=255)
                                                                                           gTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCAGTGAGCGTCACACTGAGGTTTTTCGCCAGACGCCACTGCTGTCCCGCGCTGATGTGCCCAACTTTTGACCCTTCGGTCGCGTTCGGTTGCCCTACTCGTACTGTGAGCCAGAGTTGCCCTGCGCTCTCCGGCTGCGGTCGTTCAGGAAGTTCAATCACATGTTTAc                                                                                                                                    >  2:356074/1‑251 (MQ=255)
                                                                                             tAAATTGCCAACGCTTATTCCCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACCCTGGGGTGTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTTCCCAACTTCTTACCATGCGGTCGCGGTCGGTTGCACTACGCGTACTGTGCGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTGGTTCAGGCAGTTCAATCAACTGTTTAccc                                                                                                                                  >  2:131979/1‑250 (MQ=255)
                                                                                              aaaTTGCCAACACTTATTACCCAGCTCGATGCAAAAATCCATTTCGGTGGGGGTCAGATGCGGGATGGCGTGGGACGCGGAGGGGAGCGTCACAATGCGGTTTTCCGCCAGGCGCCACTTCTACCAGGCGCTGAGGTGCCCAACTTCTGACCATGCGGTAGCGTTCGGTTGCGCTACGCCTACTGTTAGCCAGAGGTGCCCGGCGCTCTCCGGCTGCGGTAGTTCCGGCAGTTCAATAAACTGTTTACCtt                                                                                                                                 >  1:385435/1‑251 (MQ=255)
                                                                                               aaTTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTtg                                                                                                                                >  1:523402/1‑251 (MQ=255)
                                                                                                aTTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGGTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTTAATCAACTATTTACCTtgt                                                                                                                               >  1:270397/1‑251 (MQ=255)
                                                                                                      aaCGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCg                                                                                                                         >  1:177678/1‑251 (MQ=255)
                                                                                                        cGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCCGGCGCTGATGTGCCCAAATTCTGACCATGCGGTCGCGGTCGGTTGCAATACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGGTCAATCAACTGTTTACCTTGTGGCGCGAc                                                                                                                       >  1:39474/1‑251 (MQ=255)
                                                                                                             attaCCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATcc                                                                                                                   >  1:99380/1‑250 (MQ=255)
                                                                                                               taCCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCag                                                                                                                 <  1:322134/250‑1 (MQ=255)
                                                                                                                 cccAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAgg                                                                                                              >  1:193090/1‑251 (MQ=255)
                                                                                                                  ccAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGCTGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCCGACGCCCCTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTCCTGTGAGCCCGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCCGTTCAATCAACTGCTTTCCTTGTTGAGCGACCTCCAGAGGc                                                                                                             >  2:426038/1‑251 (MQ=255)
                                                                                                                  ccAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGGCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGc                                                                                                             >  2:157545/1‑251 (MQ=255)
                                                                                                                    aGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCAc                                                                                                           >  1:480426/1‑251 (MQ=255)
                                                                                                                       tCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCCCTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGGCGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTc                                                                                                        >  2:174306/1‑251 (MQ=255)
                                                                                                                       ccGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTc                                                                                                        <  1:321481/250‑1 (MQ=255)
                                                                                                                              aaaaaTCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCtt                                                                                                 >  1:86064/1‑251 (MQ=255)
                                                                                                                                  aTCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGCGGAACGTCACACTGAGGTTTTCCGCCAGACGACACTGATGCCCGGAGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGTCAGAGTTGCCCGGCGCTCTCCCGCTGCGGTAGTTCCGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGATATCAAGAGGCACTTCCCCGCTTGcc                                                                                              >  1:288502/1‑250 (MQ=255)
                                                                                                                                     cATTGCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCg                                                                                          <  1:228387/251‑1 (MQ=255)
                                                                                                                                             tggtggTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGGTCAATCAACTGATTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTAcc                                                                                   >  1:475904/1‑250 (MQ=255)
                                                                                                                                              ggtggtCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCTGCTTACCAt                                                                                 <  2:587249/251‑1 (MQ=255)
                                                                                                                                                                ggCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGt                                                               <  1:295639/251‑1 (MQ=255)
                                                                                                                                                                ggCGTGGGACACTGCGGGGAGCGTCCCATTGAGGTTTTCCGACAGAGGCCACGGCTGCCAGGCGCTGATGGGCCCAACTTCTGACCATGCGGGCGTGTTCGGTTGCCCTACGGGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGTAGCGACATCCAGAGGCAGTTCACCGCTGGCCAGCGGCTTACCATCCAGCGCCACCATCCAGt                                                               <  2:433237/251‑1 (MQ=255)
                                                                                                                                                                 gCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTg                                                              <  1:262019/251‑1 (MQ=255)
                                                                                                                                                                  cGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGc                                                             >  1:170409/1‑251 (MQ=255)
                                                                                                                                                                      ggACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTGCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGa                                                         <  1:284490/251‑1 (MQ=255)
                                                                                                                                                                          ggggcggGGAGCTGCACACTGGTGTTATCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTGTTGACCATGCGGTCGCGTTCGGTTGCACTACGCGGACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGCGCCTTCATCCAGTGCAGGAGCtt                                                     <  2:117679/249‑2 (MQ=255)
                                                                                                                                                                          gcggcgGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTc                                                     <  2:291315/251‑1 (MQ=255)
                                                                                                                                                                             gcggGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGtt                                                  <  2:439392/251‑1 (MQ=255)
                                                                                                                                                                              tgacgcgCGTCACACTGAGGTTTTCAGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGGTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGAGCAGGAGCTCGtt                                                  >  2:380543/7‑250 (MQ=255)
                                                                                                                                                                                        cccacTGAGGTTTTCCGCCAGACTCCTCTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGGTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCGTACCATCCAGCGCCGCCATCCAGTGCATGAGCTCGTTATCGCTATGAc                                       <  2:611024/249‑1 (MQ=255)
                                                                                                                                                                                           acTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCAGCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGa                                    <  1:566321/251‑1 (MQ=255)
                                                                                                                                                                                            cTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGa                                    >  1:432631/1‑250 (MQ=255)
                                                                                                                                                                                               gAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAAc                                  >  1:467248/1‑250 (MQ=255)
                                                                                                                                                                                                aGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAAc                                  >  1:366723/1‑249 (MQ=255)
                                                                                                                                                                                                aGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTGGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAAGCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAg                                <  1:23539/251‑1 (MQ=255)
                                                                                                                                                                                                 ggTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAgg                               >  1:625175/1‑251 (MQ=255)
                                                                                                                                                                                                    tttCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTa                             <  1:585338/250‑1 (MQ=255)
                                                                                                                                                                                                       ccGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCg                         <  1:57282/251‑1 (MQ=255)
                                                                                                                                                                                                        cGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACCGGAATTCGc                        >  2:449739/1‑251 (MQ=255)
                                                                                                                                                                                                         gCCAGACGCCACTGCTGCCAGGCGCTGATGGGCCCAACTTCTGACCACGCGGTCGCGTTCGGTTGCACTACGCGAACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGGGAGACATCCAGAGGCACTTCAACGCTTCCCAGCGGCTTACCATCCAGCGCCACCATCCAGTACAGGAGCTCGCTATCGCTATGACGTAACGGGTAGTCGCt                       >  2:260039/1‑251 (MQ=255)
                                                                                                                                                                                                         gCCAGACGCCACTGCTGCAAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCt                       >  1:163611/1‑251 (MQ=255)
                                                                                                                                                                                                           cAGACGCCACTGCTGCCATGCGCTGATGTGCCCAACTGCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTgg                     <  2:425319/251‑1 (MQ=255)
                                                                                                                                                                                                                      cgcggcCAGGCGCTTATGTGCCCAACTTCTGACCAGGCGGTCGTGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGCTATCGCTATGACGGAACAGGTACTCGCTGGTCACTTCGAt           <  2:497204/247‑1 (MQ=255)
                                                                                                                                                                                                                            cAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTg     >  1:215811/1‑251 (MQ=255)
                                                                                                                                                                                                                              ggCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTACCTTCGATGGTTTGcc   >  2:88204/1‑251 (MQ=255)
                                                                                                                                                                                                                              ggCGCTGATGTGCCCAACTGTTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGcc   <  1:380543/251‑1 (MQ=255)
                                                                                                                                                                                                                               gcgcTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGAATTCGCTGGTCACTTCGATGGTTTGccc  >  1:412896/1‑251 (MQ=255)
                                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                            
CCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCGGCTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCC  >  NC_000913/363935‑364406

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