Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 65,381 T→C T134A (ACT→GCT)  polB ← DNA polymerase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091365,3810TC97.7% 133.9 / ‑4.1 44T134A (ACT→GCT) polBDNA polymerase II
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (15/27);  total (16/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

GTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGC  >  NC_000913/65136‑65626
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
gTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCAt                                                                                                                                                                                                                                                   <  1:171719/251‑1 (MQ=255)
     gtgCAAACCAGGGGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATAGCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGc                                                                                                                                                                                                                                              <  2:229418/249‑1 (MQ=255)
        cAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGTCATTAACGATAGCGCCATTGTGCat                                                                                                                                                                                                                                            >  1:355551/1‑250 (MQ=255)
         aaaCCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCatat                                                                                                                                                                                                                                          <  1:81881/251‑1 (MQ=255)
         aaaCCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCatat                                                                                                                                                                                                                                          <  1:241127/251‑1 (MQ=255)
         aaaCCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGAGGCGGGTTCCGACGGGACTTAAAGATAACGCAATTGTGCatat                                                                                                                                                                                                                                          >  1:64069/1‑251 (MQ=255)
              aGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGGAGCGGTTTCCGACGGGCATTAACGATATCGCCATTGTGCctattcccc                                                                                                                                                                                                                                     >  1:461315/1‑246 (MQ=255)
                              cagtagcTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGGGGTTTCAATATCAATAGAAAACCACTTGAGCGGCGGACGAAAGGCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCCTTGTGCCTATCACCCTCGAccccacccggtg                                                                                                                                                                                                                     >  1:4662/7‑242 (MQ=255)
                              cagtagcTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGGCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGTCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCggtg                                                                                                                                                                                                                     <  2:355551/245‑1 (MQ=255)
                                   gcTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAAGCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCggtgaggtg                                                                                                                                                                                                                <  1:283422/250‑1 (MQ=255)
                                                            cattacaaaatcgcgcgcggtggCGTCGCCATTCTCTGGACCCAGCATATAAACGCTGCGCTGCCCGCAGCCTTCCAGGCCTATGCAGTCCAGCTCACCGTGGCGGGTGGTTGCAATAGCGATAGAAGCCCAGTTGAGCGGCGGACGATAGTCGGGATGCTGTTTCAGACGGGCATTAACGAGAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGtt                                                                                                                                                                                       <  2:132493/230‑1 (MQ=255)
                                                               ttcaaaaTCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGGGGTTTCAATATCAATAGAAACCCACCTTAGCGGCGGACGATAGGCGGGATGCGGTTTTAGACGGGCATTAAAGCAAGCGCCATTTTGCCTATCAACCTCTACCCACACCGGGGGGGGGAATAACCGCTCCCTCAGATAACGTTCTg                                                                                                                                                                                    >  1:83784/5‑251 (MQ=255)
                                                                       cgaGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCgggcg                                                                                                                                                                             <  1:32868/248‑4 (MQ=255)
                                                                            gcggaggCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATc                                                                                                                                                                       <  1:92475/249‑1 (MQ=255)
                                                                                aggCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGcc                                                                                                                                                                   <  1:366417/251‑1 (MQ=255)
                                                                                      cgccttacgttggCCCCACAATATAAACGATGCGCTGCCCGCAGCCTTACAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCTCATCGGCCTCGTAg                                                                                                                                                             <  2:330722/240‑1 (MQ=255)
                                                                                                            ataATCGATGCGCTGCCCGCAGCCGTCCAGGCCGAGGCAGTACAGCTCCCCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGca                                                                                                                                                                           <  2:117940/215‑1 (MQ=255)
                                                                                                               aaCGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCCCACCgggg                                                                                                                                                                                                                     >  1:194902/1‑168 (MQ=255)
                                                                                                               aaCGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCggtg                                                                                                                                                                                                                     <  2:194902/170‑1 (MQ=255)
                                                                                                                      cgcTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTc                                                                                                                             <  2:345117/251‑1 (MQ=255)
                                                                                                                         tGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGGGGACGATAGTCGGGCTGCGGTTTTAGGCGGGCATTAACGGTAGGGCCATTTTGCATATCACCCTCGACCCCCACCGGGGAGGTGATAAACCGATCCCTCAGATAGCGTTCTTGCGAGGGCACATCGGGCTCGCTGACGGGAAAGCCACCCTCAcaccggccgtttttgtt                                                                                                                          >  1:93031/1‑235 (MQ=255)
                                                                                                                                     tCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAAt                                                                                                                <  1:375269/249‑1 (MQ=255)
                                                                                                                                       cAGGCCGATGCAGTACAGCTCACCGTGGCGGGAGGTTTTAATATCCATAGAAACCCAATTGAGCGGCGGACGATACTGGGGATGAGGATTCACACAGGCATTTACGATAGCGCAATTGTGCATATCAACCTCGACCCACCCCGGGGAGGGGATAAAACGCCCCATCAGAAAACGTTCTGGCGGGCGCACATCGGCCCCGAAAACCGGAAAGCCCCACTCACGCCAGCGCTTTTCGTAATTTCTCAAtttgg                                                                                                            >  1:30244/1‑248 (MQ=255)
                                                                                                                                                        gCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTaa                                                                                           <  1:146973/251‑1 (MQ=255)
                                                                                                                                                                tGGCGGGTGGTTTCAATCTCAATAGAACCCCACTTGAACGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCCGCCTCGGAGACGGGAAAGCCACCTTCAAGCAGGCGGTTTTCGTAATTCCTCCATTGGCGACGGGCGCGGCAGTAAAAGccca                                                                                    >  1:397458/1‑248 (MQ=255)
                                                                                                                                                                    gggTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACAcc                                                                              <  1:361300/251‑1 (MQ=255)
                                                                                                                                                                            tcaatatcaataGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAAAGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACCGTAAAAGCCATACACCGGCTGGCg                                                                      >  1:221559/1‑251 (MQ=255)
                                                                                                                                                                                 gtcaataGAAGCCCTCTGGAGCGGCGGACGATAGTCGGGATGAGGTTGCAGCCGGGCATTAGCGATAGCGCAATTGTGCCGAGCACCCTCGACCCCCACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTTTGGCGGGCGCGCATCGGCCTCGTAGGCGGTAGCGCCACCTTCACGCAGGCGGTTTTCGTAATTCAGCAATTGGCGATGGGCGCGACTGTAAAGGCCATACACCGGCTGGCGGTGa                                                                  <  2:27931/246‑1 (MQ=255)
                                                                                                                                                                                    aataGAAACCCTCTGCAGCTGCGGCCGATAGTCGGGATGCGGTTTCAGACGCTCACTAACGATAGCGCCATTGTGCATATCCCCCTCGACCCTCACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTACCGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGCCTGGCGGTGAAAAt                                                              <  2:57941/250‑1 (MQ=255)
                                                                                                                                                                                       aGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCt                                                           >  1:45993/1‑251 (MQ=255)
                                                                                                                                                                                        gAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGGAAATCCtt                                                          >  1:285959/1‑251 (MQ=255)
                                                                                                                                                                                                cTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCgg                                                                                                                                                                     >  2:77632/1‑137 (MQ=255)
                                                                                                                                                                                                cTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCgg                                                                                                                                                                     <  1:77632/137‑1 (MQ=255)
                                                                                                                                                                                                  tGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAg                                                 <  1:423131/250‑1 (MQ=255)
                                                                                                                                                                                                   gAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCTCCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGTCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGc                                                <  1:208596/250‑1 (MQ=255)
                                                                                                                                                                                                          ggACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCCCATCGGCCTCGTAGACGGTCACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGGGGCAGTAAAGGCCATACACCGGCTGGCGGGGAAAATCCTTTaaacaccacggtgtt                                          >  2:5153/1‑236 (MQ=255)
                                                                                                                                                                                                               aTAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGAACATAGGCCTCGTAGACGGTAACGCCACCTTCCCGACGGCGCTTTTCGTAATTACTCAATTAGCGATGGGCGCCACAGGAAAGGCAAAACAACCGCCGGGCGGGAAAAACCCTTTAACCCACCGGGGTCAGGCgg                                    >  2:460820/1‑249 (MQ=255)
                                                                                                                                                                                                                   ccggtaTGCGGTTTCAGAAGGGCATTAACGATAGCGCCTTTGGGCATATCACCCTCGACCCCCCCCGGTGAGGTGATAATCCGCTCCATCATATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCGTCACGCAGGCGCGTTTCGTAATTCATCAATTGGCGATGGTCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGcc                               <  2:359698/246‑1 (MQ=255)
                                                                                                                                                                                                                      gtagtcGGTTTCAGACGGGTATTAACGATAGCGCAATTGTGCAGATCCCCCCCGACCCACACCGGTGAGGTGATAAACTGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAAGTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCACCGGTGTCAGGCGAACGCCTTg                            <  2:4662/246‑1 (MQ=255)
                                                                                                                                                                                                                               ttCAGACGGGCATTAACGATAGCGCCATTGTTCATACCACCCGCGACCCACCCCGGTGAGGTGAGAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCCGGCGCTGTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCGGGCGGTGAAAATCCGTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCt                    <  2:67694/250‑1 (MQ=255)
                                                                                                                                                                                                                               ttCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTg                   >  1:352382/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ggCATTAATGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGt                           <  1:149253/235‑1 (MQ=255)
                                                                                                                                                                                                                                               cGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAgcg   <  2:45993/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gATAGCGCCCTTGTGCTTATCACCCGCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAgcgc  <  2:167010/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
GTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGC  >  NC_000913/65136‑65626

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: