Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,859,142 G→A 66.9% A684T (GCC→ACC)  mutS → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,859,1420GA66.9% 345.2 / 278.1 283A684T (GCC→ACC) mutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (47/46);  new base A (90/98);  total (137/146)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

Only 100 of 286 total aligned reads displayed.

CTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACC  >  NC_000913/2858994‑2859291
                                                                                                                                                    |                                                                                                                                                     
cTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGc                                                                                                                                                      >  1:1888639/1‑150 (MQ=255)
     tCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCg                                                                                                                                                 >  1:2584988/1‑150 (MQ=255)
        gCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAAt                                                                                                                                              >  1:3879020/1‑150 (MQ=255)
        gCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAAt                                                                                                                                              >  1:3121978/1‑150 (MQ=255)
        gCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAAt                                                                                                                                              >  2:277492/1‑150 (MQ=255)
          aGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATAc                                                                                                                                            <  1:650798/150‑1 (MQ=255)
            cTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAg                                                                                                                                          >  1:1273560/1‑150 (MQ=255)
            cTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAg                                                                                                                                          >  2:42473/1‑150 (MQ=255)
             tATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGt                                                                                                                                         >  2:1560946/1‑150 (MQ=255)
                gTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTg                                                                                                                                      <  2:926371/150‑1 (MQ=255)
                gTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTg                                                                                                                                      <  1:3769393/150‑1 (MQ=255)
                    cGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGgtgt                                                                                                                                  >  2:3002402/1‑150 (MQ=255)
                      gCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTa                                                                                                                                >  1:336647/1‑150 (MQ=255)
                          aaaaaGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATgg                                                                                                                            <  2:856816/150‑1 (MQ=255)
                           aaaaGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGa                                                                                                                           <  1:3217169/150‑1 (MQ=255)
                              aGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATga                                                                                                                        <  1:1636370/150‑1 (MQ=255)
                              aGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATga                                                                                                                        <  1:169620/150‑1 (MQ=255)
                                  gagaTTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATc                                                                                                                    >  2:4089675/1‑150 (MQ=255)
                                  gagaTTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATc                                                                                                                    >  2:1775769/1‑150 (MQ=255)
                                  gagaTTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATc                                                                                                                    >  2:1623510/1‑150 (MQ=255)
                                  gagaTTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATc                                                                                                                    >  2:2913691/1‑150 (MQ=255)
                                  gagaTTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATc                                                                                                                    >  1:216146/1‑150 (MQ=255)
                                  gagaTTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATc                                                                                                                    <  2:3881003/150‑1 (MQ=255)
                                    gaTTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCgg                                                                                                                  >  2:108628/1‑150 (MQ=255)
                                     aTTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCggg                                                                                                                 <  2:2224280/150‑1 (MQ=255)
                                            cTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGac                                                                                                          >  1:308377/1‑149 (MQ=255)
                                            cTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGaa                                                                                                          >  2:3389657/1‑150 (MQ=255)
                                             tATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAAc                                                                                                         <  1:2377054/150‑1 (MQ=255)
                                             tATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAAc                                                                                                         >  2:1734932/1‑150 (MQ=255)
                                               tcgatcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGt                                                                                                       >  1:2327779/1‑150 (MQ=255)
                                               tcgatcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGt                                                                                                       <  1:1346471/150‑1 (MQ=255)
                                                 gatcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTcc                                                                                                     >  1:3858454/1‑150 (MQ=255)
                                                  atcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCa                                                                                                    >  1:1993256/1‑150 (MQ=255)
                                                  atcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCa                                                                                                    >  1:4107588/1‑150 (MQ=255)
                                                  atcgCATCTTGACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCa                                                                                                    <  2:1925365/150‑1 (MQ=255)
                                                   tcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCAc                                                                                                   <  1:993610/150‑1 (MQ=255)
                                                   tcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCAc                                                                                                   <  1:2952908/150‑1 (MQ=255)
                                                       aTCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGTTATAGGGCGTGGAACTTCCACCTAc                                                                                               >  2:2409149/1‑150 (MQ=255)
                                                       aTCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTAc                                                                                               <  1:3640119/150‑1 (MQ=255)
                                                             aCCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGt                                                                                         >  1:2191298/1‑150 (MQ=255)
                                                               ccGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCt                                                                                       >  2:1114626/1‑150 (MQ=255)
                                                               ccGCGTAGGCGCGGCAGATGACCTGGCGTCCGGCCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCt                                                                                       >  2:3240421/1‑150 (MQ=255)
                                                                 gcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGt                                                                                     >  2:2116931/1‑150 (MQ=255)
                                                                 gcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGt                                                                                     >  2:1590709/1‑150 (MQ=255)
                                                                  cgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTc                                                                                    <  2:1177189/150‑1 (MQ=255)
                                                                  cgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTc                                                                                    <  2:267118/150‑1 (MQ=255)
                                                                  cgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTc                                                                                    <  1:579967/150‑1 (MQ=255)
                                                                   gTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCg                                                                                   >  1:102172/1‑150 (MQ=255)
                                                                       gcgcGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTgg                                                                               >  1:444541/1‑150 (MQ=255)
                                                                        cgcgGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggc                                                                              <  1:3553417/150‑1 (MQ=255)
                                                                        cgcgGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACAGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggc                                                                              >  1:3293399/1‑150 (MQ=255)
                                                                         gcgGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcg                                                                             >  2:3468370/1‑150 (MQ=255)
                                                                          cgGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgt                                                                            <  1:70145/150‑1 (MQ=255)
                                                                            gCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGATATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgg                                                                          <  2:1890511/150‑1 (MQ=255)
                                                                             cAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtggg                                                                         <  2:377811/150‑1 (MQ=255)
                                                                              aGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgggc                                                                        >  2:552520/1‑150 (MQ=255)
                                                                                aTGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgggcgt                                                                      <  1:3389657/150‑1 (MQ=255)
                                                                                  gtccTGGCGTCCGGGCGCTCCACCTTGATGGTGGTGTTGACTTAAACCGCCAAAATTTTACCTTACACCAACGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGCACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTgc                                                                    >  1:1007435/3‑150 (MQ=255)
                                                                                     cTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGcgcg                                                                 <  1:485173/150‑1 (MQ=255)
                                                                                          gTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGaaaa                                                            <  1:856074/150‑1 (MQ=255)
                                                                                          gTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGaaaa                                                            <  2:1789159/150‑1 (MQ=255)
                                                                                               ggCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTgg                                                       <  2:3166275/150‑1 (MQ=255)
                                                                                               ggCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTgg                                                       <  2:934021/150‑1 (MQ=255)
                                                                                               ggCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTgg                                                       >  1:2736430/1‑150 (MQ=255)
                                                                                                   cTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGaa                                                   <  2:2549060/150‑1 (MQ=255)
                                                                                                     cAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATa                                                 <  2:3554875/150‑1 (MQ=255)
                                                                                                      aaCCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATaa                                                >  2:709237/1‑150 (MQ=255)
                                                                                                            tATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTaa                                          <  2:2928401/150‑1 (MQ=255)
                                                                                                            tATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTaa                                          <  2:125236/150‑1 (MQ=255)
                                                                                                             aTGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAg                                         <  2:321602/150‑1 (MQ=255)
                                                                                                             aTGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAg                                         <  2:3145990/150‑1 (MQ=255)
                                                                                                                gtggAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCa                                      >  1:3949398/1‑150 (MQ=255)
                                                                                                                gtggAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCa                                      >  1:1136737/1‑150 (MQ=255)
                                                                                                                  ggAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCAtt                                    <  2:3684672/150‑1 (MQ=255)
                                                                                                                  ggAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCAtt                                    <  1:526824/150‑1 (MQ=255)
                                                                                                                     gaTGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGAc                                 <  1:1114527/150‑1 (MQ=255)
                                                                                                                      aTGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACg                                >  1:2068914/1‑150 (MQ=255)
                                                                                                                        gACTGAAACCGCCAATATTTTACATAACCCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGtt                              <  2:644590/150‑1 (MQ=255)
                                                                                                                        gACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGtt                              <  2:2349450/150‑1 (MQ=255)
                                                                                                                           tGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTAtt                           >  1:3014572/1‑150 (MQ=255)
                                                                                                                           tGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTAtt                           >  1:4210052/1‑150 (MQ=255)
                                                                                                                              aaCCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGc                        <  1:950462/150‑1 (MQ=255)
                                                                                                                              aaCCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGc                        >  1:2969765/1‑150 (MQ=255)
                                                                                                                              aaCCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGc                        <  2:3722913/150‑1 (MQ=255)
                                                                                                                               aCCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCt                       >  1:383247/1‑150 (MQ=255)
                                                                                                                               aCCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCt                       >  1:345035/1‑150 (MQ=255)
                                                                                                                                ccGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTa                      <  1:1814967/150‑1 (MQ=255)
                                                                                                                                ccGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTa                      <  2:3619080/150‑1 (MQ=255)
                                                                                                                                 cGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCTGAAAATCTGGCGAATTAGATTAAGGCATTGACGTTATTTGCTAc                     <  1:3240334/150‑1 (MQ=255)
                                                                                                                                  gCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTAcc                    <  1:157090/150‑1 (MQ=255)
                                                                                                                                   ccAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTAccc                   >  2:1476591/1‑150 (MQ=255)
                                                                                                                                    cAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTTCCCa                  <  1:3434115/150‑1 (MQ=255)
                                                                                                                                    cAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCa                  >  2:1314346/1‑150 (MQ=255)
                                                                                                                                    cAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCa                  <  1:2828033/150‑1 (MQ=255)
                                                                                                                                          tttACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTAttt            >  2:2868862/1‑150 (MQ=255)
                                                                                                                                              cATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAg        <  2:184775/150‑1 (MQ=255)
                                                                                                                                              cATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAg        <  1:1413675/150‑1 (MQ=255)
                                                                                                                                               aTAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGc       >  2:44082/1‑150 (MQ=255)
                                                                                                                                               aTAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGAGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGc       >  2:2492174/1‑150 (MQ=255)
                                                                                                                                                    gCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGAcc  <  2:24611/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACC  >  NC_000913/2858994‑2859291

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