Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,306 C→G 15.7% G73R (GGC→CGC)  aspS ← aspartyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,3060CG15.7% 362.4 / 59.8 141G73R (GGC→CGC) aspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (47/70);  new base G (10/12);  total (59/82)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-02

Only 100 of 144 total aligned reads displayed.

TGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTC  >  NC_000913/1950208‑1950405
                                                                                                   |                                                                                                   
tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGccc                                                                                                    >  2:3668676/1‑100 (MQ=255)
 gcgcGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCg                                                                                                   <  2:38662/100‑1 (MQ=255)
 gcgcGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCg                                                                                                   <  2:38619/100‑1 (MQ=255)
   gcgGTTGATGATAGTCCGCGAGGGCGCCAGAACTTCGATTTCGCCGGTCGCCCTATCGCGGTTAACATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCg                                                                                                   <  1:3018778/98‑1 (MQ=255)
   gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                                 <  1:3747486/100‑1 (MQ=255)
   gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                                 <  1:815167/100‑1 (MQ=255)
   gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                                 >  2:1444977/1‑100 (MQ=255)
    cgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGa                                                                                                <  2:1891247/100‑1 (MQ=38)
    cgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGa                                                                                                <  2:2356171/100‑1 (MQ=255)
    cgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGa                                                                                                <  1:1845594/100‑1 (MQ=255)
    cgGTTGATGATAGACAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGa                                                                                                >  2:2561437/1‑100 (MQ=255)
     ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                               <  1:3563502/100‑1 (MQ=38)
     ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                               <  2:3343654/100‑1 (MQ=38)
     ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                               <  2:3045737/100‑1 (MQ=38)
     ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                                               >  2:553780/1‑100 (MQ=255)
     ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                                               <  1:2815395/100‑1 (MQ=255)
     ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                                               <  2:2193568/100‑1 (MQ=255)
      gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAcc                                                                                              <  2:3538206/100‑1 (MQ=255)
      gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAcc                                                                                              >  1:2182763/1‑100 (MQ=255)
      gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAcc                                                                                              <  2:3683584/100‑1 (MQ=255)
      gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGGCACGGGCACCGTAC‑GTGCGCGTGAcg                                                                                              <  1:1350704/100‑2 (MQ=38)
       ttGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcga                                                                                             <  1:574534/100‑3 (MQ=37)
       ttGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcga                                                                                             <  1:1876422/100‑3 (MQ=37)
       ttGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcga                                                                                             <  2:101669/100‑3 (MQ=37)
        tgatgcTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTg                                                                                            <  1:712976/100‑1 (MQ=255)
          atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGaa                                                                                          >  2:2029209/1‑100 (MQ=255)
          aggatAGTCAGCGGGGACGCCAGCACTTCGAATTCGCCGGTCACCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGaa                                                                                          >  2:2571215/3‑100 (MQ=255)
             atAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGc                                                                                       <  1:2422140/100‑1 (MQ=255)
             atAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGc                                                                                       >  1:1936426/1‑100 (MQ=255)
             atAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGc                                                                                       <  2:212395/100‑1 (MQ=255)
             atAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGc                                                                                       <  1:2970790/100‑1 (MQ=255)
                gTCAGCGAGGACGCCAGCACTTCTATTTCGCCGGTCGCCATATCGCGGTTACTATTGCTTTCGTCACGCGCC‑CGTAAGGTGCTCCTGACTTGatttctga                                                                                    >  1:3059092/1‑93 (MQ=255)
                 tCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCaaga                                                                                   >  1:2186998/1‑97 (MQ=255)
                  cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                                  <  1:1637621/100‑1 (MQ=255)
                  cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                                  <  2:1213758/100‑1 (MQ=255)
                   aGCGAGGCCGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTCAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACt                                                                                 <  2:2274328/100‑1 (MQ=255)
                   aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGa                                                                                    <  1:508121/97‑1 (MQ=255)
                   aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACt                                                                                 >  2:2736556/1‑100 (MQ=255)
                    gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTc                                                                                >  2:2897084/1‑100 (MQ=255)
                     cGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                                               <  2:1165447/100‑1 (MQ=255)
                      gAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCAt                                                                              >  2:3574931/1‑100 (MQ=255)
                       aGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCAtt                                                                             <  1:553774/100‑1 (MQ=255)
                        ggACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                            >  2:3526701/1‑100 (MQ=255)
                         gACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAc                                                                           <  2:2383732/100‑1 (MQ=255)
                          aCGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACg                                                                          >  1:3251039/1‑100 (MQ=255)
                             ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                                       <  2:902932/100‑1 (MQ=255)
                             ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                                       >  1:920042/1‑100 (MQ=255)
                                gCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTc                                                                    <  1:1372410/100‑1 (MQ=255)
                                   cttcgccgtCGCCCGTCGCCATACTCCGGTTACTATTTTTTTCGTCACGCGGC‑CGTACGGTGCCCGTGACCTCAATGCAGAACTCATTACGCAGTCCaca                                                                 <  1:2276549/92‑3 (MQ=255)
                                   cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTTTTTTTTTTTA‑CATACGGTGCCCGTGACCTGAATGCACAACTCATTACGCAGTTCa                                                                   >  1:1731022/1‑98 (MQ=255)
                                   cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                                 >  1:424428/1‑100 (MQ=255)
                                    ttCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                <  1:678617/100‑1 (MQ=255)
                                     tCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAgg                                                               <  2:2667489/100‑1 (MQ=255)
                                     tCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAgg                                                               >  2:1699947/1‑100 (MQ=255)
                                     tCGATTTCACCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAgg                                                               <  2:690196/100‑1 (MQ=255)
                                      cGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGc                                                              >  2:1813498/1‑100 (MQ=255)
                                        aTTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCa                                                            >  2:3356495/1‑100 (MQ=255)
                                             gCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCttt                                                       <  1:2858998/100‑1 (MQ=255)
                                               cGGTCGCCATCTCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGGCCTGCCTGCAGAACTCATTACGCAGTTCGGAGGCCAGTTTTa                                                      <  1:3570033/99‑1 (MQ=255)
                                               cGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTaa                                                     <  2:3344358/100‑1 (MQ=255)
                                                ggTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAc                                                    <  2:1342093/100‑1 (MQ=255)
                                                    gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGt                                                >  2:3324424/1‑100 (MQ=255)
                                                      cATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTcc                                              >  1:1895510/1‑100 (MQ=255)
                                                      cATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCATcc                                              >  1:1896465/1‑100 (MQ=255)
                                                         atCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGTTCTGAATGCATAACTCATTACGCAGTTCGGAGGCCAGCTTTAACGCGTCCGCa                                           >  1:2681821/1‑100 (MQ=255)
                                                          tCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCAc                                          >  1:3720289/1‑100 (MQ=255)
                                                           cgcgGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACGCATTAAGCAGTTCACAGGCCAGCTTTAACGCGTCAGCACg                                         >  1:2881542/1‑100 (MQ=255)
                                                              ggTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATc                                      >  2:2050190/1‑100 (MQ=255)
                                                                ttAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCg                                    <  2:535987/100‑1 (MQ=255)
                                                                  aaTATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCgg                                   >  2:1788662/1‑99 (MQ=255)
                                                                  aaTATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCgg                                   >  1:3019122/1‑99 (MQ=255)
                                                                   atatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGAt                                 >  2:2887612/1‑100 (MQ=255)
                                                                   atatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGAt                                 <  1:3572309/100‑1 (MQ=255)
                                                                    tatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                                >  2:3761150/1‑100 (MQ=255)
                                                                    tatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                                <  2:3576750/100‑1 (MQ=255)
                                                                     atTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                               <  1:974136/100‑1 (MQ=255)
                                                                      tttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                                <  1:3293854/98‑1 (MQ=255)
                                                                      tttgtttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGAt                                 <  1:3534970/97‑1 (MQ=255)
                                                                      ttgttttCGTCACGCGCA‑CGTACGGTGCCCGGGACCTGGTTGCAGAACGCAGTACGCAGCTCAGAGGCCAGCTTTAACGCGTCCGCACGAGCCGGAt                                 <  1:1037275/94‑1 (MQ=255)
                                                                       ttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaa                             >  1:803554/1‑100 (MQ=255)
                                                                        tttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGACTGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa                            <  1:2398428/100‑1 (MQ=255)
                                                                        tttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa                            <  2:2812728/100‑1 (MQ=255)
                                                                           caggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATAGAAAAAt                         >  1:926400/4‑100 (MQ=37)
                                                                            tggTCACGCGCA‑CGTACGGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAGGCGTCCGCACGATCCGGTGCGAAGAAt                         >  1:1642009/3‑99 (MQ=255)
                                                                            tCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATa                        <  2:2341904/100‑1 (MQ=255)
                                                                             ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                       >  2:2520156/2‑100 (MQ=38)
                                                                             ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                       >  1:691459/2‑100 (MQ=38)
                                                                              gTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                      >  2:1557798/1‑100 (MQ=38)
                                                                              gTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                      <  2:941600/100‑1 (MQ=255)
                                                                               tCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                     <  2:1880163/100‑1 (MQ=255)
                                                                                cACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg                    >  1:1029999/1‑100 (MQ=255)
                                                                                cACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg                    <  1:35119/100‑1 (MQ=255)
                                                                                 aCGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGc                   <  2:2266387/100‑1 (MQ=255)
                                                                                 aCGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGc                   >  2:1536967/1‑100 (MQ=255)
                                                                                       a‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATa             <  2:937536/100‑1 (MQ=255)
                                                                                           tACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCt          >  2:2130679/1‑100 (MQ=255)
                                                                                           tACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACATGCACGATAc            >  1:2493611/1‑98 (MQ=255)
                                                                                               gTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTcgc      <  1:1909927/100‑1 (MQ=255)
                                                                                                tGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTcgcg     <  1:1103822/100‑1 (MQ=255)
                                                                                                   ccGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTc  <  2:1154573/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTC  >  NC_000913/1950208‑1950405

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