Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,717,170 C→A 5.9% 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.9% 450.7 / 10.6 131L139L (CTC→CTA) ‡unguracil‑DNA‑glycosylase
Reads supporting (aligned to +/- strand):  ref base C (55/68);  new base A (3/5);  total (58/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

TGAATATGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCACCGCAT  >  NC_000913/2717022‑2717314
                                                                                                                                                    |                                                                                                                                                
tGAATATGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTc                                                                                                                                                  <  3:104198/149‑1 (MQ=255)
  aaTATGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCgg                                                                                                                                                <  4:790480/149‑1 (MQ=255)
   atatGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCgg                                                                                                                                                >  2:140716/1‑148 (MQ=255)
   atatGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGc                                                                                                                                               >  4:441747/1‑149 (MQ=255)
    tatGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGc                                                                                                                                               >  2:592373/1‑148 (MQ=255)
     atGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTg                                                                                                                                             <  3:441747/149‑1 (MQ=255)
     atGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCGAGGCTg                                                                                                                                             >  3:108854/1‑149 (MQ=255)
        tataAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTggg                                                                                                                                           <  1:140716/148‑1 (MQ=255)
         ataAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGa                                                                                                                                          >  2:89781/1‑148 (MQ=255)
          taAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGaaa                                                                                                                                        >  3:726040/1‑149 (MQ=255)
             agagCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAAcc                                                                                                                                      >  2:610299/1‑148 (MQ=255)
             agagCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAAcc                                                                                                                                      >  1:42114/1‑148 (MQ=255)
                    gAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGa                                                                                                                              >  3:835052/1‑149 (MQ=255)
                     aaaaTACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGa                                                                                                                              >  1:489527/1‑148 (MQ=255)
                        aTACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATaa                                                                                                                           <  2:911931/148‑1 (MQ=255)
                         tACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAg                                                                                                                         >  3:909991/1‑149 (MQ=255)
                            tATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAg                                                                                                                         >  4:604544/1‑146 (MQ=255)
                            tATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAg                                                                                                                         <  3:604544/146‑1 (MQ=255)
                            tATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGa                                                                                                                      >  4:402683/1‑149 (MQ=255)
                            tATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGa                                                                                                                      >  3:281402/1‑149 (MQ=255)
                             aTTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCGAGGCTggcgtg                                                                                                                                        <  1:816608/130‑5 (MQ=255)
                             aTTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCGAGGCTggcgtg                                                                                                                                        >  2:816608/1‑126 (MQ=255)
                              ttCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATc                                                                                                                    >  4:756361/1‑149 (MQ=255)
                                 cGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAg                                                                                                                  >  2:250857/1‑148 (MQ=255)
                                 cGGACTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAg                                                                                                                  <  1:349081/148‑1 (MQ=255)
                                  gggCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGcc                                                                                                                <  4:425342/149‑1 (MQ=255)
                                      ttCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGa                                                                                                             <  1:167390/148‑1 (MQ=255)
                                       tCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGAtt                                                                                                           <  4:899278/149‑1 (MQ=255)
                                       tCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGAt                                                                                                            >  2:965397/1‑148 (MQ=255)
                                         aCCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTaa                                                                                                         <  4:736379/149‑1 (MQ=255)
                                            cGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCa                                                                                                      <  3:612425/149‑1 (MQ=255)
                                             gCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTc                                                                                                                                   <  1:611643/119‑1 (MQ=255)
                                             gCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTc                                                                                                                                   >  2:611643/1‑119 (MQ=255)
                                              cccTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGc                                                                                                    >  4:373532/1‑149 (MQ=255)
                                               ccTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTATGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATaaa                                                                                                                          >  4:738887/1‑126 (MQ=255)
                                                cTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCa                                                                                                   >  1:127227/1‑148 (MQ=255)
                                                  aaTCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATc                                                                                                 <  2:5319/148‑1 (MQ=255)
                                                     cATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATcgcg                                                                                              >  2:348544/1‑148 (MQ=255)
                                                      aTGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGaa                                                                                            <  4:338251/149‑1 (MQ=255)
                                                            aTCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCgtg                                                                                      <  4:491449/149‑1 (MQ=255)
                                                            aTCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCgt                                                                                       <  2:489527/148‑1 (MQ=255)
                                                               ttGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATCCTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCgtggtg                                                                                   >  4:697913/1‑149 (MQ=255)
                                                               ttGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCGAGGCTggcgtg                                                                                                                                        <  3:775387/96‑5 (MQ=255)
                                                               ttGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCGAGGCTggcatg                                                                                                                                        >  4:775387/1‑92 (MQ=255)
                                                                tGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTg                                                                                                                                             <  2:380630/90‑1 (MQ=255)
                                                                tGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTg                                                                                                                                             >  1:380630/1‑90 (MQ=255)
                                                                tGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCgtggtg                                                                                   >  1:890673/1‑148 (MQ=255)
                                                                 gAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAg                                                                                                     <  4:232478/129‑1 (MQ=255)
                                                                 gAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGgtgt                                                                                  <  1:898771/148‑1 (MQ=255)
                                                                 gAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGtt                                                                                 <  3:620028/149‑1 (MQ=255)
                                                                  aaaGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGtt                                                                                 >  2:706738/1‑148 (MQ=255)
                                                                    aGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGtttt                                                                               <  1:111493/148‑1 (MQ=255)
                                                                    aGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGtttt                                                                               <  2:912156/148‑1 (MQ=255)
                                                                    aGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGtttt                                                                               <  2:319299/148‑1 (MQ=255)
                                                                        gggCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTttgtt                                                                           <  1:592373/148‑1 (MQ=255)
                                                                           cgcgTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTtgtg                                                                        <  2:890673/148‑1 (MQ=255)
                                                                           cgcgTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTgg                                                                       <  4:588148/149‑1 (MQ=255)
                                                                                  gggCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGAt                                                                   <  2:274280/146‑1 (MQ=255)
                                                                                   ggCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGc                                                                >  2:930677/1‑148 (MQ=255)
                                                                                      gTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATgc                                                            >  3:457602/1‑149 (MQ=255)
                                                                                                tCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAaaagaaa                                                  <  4:707515/149‑1 (MQ=255)
                                                                                                tCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAaaagaa                                                   <  1:222536/148‑1 (MQ=255)
                                                                                                tCAATACGGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCGAGTCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAaaagaa                                                   >  2:287685/1‑148 (MQ=255)
                                                                                                 cAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAaaagaaag                                                 <  4:281402/149‑1 (MQ=255)
                                                                                                  aaTACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAAc                                                                                                                                       <  4:489246/62‑1 (MQ=255)
                                                                                                  aaTACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAAc                                                                                                                                       >  3:489246/1‑62 (MQ=255)
                                                                                                  aaTACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAgg                                                >  4:25734/1‑149 (MQ=255)
                                                                                                    tACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAggg                                               <  1:851044/148‑1 (MQ=255)
                                                                                                     aCTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGc                                             <  4:815903/149‑1 (MQ=255)
                                                                                                       tgtgTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCg                                            >  2:909271/1‑148 (MQ=255)
                                                                                                          gtTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATta                                        <  3:290553/149‑1 (MQ=255)
                                                                                                            tGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATtat                                       <  1:352361/148‑1 (MQ=255)
                                                                                                              aCGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAg                                     <  1:138770/148‑1 (MQ=255)
                                                                                                              aCGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGATGGGAAACCTTCACCGATAAAGTGATCAGCCTAATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGa                                    <  4:681884/149‑1 (MQ=255)
                                                                                                               cGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGa                                    <  1:801668/148‑1 (MQ=255)
                                                                                                                 gTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTgg                                                                       >  3:71169/1‑111 (MQ=255)
                                                                                                                 gTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATaa                                 <  3:49445/149‑1 (MQ=255)
                                                                                                                 gTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATa                                  >  1:196081/1‑148 (MQ=255)
                                                                                                                 gTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGACTATAGATa                                  <  2:724681/148‑1 (MQ=255)
                                                                                                                   aCGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAg                                <  2:42114/148‑1 (MQ=255)
                                                                                                                   aCGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGc                               <  3:383256/149‑1 (MQ=255)
                                                                                                                    cgcgCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGc                               <  1:909271/148‑1 (MQ=255)
                                                                                                                      cgcAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCaa                             <  2:895616/148‑1 (MQ=255)
                                                                                                                              aGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCcatcat                    <  3:939494/149‑1 (MQ=255)
                                                                                                                                gcgcATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGt                  <  4:726040/149‑1 (MQ=255)
                                                                                                                                  gcATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATGAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTa                 <  1:930198/148‑1 (MQ=255)
                                                                                                                                  gcATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTa                 <  2:239103/148‑1 (MQ=255)
                                                                                                                                   cATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATtat                                       <  1:52204/125‑1 (MQ=255)
                                                                                                                                   cATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATtat                                       >  2:52204/1‑125 (MQ=255)
                                                                                                                                     ttCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATg                                                             <  1:560636/101‑1 (MQ=255)
                                                                                                                                     ttCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCaaaa                                                      <  3:343534/108‑1 (MQ=255)
                                                                                                                                      tCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGaa            <  4:835052/149‑1 (MQ=255)
                                                                                                                                       cccACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATTGATAAGCAACGCCATCATGTACTGaaa           <  3:639717/149‑1 (MQ=255)
                                                                                                                                         cACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAg          >  1:893291/1‑148 (MQ=255)
                                                                                                                                         cACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGc         >  3:254318/1‑149 (MQ=255)
                                                                                                                                          aCGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCa        >  3:746383/1‑149 (MQ=255)
                                                                                                                                           tcccAGCCGAGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTg              <  1:748111/140‑1 (MQ=255)
                                                                                                                                           cGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCAc       <  4:537823/149‑1 (MQ=255)
                                                                                                                                           cGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCAc       <  3:753457/149‑1 (MQ=255)
                                                                                                                                               aGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCACCGc    <  2:127227/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TGAATATGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCACCGCAT  >  NC_000913/2717022‑2717314

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