Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,717,170 C→A 5.6% L139L (CTC→CTA) ‡ ung → uracil‑DNA‑glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,717,1700CA 5.6% 255.9 / 5.3 71L139L (CTC→CTA) ‡unguracil‑DNA‑glycosylase
Reads supporting (aligned to +/- strand):  ref base C (31/36);  new base A (2/2);  total (33/38)
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.77e-01

AATATGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCACCGCA  >  NC_000913/2717024‑2717313
                                                                                                                                                  |                                                                                                                                               
aaTATGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCgg                                                                                                                                               <  2:148052/149‑1 (MQ=255)
 atatGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGc                                                                                                                                              >  4:468982/1‑149 (MQ=255)
  tatGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCt                                                                                                                                             >  3:332117/1‑149 (MQ=255)
   atGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTg                                                                                                                                            <  2:451030/149‑1 (MQ=255)
   atGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTg                                                                                                                                            <  4:169509/149‑1 (MQ=255)
      tataAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGa                                                                                                                                         >  2:300016/1‑149 (MQ=255)
        taAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGaaa                                                                                                                                       >  4:463776/1‑149 (MQ=255)
                  gAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGa                                                                                                                             <  1:596681/149‑1 (MQ=255)
                      aTACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCGAGGCTggcgtg                                                                                                                                       <  2:58173/135‑5 (MQ=255)
                      aTACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCGAGGCTggcgtg                                                                                                                                       >  1:58173/1‑131 (MQ=255)
                        aCTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGt                                                                                                                       >  4:579740/1‑149 (MQ=255)
                            ttCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATc                                                                                                                   <  3:23427/149‑1 (MQ=255)
                            ttCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATc                                                                                                                   <  3:541429/149‑1 (MQ=255)
                              ccGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAg                                                                                                                 <  3:468982/149‑1 (MQ=255)
                               cGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGc                                                                                                                >  3:497638/1‑149 (MQ=255)
                                  gCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTg                                                                                                             <  4:95481/149‑1 (MQ=255)
                                   cTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCGAGGCTggc                                                                                                                                          <  3:89960/119‑2 (MQ=255)
                                   cTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCGAGGCTggc                                                                                                                                          >  4:89960/1‑118 (MQ=255)
                                    ttCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGAt                                                                                                           >  3:558294/1‑149 (MQ=255)
                                       aCCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGAt                                                                                                                    <  1:575493/137‑1 (MQ=255)
                                       aCCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGAt                                                                                                                    >  2:575493/1‑137 (MQ=255)
                                        cccGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAAc                                                                                                       >  3:512783/1‑149 (MQ=255)
                                        cccGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAAc                                                                                                       <  1:455066/149‑1 (MQ=255)
                                        cccGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAAc                                                                                                       <  2:21988/149‑1 (MQ=255)
                                        cccGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAAc                                                                                                       <  2:512929/149‑1 (MQ=255)
                                               tAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATaaa                                                                                                                         <  3:194777/124‑1 (MQ=255)
                                               tAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATaaa                                                                                                                         >  4:194777/1‑124 (MQ=255)
                                                 aTCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATcgc                                                                                              >  1:510113/1‑149 (MQ=255)
                                                  tCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATcgcg                                                                                             >  2:41886/1‑149 (MQ=255)
                                                    aTGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCAtgtgatc                                                                                           <  3:463776/149‑7 (MQ=255)
                                                          aTCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCgtg                                                                                     <  4:332117/149‑1 (MQ=255)
                                                           tCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCgtgg                                                                                    <  2:303695/149‑1 (MQ=255)
                                                           tCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCgtgg                                                                                    <  4:207394/149‑1 (MQ=255)
                                                           tCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCgtgg                                                                                    <  4:381588/149‑1 (MQ=255)
                                                             ttGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCgtggtg                                                                                  >  3:513608/1‑149 (MQ=255)
                                                             ttGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCgtggtg                                                                                  >  2:423607/1‑149 (MQ=255)
                                                               gAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGtt                                                                                <  3:579740/149‑1 (MQ=255)
                                                                 aaGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGtttt                                                                              <  4:512783/149‑1 (MQ=255)
                                                                  aGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGttttt                                                                             <  3:253114/149‑1 (MQ=255)
                                                                      gggCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGtttt                                                                              <  1:169972/144‑1 (MQ=255)
                                                                      gggCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGtttt                                                                              >  2:169972/1‑144 (MQ=255)
                                                                        gcgcGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTtgtg                                                                       <  4:558294/149‑1 (MQ=255)
                                                                         cgcgTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTgg                                                                      >  4:239983/1‑149 (MQ=255)
                                                                             tCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTgg                                                                      >  2:155212/1‑145 (MQ=255)
                                                                             tCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTgg                                                                      <  1:155212/145‑1 (MQ=255)
                                                                                  gCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCAt                                                             <  4:88031/149‑1 (MQ=255)
                                                                                    gTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGaaa                                                                                                                                       <  3:420272/73‑1 (MQ=255)
                                                                                    gTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGaaa                                                                                                                                       >  4:420272/1‑73 (MQ=255)
                                                                                     ttCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATgcg                                                          <  2:421368/149‑1 (MQ=255)
                                                                                          cTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCaaaa                                                     <  3:555941/149‑1 (MQ=255)
                                                                                             cTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAaaagaa                                                  >  2:324318/1‑149 (MQ=255)
                                                                                                aaTACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAgg                                               >  1:166084/1‑149 (MQ=255)
                                                                                                aaTACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCACAAAAGAAAgg                                               >  1:494637/1‑149 (MQ=255)
                                                                                                aaTACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCCGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGCTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAgg                                               >  2:90791/1‑149 (MQ=255)
                                                                                                      gtgtTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGAt                                         >  2:572022/1‑149 (MQ=255)
                                                                                                       tgtTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGAtt                                        >  1:432024/1‑149 (MQ=255)
                                                                                                        gtTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATta                                       <  1:423607/149‑1 (MQ=255)
                                                                                                           gACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCa                                                                                                     >  4:403191/1‑84 (MQ=255)
                                                                                                           gACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCa                                                                                                     <  3:403191/84‑1 (MQ=255)
                                                                                                             cGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGAt                                  <  1:264446/149‑1 (MQ=255)
                                                                                                               gTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATaa                                >  1:334440/1‑149 (MQ=255)
                                                                                                                tACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAg                               >  4:8625/1‑149 (MQ=255)
                                                                                                                       aGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGcc                        <  2:334440/149‑1 (MQ=255)
                                                                                                                                gcATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTAc               <  1:300016/149‑1 (MQ=255)
                                                                                                                                   ttCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGa            <  4:497638/149‑1 (MQ=255)
                                                                                                                                      ccACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAg         >  1:157794/1‑149 (MQ=255)
                                                                                                                                       cACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGc        >  1:498888/1‑149 (MQ=255)
                                                                                                                                          gCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCAcc     <  1:41886/149‑1 (MQ=255)
                                                                                                                                           ccAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAACCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTtgttgt                                                                        >  4:345337/1‑81 (MQ=255)
                                                                                                                                           ccAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAACCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTtgttgt                                                                        <  3:345337/81‑1 (MQ=255)
                                                                                                                                             aGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCACCGCa  >  1:159298/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
AATATGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCACCGCA  >  NC_000913/2717024‑2717313

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