Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 65,594 A→C F63V (TTT→GTT)  polB ← DNA polymerase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091365,5940AC95.8% 126.5 / ‑3.4 48F63V (TTT→GTT) polBDNA polymerase II
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (15/31);  total (16/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

GGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTT  >  NC_000913/65349‑65824
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
gggAGGCGGTTTCAGACGGGCATTATAGATAGCACCATTGTGCATACCACCCCCGACCCCCACCGGTGAGGTGATAAACCGCTCCATCAGTGAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCTCCTGCACGCAGGCGCTTTTCGTAATTCAGCAATGGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCTGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCtt                                                                                                                                                                                                                                           <  2:32553/251‑1 (MQ=255)
   aTGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGGGAAAATCCTTTAACGCCCGCGGTGTCAAGCGAACGCCTTGtt                                                                                                                                                                                                                                        >  1:701041/1‑251 (MQ=255)
            cAGACGGGCAGTAACGATAGCGCCATGTTGCATATCACCCGCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCa                                                                                                                                                                                                                               <  2:123459/251‑1 (MQ=255)
                       tAACGCTAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGCGATACACCGATCCATCAGATCGCGATCTGGCGGGCGCACATCGGACTCGTAGACGGTTACGCCACCTTCAAGCAGGCGCTTTTCGTAATTCCACCATTGGCCAGGGGCGCGAACGGAAAGGCCATACAACGGCCGGCGGGGTAAACTCCTTAAACGCCAGGGGTGGAGGCGAAAGCCttatttcacacgccaaaaagcac                                                                                                                                                                                                                      >  2:668198/1‑228 (MQ=255)
                            aTAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATAAACCCGCTGGCGGTGGAAATCCTTTAAAGCAAAAGGATGCAAGCGAACGCCTTGTTAACCCTGCAAAataatgtgagctgg                                                                                                                                                                                                               >  1:551647/1‑240 (MQ=255)
                                cgccttggcacatcacctccTCGACCCACCCCGGTTAAGTGATAACCCGCTCCATCAGATAGCGTTCTGGCGGGCTCACATCGGCCTCTTAGACGGTCACGCCTCCTTCACGCACGCTCATTTCGTAATTCATCAATTGGCTATGGGGGCGACAGTAAAGGCCATCCACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGGCAGGCGAACTCCTTGTTCACCCTGTAAAATATGCTGTGCGCggtg                                                                                                                                                                                                            <  2:498561/232‑3 (MQ=255)
                                         gCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATcg                                                                                                                                                                                                  <  1:531674/251‑1 (MQ=255)
                                           atatCCCCCGCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCTCCTTCACGCAGGCGCTGTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATcggc                                                                                                                                                                                                <  2:575412/251‑1 (MQ=255)
                                           atagcACCCCCGTCCCACACCGGTGGGGTGATAAACCGTGCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCGGAGCGCGGTGAACCTGATcggc                                                                                                                                                                                                <  2:8225/247‑1 (MQ=255)
                                                         cctaacCGTTGAGGTGATAAACCTCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCTAACGCCGTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACggca                                                                                                                                                                                 <  2:577322/247‑4 (MQ=255)
                                                                 gtgaggtgaTAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGGAAAAAACCCACGGGctct                                                                                                                                                                         >  1:307017/1‑251 (MQ=255)
                                                                   gaggtgaTAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCtt                                                                                                                                                                        <  1:64270/250‑1 (MQ=255)
                                                                      gtgaTAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGAAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACCCCAGCGGTGTCAGGCGAACGACTTGTTCACCCTGCAAAATATGCTTAGGGCGGGGAACCTGACCGGCGGGAATAAAcaacccaggcgtttgtcg                                                                                                                                                                    >  1:403798/1‑234 (MQ=255)
                                                                            aaCCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCtt                                                                                                                                                                                                                                           <  2:108316/175‑1 (MQ=255)
                                                                            aaCCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGGGTCAGGCGAACGCCtt                                                                                                                                                                                                                                           >  1:108316/1‑175 (MQ=255)
                                                                                    cATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCTCCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGTCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCt                                                                                                                                                      <  1:80995/251‑1 (MQ=255)
                                                                                    cATCAGATAGCGTCCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCGTCACGCAGGCGCTTTTCGTACTTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCTGTGCAAGCGTAACCt                                                                                                                                                      <  2:522529/251‑1 (MQ=255)
                                                                                     aTCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTg                                                                                                                                                     <  1:154289/251‑1 (MQ=255)
                                                                                        aGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCaa                                                                                                                                                  <  2:701041/251‑1 (MQ=255)
                                                                                                                 cGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCACGCGTAACCTGCAACGGCCCGTTGTCCGTCGccagccag                                                                                                                         <  1:186361/251‑1 (MQ=255)
                                                                                                                     cTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGGCTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGGTGTCCGTCGCCAGCCAGAAgg                                                                                                                     >  1:174477/1‑251 (MQ=255)
                                                                                                                       cGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAAGCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGga                                                                                                                    <  2:157707/250‑1 (MQ=255)
                                                                                                                              ggTAACGCCACCTTCACGCAGGCGTTTATCGTAATTCATCAAGTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCCCCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTg                                                                                                           <  2:488484/251‑1 (MQ=255)
                                                                                                                                   cGCCCCCTTCACGCAGGCTCTTTTCTCAATTCATCAATTGGCGACGGGAGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCt                                                                                                     <  2:174477/251‑1 (MQ=255)
                                                                                                                                      cACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGc                                                                                                  <  1:139121/251‑1 (MQ=255)
                                                                                                                                                cAGGCGCTCTTCGTATTTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGGGGAAACTCCAACGGCCCGTTGTCCGTCGCCCGCCCGAAGGAGAATTCTGTCCCCTGCGGGGGGTccc                                                                                       >  1:379234/1‑251 (MQ=255)
                                                                                                                                                 aGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCg                                                                                      <  1:166601/251‑1 (MQ=255)
                                                                                                                                                     gctgttCGTAACTCAACAATTGGGTGTGTGCGGGACAGTAAAGGCCCTACACCGGCTGGCGGGTAAAAGCCTTTAACGCCAGCGGTGTCAGGCTAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGAAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTTCGCCGCCAGCCAGAAGGAGACTTCTGTCCTGTGCGTGGTTTCCCGCCAg                                                                                  <  2:379666/247‑1 (MQ=255)
                                                                                                                                                        tttCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAgtg                                                                                <  1:261414/250‑1 (MQ=255)
                                                                                                                                                          tCGTAATTCACCAATTGGCGATGGGCGCGACAGTAAAGGCCATACCCCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCCCCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGGCCCTTGCGGGTTGTCCCGCCAGTGTc                                                                             <  1:194059/250‑1 (MQ=255)
                                                                                                                                                              aaTTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAAAGCCTTGTTCCCCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGCCTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGtt                                                                       <  1:206250/251‑1 (MQ=255)
                                                                                                                                                                ttCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCCCGGACCCTTGCGGTGCCAACGGAACCTGCAACGGCCCCTTGTCCGTCCCCAGCCCGAAGGAGACTTCTGTCCCTTTCGGGGGGGCCCGcccggggcggggta                                                                        >  1:677838/1‑238 (MQ=255)
                                                                                                                                                                ttCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAAAAAACGCCCCGGACTCTTGCGGTGCAAACGTAACCTGCAACGGCCCGGTGTCCGTCGCCCGCCCGAAGGGAGCCTCTGTCCCCTGCGGGGGTGCCCGCCCGTGTCggggtt                                                                      >  1:276447/1‑247 (MQ=255)
                                                                                                                                                                      aaTTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGGATAAACCCCCCGGACTCTTTCGGGGCCAGCGTAACCCGCAACGGCCCCGTGTCCGTCGCCCGCCCCAAGGAGACCTCTGTCCCTTTCGGGGGGTcccccccgagtccgggttagaaga                                                                  >  1:688177/1‑229 (MQ=255)
                                                                                                                                                                              aTGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGGATAAAAGCCCCGGACTCTTGCGGTGCAAGCGTAAACTGCAACCGCCCGTTGTCCGTCGCCCGCCACAAGGAGAATTCTGTCCCTTGCGGGGGGTCCCGCCCGTGGCCGGGTTAGATaaaaacatgcc                                                       >  1:293315/1‑244 (MQ=255)
                                                                                                                                                                                     ggACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCCCCCTGCAAAATATGCTGAGCGCGGGTACCCTGATCGGCGGGAATAAACGCCACGGACTTTTGCGGGGCAAGCGTAACCTGCAACGCCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGCCCCTTGCGGGTTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCAc                                                <  1:264223/250‑1 (MQ=255)
                                                                                                                                                                                         aGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTg                                            <  1:115244/251‑1 (MQ=255)
                                                                                                                                                                                            aagggccctacaccgtcgtgCGGGGAAAACCCTTTACCGCCAGGGGTGTCAGTCGAACGCCTTGTTCACCCTGCAAAATAGGCGGAGCGCGGGGACCCTGATCGGCGGAAATAACCGCCACGGACTCTTGCGGTGCAAGCTTACCCGTCACCGGCCCTGTGTCCTTCGCCAGCCAGAAGGAGTCTTCTTCCCCTTGCGGGGTCCCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCTACGCTgta                                          <  2:701050/231‑3 (MQ=255)
                                                                                                                                                                                              aGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCCGCGGTGTCAGGCGACCGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGACAAAACCACACGGGCTCTTAGGGTGCGAACGTAACACGCAACGGGCCGCTTTCCCTCGCCAGCCCGAAGGAGACTTCTGTACCTTGCGGGGTGTCCCGACAGGGTCGGGGGAAAATAAAAACTGCCCGCGCCAccccgtcatt                                       >  1:709255/1‑242 (MQ=255)
                                                                                                                                                                                                    aacacCGGCTGGCGTGAAAAACCCTTTAACGCCAGCTGTGTCAGGCGATCGCCTTGTTCACCCTGAAAAATAGCCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCGTAAACGTCCCGGTGGCCGTCGCCAGCCATAAGGAGACTGCTGTCCCTTGCGGTGTGTCCCGCCAGTGTCGGGTTAAGATAAAGCCTGCCTGCGCCACGCTGTAATTCCATCa                               <  2:403798/250‑1 (MQ=255)
                                                                                                                                                                                                     cccccgggggggggTGAAAATCCTTAAACGCCAGCGGGTTCAGGCGAACCCCTTGTTCCCCCTGAAAAATATGCTGGGCGGGGGGAACCTGACCGGGGGGAATAACCGCCAGGGACTCTTGCGGTGAAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGCCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGTAATTCCATCaa                              <  2:138616/239‑1 (MQ=255)
                                                                                                                                                                                                     acacCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGTAATTCCATCa                               <  1:533494/250‑1 (MQ=255)
                                                                                                                                                                                                     acacCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGTAATTCCATCa                               <  1:484242/250‑1 (MQ=255)
                                                                                                                                                                                                            cTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACACCCCGGACTCTTTCGGTGCCAACGTAACCTGCAACGGCCCGTTGTCCCTCCCCCGCCAGAAAGAGACTTTTTTCTCCTTCGGGGGGTCCCCCCCGGGTCCGGGTTAAATTAAAACCTCCTGCCCCAccccttaaattccccaaaaaaac                          >  2:680083/1‑228 (MQ=255)
                                                                                                                                                                                                                 ggTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACAGCCCGTTGTCCGTCGCCAGCAAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGGTAAGATAAAAACCGCCTGCGCCCCGCTGTAATTCCCTTCAAAAAACAGGCtt                  >  1:380146/1‑251 (MQ=255)
                                                                                                                                                                                                                                tAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGAAAACGTAACCCGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCCGTGTCGGGGTAAAATAAAACCTGCCTGCGCCACGCGGTAATTTCCACAAAAAACCAGGCTTTATTATAGcccgg     >  1:326304/1‑248 (MQ=255)
                                                                                                                                                                                                                                tAACGCCAGCGGTGTCAGGCGAACGCCTTGTTAACCCTGCAAAATATGCGGAGCGCGGGGAACCTGACCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGTAATTCCATCAAAAAACCAGGCTTGATTATAGCCTg      <  1:483051/250‑1 (MQ=255)
                                                                                                                                                                                                                                      cAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGACCGGCGGGAATAAACGCCACGGACTCTGGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGTAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGttt  <  1:655918/248‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
GGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTT  >  NC_000913/65349‑65824

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