Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,302 Δ1 bp 17.7% coding (221/1773 nt) aspS ← aspartyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,3010G.17.7% 328.9 / 77.3 124coding (222/1773 nt)aspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (50/50);  new base . (10/12);  total (60/64)
Fisher's exact test for biased strand distribution p-value = 8.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-02

Only 100 of 126 total aligned reads displayed.

AACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTC  >  NC_000913/1950202‑1950399
                                                                                                    |                                                                                                  
aaCATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACg                                                                                                    <  2:2622729/100‑1 (MQ=255)
  cATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCAGTCGCCATATCGCAGTTAATATTTTTTTCGTCACGCGCACGTACg                                                                                                    >  1:432854/1‑98 (MQ=255)
   aTCTGCGCGGTTGATGATAGTTAGCGATGACGTTATTATTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGt                                                                                                  <  1:1096280/99‑1 (MQ=255)
      tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGccc                                                                                              >  2:3824012/1‑100 (MQ=255)
      tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGccc                                                                                              >  1:491365/1‑100 (MQ=255)
      tGCGCGGTTGATGATAGTCAGCGAGGACGCAAGCACCTCGATTTCGCCGATCGCCATATCGCGGTTAATATGTTTGTCGACAGGCGCACGTACGGTg                                                                                                 >  2:3050177/1‑97 (MQ=255)
        cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTTTTTATGTGCAAGTACGGTGCACGt                                                                                            >  2:893231/1‑100 (MQ=255)
        cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGt                                                                                            <  2:465288/100‑1 (MQ=255)
         gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTg                                                                                           >  1:426123/1‑100 (MQ=255)
           ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTTTTCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAc                                                                                         <  1:1030199/100‑1 (MQ=37)
           ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAc                                                                                         <  2:1219869/100‑1 (MQ=38)
           ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAc                                                                                         <  2:1201011/100‑1 (MQ=38)
           ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAc                                                                                         <  1:2795032/100‑1 (MQ=38)
            gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGTCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcg                                                                                        <  1:3200376/100‑2 (MQ=37)
            gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcg                                                                                        <  1:693172/100‑2 (MQ=38)
            gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcg                                                                                        <  1:1520990/100‑2 (MQ=38)
            gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGAcc                                                                                        <  1:3339480/100‑1 (MQ=255)
            gTTGATGAGAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAc                                                                                         <  2:2721616/99‑1 (MQ=38)
             ttGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGGCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCt                                                                                       <  2:1752485/100‑1 (MQ=255)
              tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTg                                                                                      <  1:3573579/100‑1 (MQ=255)
              tgatgaTAGACATCGAGGATGCCAGCATTTTGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGTGCACCGTACGTGCGCGTgtcgaa                                                                                      <  2:3571268/100‑6 (MQ=255)
                atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGTCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcgaaaa                                                                                    <  1:1179590/100‑6 (MQ=25)
                 tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAAt                                                                                   >  1:2456966/1‑100 (MQ=255)
                    tAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCa                                                                                >  2:3570097/1‑100 (MQ=255)
                      gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGAGGTTAATATTTTTTTCGTCACGCGCACGTACGGTTCCCGTGACCTGAATGCAg                                                                               >  2:3508311/1‑99 (MQ=255)
                        cAGCGTGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                            <  1:3663895/100‑1 (MQ=255)
                        cAGCGAGGACGCCAGCACTTCGATTTCTCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                            <  1:636010/100‑1 (MQ=255)
                        cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGTTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                            <  1:1594674/100‑1 (MQ=255)
                         aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACt                                                                           >  2:3530363/1‑100 (MQ=255)
                         aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACt                                                                           <  1:160845/100‑1 (MQ=255)
                          gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTc                                                                          <  2:67848/100‑1 (MQ=255)
                           cGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                                         >  1:2942037/1‑100 (MQ=255)
                             agtaCGGCAGCACTTTGATGTCGCCTGTCGCCATCTCGCGGTTATTATTTTTTTCGTCACGCGCACGTAAGGTGCCCGTGCCCAGAATGCAGaca                                                                            <  2:1347822/92‑3 (MQ=255)
                               gACGCCAGTACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                      <  2:440467/99‑1 (MQ=255)
                                 cGCCAGCACTTTGATTTTGCCTGTCTTCATTTTGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGc                                                                   <  1:653240/100‑1 (MQ=255)
                                 cGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGc                                                                   >  1:3576122/1‑100 (MQ=255)
                                 cGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGc                                                                   >  1:1337159/1‑100 (MQ=255)
                                  gCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGACCGTAACCTAAATGCAGAACTCATTACGCa                                                                  >  1:3195188/1‑100 (MQ=255)
                                   ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                                 >  2:2140644/1‑100 (MQ=255)
                                   ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                                 >  2:3703436/1‑100 (MQ=255)
                                    cAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGt                                                                >  1:2144655/1‑100 (MQ=255)
                                      gCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTTTTTTTTTTCACGTACGGTGCCCGTTACCTGAATGCAGAACTCATTACGCAGTTc                                                              >  2:1947241/1‑100 (MQ=255)
                                      gCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTc                                                              <  1:2478539/100‑1 (MQ=255)
                                        aCTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCa                                                             >  1:2822667/1‑99 (MQ=255)
                                         cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTTTTTTTTTTCACGTAAGATGACCGTGACCTGATTGCAGAACTTATTACGCAGTTCag                                                            >  2:2246711/1‑99 (MQ=255)
                                         cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                           <  1:3436484/100‑1 (MQ=255)
                                         cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                           >  1:2186504/1‑100 (MQ=255)
                                         cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                           >  1:736490/1‑100 (MQ=255)
                                          ttCGATTTCGCCGGTCGCCATATCGCGGTTCATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTTagag                                                          <  2:598529/100‑1 (MQ=255)
                                           tCGGTTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCAGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                          <  2:300940/99‑1 (MQ=255)
                                             gATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAAAACTCATTACGCAGTTCAGAGGc                                                        >  1:126155/1‑99 (MQ=255)
                                              aTTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCTTCACTCTCACTTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGcc                                                       <  1:1263425/99‑1 (MQ=255)
                                               tttCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGATGCCCGTGACCTGAATGCAGAAATCATTACGCAGTTCAGAGGCCAg                                                     >  1:1034429/1‑100 (MQ=255)
                                               tttCGCCGGTCGCCATATCGCGGATAATATTTTTTTCATGACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAATCATTAGGCAGTTGAGAGGACa                                                      >  2:3793893/1‑99 (MQ=255)
                                                 tCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACTCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCCGCt                                                   >  2:580510/1‑100 (MQ=255)
                                                  cGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCtt                                                  >  2:1994681/1‑100 (MQ=255)
                                                   gCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCttt                                                 >  1:13455/1‑100 (MQ=255)
                                                    ccGGTCGCCATATCGCGGTTAATATTTTTTTTTTTATTCGCAAGTACGGAGCCCACGACCTGAATGCAGAAATCTGTACGTAGTTCAGATACCAGCtt                                                  >  2:1919062/1‑98 (MQ=255)
                                                      ggTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTa                                                <  2:314311/98‑1 (MQ=255)
                                                      ggTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGAAGAAATCATTACGCAGTTCAGAGGCCAGCTTTaa                                               >  1:2184854/1‑99 (MQ=255)
                                                       gTCGCCATATCGCGGTTAATATTTTTTTCGTTTTGCGTACGTACTGTGCTCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                             <  1:144851/100‑1 (MQ=255)
                                                       gTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACTTACGGTGTCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                             <  1:2113084/100‑1 (MQ=255)
                                                       gTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                             <  1:876219/100‑1 (MQ=255)
                                                       gTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                             <  1:2631018/100‑1 (MQ=255)
                                                           ccATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTc                                         >  2:3798042/1‑100 (MQ=255)
                                                            cATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTcc                                        >  1:2724517/1‑100 (MQ=255)
                                                                tCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCa                                     <  2:1022597/99‑1 (MQ=255)
                                                                   cgGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGAt                                 >  1:1259593/1‑100 (MQ=255)
                                                                       tAATATTTTTTTCGTCACGCGCTCGTACGGTGCCCGTGATCTGAATGCAGAACTTATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCgg                             <  1:3367702/100‑1 (MQ=255)
                                                                        aaTATTTTTTTTGTCACGCGCACGTACGTTGCCCGTGACATGAATACAGAACTCATTACGCAGTTCAGATGCGAACTTTAACGCGTCGTCGCGATccag                             >  2:58414/1‑97 (MQ=255)
                                                                          tatTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                          <  2:2695714/100‑1 (MQ=255)
                                                                          tatTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAAAGGCAAGCTTTAACGCGTCCGAACGATCCGGATc                          >  1:1748274/1‑100 (MQ=255)
                                                                           atTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGTAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                         <  1:366397/100‑1 (MQ=255)
                                                                           atTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                         >  1:2082590/1‑100 (MQ=255)
                                                                                tttCGTCACGCGCACGTACTGTGCCCGTGACCTGAATGCAGAATTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaaa                    <  1:3698019/100‑1 (MQ=255)
                                                                                ttacGTCACTCGCAAGTTAGGTTATCGTGATCTGAATGCTGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaa                     <  1:2210351/96‑1 (MQ=255)
                                                                                 ttCGTCACGCGCACCTACGGTGCCTATTACCTGAATGCAGAACTCTTTACGCAGTTCAGAGGCTATCTTTAACTCGTCCGCATGATCCGAATCGAAAAAt                   <  2:881015/100‑1 (MQ=255)
                                                                                  tcttCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATa                  <  1:856131/97‑1 (MQ=38)
                                                                                  aggTCACGGGCACCGTACGTGCGCGTCACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATa                  >  1:2384458/3‑100 (MQ=37)
                                                                                   ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTTACGCGTCCGCACGAACCGGATCGAAAAATAc                 >  2:1758529/2‑100 (MQ=37)
                                                                                   ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                 >  1:2398707/2‑100 (MQ=38)
                                                                                   ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                 >  1:3298204/2‑100 (MQ=38)
                                                                                   ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATAGAAAAATAc                 >  1:742235/2‑100 (MQ=37)
                                                                                   cGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                 >  2:386908/1‑100 (MQ=255)
                                                                                    gTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                >  1:3363193/1‑100 (MQ=38)
                                                                                    gTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                >  1:1645130/1‑100 (MQ=38)
                                                                                    gTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                <  1:2336015/100‑1 (MQ=255)
                                                                                    gTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCAGATCGAAAAATAcc                >  1:1938298/1‑100 (MQ=255)
                                                                                     tCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATaaga               >  2:856121/1‑97 (MQ=38)
                                                                                     tCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt               >  1:3383599/1‑100 (MQ=255)
                                                                                      cACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg              >  1:1040248/1‑100 (MQ=255)
                                                                                      cACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg              >  2:1985561/1‑100 (MQ=255)
                                                                                       aCGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg              <  2:2317905/99‑1 (MQ=255)
                                                                                        cgcgCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCa            >  2:3903000/1‑100 (MQ=255)
                                                                                         gtgcACGTACGGTGGCCGTGACATGAATGCAAACTTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATTCGGATCGAAAAATACCTGCAc           <  2:1624416/98‑1 (MQ=255)
                                                                                          cgcACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCAc           >  1:3189267/1‑99 (MQ=255)
                                                                                            caccGTACGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGAt        >  2:2155971/4‑100 (MQ=255)
                                                                                             aCGTACGGTGCCCGTGAGCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATa       <  2:1432012/100‑1 (MQ=255)
                                                                                                   cGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTc  >  1:180759/1‑100 (MQ=255)
                                                                                                   cGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCACCTTTAACGCGTCCGCACCATCCGGATCGCAAAATCCCTGCACAATACCCTc  >  2:3168449/1‑100 (MQ=255)
                                                                                                    |                                                                                                  
AACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTC  >  NC_000913/1950202‑1950399

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