Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,306 C→G 17.4% 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,3060CG17.4% 266.8 / 54.4 123G73R (GGC→CGC) aspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (48/51);  new base G (9/12);  total (58/65)
Fisher's exact test for biased strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-02

Only 100 of 127 total aligned reads displayed.

TGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTC  >  NC_000913/1950208‑1950405
                                                                                                   |                                                                                                   
tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGccc                                                                                                    >  2:3824012/1‑100 (MQ=255)
tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGccc                                                                                                    >  1:491365/1‑100 (MQ=255)
  cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTTTTTATGTGCA‑AGTACGGTGCACGt                                                                                                  >  2:893231/1‑100 (MQ=255)
  cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGt                                                                                                  <  2:465288/100‑1 (MQ=255)
   gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                                 >  1:426123/1‑100 (MQ=255)
     ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTTTTCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                               <  1:1030199/100‑1 (MQ=37)
     ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                               <  1:2795032/100‑1 (MQ=38)
     ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                               <  2:1219869/100‑1 (MQ=38)
      gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcg                                                                                              <  1:1520990/100‑2 (MQ=38)
      gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcg                                                                                              <  1:1521636/100‑2 (MQ=38)
      gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcg                                                                                              <  1:693172/100‑2 (MQ=38)
      gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAcc                                                                                              <  1:3339480/100‑1 (MQ=255)
      gTTGATGAGAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                               <  2:2721616/99‑1 (MQ=38)
       ttGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGGCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCt                                                                                             <  2:1752485/100‑1 (MQ=255)
        tgatgaTAGACATCGAGGATGCCAGCATTTTGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGTGCACCGTAC‑GTGCGCGTgtcgaa                                                                                            <  2:3571268/100‑6 (MQ=255)
          atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGaa                                                                                          >  1:36905/1‑100 (MQ=255)
           tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAAt                                                                                         >  1:2456966/1‑100 (MQ=255)
            gatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTTGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATg                                                                                        <  1:2442554/100‑1 (MQ=255)
                gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGAGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTTCCCGTGACCTGAATGCAg                                                                                     >  2:3508311/1‑99 (MQ=255)
                  cAGCGTGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                                  <  1:3663895/100‑1 (MQ=255)
                  cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTAGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                                  >  1:319231/1‑100 (MQ=255)
                   aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACt                                                                                 <  1:160845/100‑1 (MQ=255)
                   aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACt                                                                                 >  2:3530363/1‑100 (MQ=255)
                    gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTc                                                                                <  2:67848/100‑1 (MQ=255)
                     cGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                                               >  1:2942037/1‑100 (MQ=255)
                       agtaCGGCAGCACTTTGATGTCGCCTGTCGCCATCTCGCGGTTATTATTTTTTTCGTCACGCGCA‑CGTAAGGTGCCCGTGCCCAGAATGCAGaca                                                                                  <  2:1347822/92‑3 (MQ=255)
                         gACGCCAGTACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                            <  2:440467/99‑1 (MQ=255)
                           cGCCAGCACTTTGATTTTGCCTGTCTTCATTTTGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGc                                                                         <  1:653240/100‑1 (MQ=255)
                           cGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGc                                                                         >  1:3576122/1‑100 (MQ=255)
                            gCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGACCGTAACCTAAATGCAGAACTCATTACGCa                                                                        >  1:3195188/1‑100 (MQ=255)
                             ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                                       >  2:3703436/1‑100 (MQ=255)
                              cAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGt                                                                      >  1:2144655/1‑100 (MQ=255)
                                gCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTTTTTTTTTTCA‑CGTACGGTGCCCGTTACCTGAATGCAGAACTCATTACGCAGTTc                                                                    >  2:1947241/1‑100 (MQ=255)
                                gCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTc                                                                    <  1:2478539/100‑1 (MQ=255)
                                  aCTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCa                                                                   >  1:2822667/1‑99 (MQ=255)
                                   cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTTTTTTTTTTCA‑CGTAAGATGACCGTGACCTGATTGCAGAACTTATTACGCAGTTCag                                                                  >  2:2246711/1‑99 (MQ=255)
                                   cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                                 >  1:2186504/1‑100 (MQ=255)
                                   cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                                 >  1:736490/1‑100 (MQ=255)
                                   cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                                 <  1:3436484/100‑1 (MQ=255)
                                    ttCGATTTCGCCGGTCGCCATATCGCGGTTCATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTTagag                                                                <  2:598529/100‑1 (MQ=255)
                                    ttCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                >  1:482757/1‑100 (MQ=255)
                                     tCGGTTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCAGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                <  2:300940/99‑1 (MQ=255)
                                       gATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAAAACTCATTACGCAGTTCAGAGGc                                                              >  1:126155/1‑99 (MQ=255)
                                        aTTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCTTCACTCTCA‑CTTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGcc                                                             <  1:1263425/99‑1 (MQ=255)
                                         tttCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGATGCCCGTGACCTGAATGCAGAAATCATTACGCAGTTCAGAGGCCAg                                                           >  1:1034429/1‑100 (MQ=255)
                                         tttCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑AGTACGGTGCCCGTGACCTGAATGCAGAACGCATTACGCAGTTCAGAGGCCAg                                                           <  1:27637/100‑1 (MQ=255)
                                         tttCGCCGGTCGCCATATCGCGGATAATATTTTTTTCATGACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAATCATTAGGCAGTTGAGAGGACa                                                            >  2:3793893/1‑99 (MQ=255)
                                           tCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACTCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCCGCt                                                         >  2:580510/1‑100 (MQ=255)
                                            cGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCtt                                                        >  2:2033753/1‑100 (MQ=255)
                                            cGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCtt                                                        >  2:1994681/1‑100 (MQ=255)
                                             gCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCttt                                                       >  1:13455/1‑100 (MQ=255)
                                              ccGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTa                                                      <  1:1270410/100‑1 (MQ=255)
                                                ggTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAc                                                    <  2:1796740/100‑1 (MQ=255)
                                                ggTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGAAGAAATCATTACGCAGTTCAGAGGCCAGCTTTaa                                                     >  1:2184854/1‑99 (MQ=255)
                                                 gTCGCCATATCGCGGTTAATATTTTTTTCGTTTTGCGTA‑CGTACTGTGCTCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                                   <  1:144851/100‑1 (MQ=255)
                                                 gTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                                   <  1:2631018/100‑1 (MQ=255)
                                                 gTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                                   <  1:876219/100‑1 (MQ=255)
                                                 gTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACGCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                                   <  1:964535/100‑1 (MQ=255)
                                                    gCCATATCGCGCTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcgcg                                                 >  2:56496/1‑99 (MQ=255)
                                                      cATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTc                                               <  2:3385700/99‑1 (MQ=255)
                                                          tCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCa                                           <  2:1022597/99‑1 (MQ=255)
                                                             cgGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGAt                                       >  1:1259593/1‑100 (MQ=255)
                                                                 tAATATTTTTTTCGTCACGCGCT‑CGTACGGTGCCCGTGATCTGAATGCAGAACTTATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCgg                                   <  1:3367702/100‑1 (MQ=255)
                                                                 tAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCgg                                   <  1:1169731/100‑1 (MQ=255)
                                                                    tatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACTTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGTGTTCGCACGATCCGGATc                                <  2:1379872/100‑1 (MQ=255)
                                                                    tatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                                <  1:3359961/100‑1 (MQ=255)
                                                                    tatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAAAGGCAAGCTTTAACGCGTCCGAACGATCCGGATc                                >  1:1748274/1‑100 (MQ=255)
                                                                     atTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGTAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                               <  1:366397/100‑1 (MQ=255)
                                                                     atTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                               >  1:2082590/1‑100 (MQ=255)
                                                                        tttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa                            <  1:2321522/100‑1 (MQ=255)
                                                                         ttttCGTCACGCGCA‑CGTACGGTGTCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCACACGATCCGGATCGaaaa                           <  2:3678095/100‑1 (MQ=255)
                                                                           ttCGTCACGCGCA‑CCTACGGTGCCTATTACCTGAATGCAGAACTCTTTACGCAGTTCAGAGGCTATCTTTAACTCGTCCGCATGATCCGAATCGAAAAAt                         <  2:881015/100‑1 (MQ=255)
                                                                            tcttCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATa                        <  1:856131/97‑1 (MQ=38)
                                                                            aggTCACGGGCACCGTAC‑GTGCGCGTCACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATa                        >  1:2384458/3‑100 (MQ=37)
                                                                             ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTTACGCGTCCGCACGAACCGGATCGAAAAATAc                       >  2:1758529/2‑100 (MQ=37)
                                                                             ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                       >  1:3298204/2‑100 (MQ=38)
                                                                             cGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                       >  2:386908/1‑100 (MQ=255)
                                                                              gTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                      >  1:1645130/1‑100 (MQ=38)
                                                                              gTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                      <  1:2336015/100‑1 (MQ=255)
                                                                              gTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCAGATCGAAAAATAcc                      >  1:1938298/1‑100 (MQ=255)
                                                                               tcacggacaccGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                     >  1:709294/11‑100 (MQ=255)
                                                                               tCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATaaga                     >  2:856121/1‑97 (MQ=38)
                                                                               tCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                     >  1:3383599/1‑100 (MQ=255)
                                                                               tCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                     <  1:2140663/100‑1 (MQ=255)
                                                                                cACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg                    >  1:1040248/1‑100 (MQ=255)
                                                                                cACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg                    >  2:1985561/1‑100 (MQ=255)
                                                                                 aCGCGCC‑CGTACGGTGCCCGTGCCCTGAATGCAGAACTCACTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGc                   <  2:3858897/100‑1 (MQ=255)
                                                                                 aCGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg                    <  2:2317905/99‑1 (MQ=255)
                                                                                 aCGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGc                   <  2:2915271/100‑1 (MQ=255)
                                                                                   gtgcA‑CGTACGGTGGCCGTGACATGAATGCAAACTTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATTCGGATCGAAAAATACCTGCAc                 <  2:1624416/98‑1 (MQ=255)
                                                                                   gcgcA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCa                  <  1:3511649/99‑1 (MQ=255)
                                                                                    cgcA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCAc                 >  1:3189267/1‑99 (MQ=255)
                                                                                      caccGTAC‑GTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGAt              >  2:2155971/4‑100 (MQ=255)
                                                                                         cGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATAc            <  2:2456943/100‑1 (MQ=255)
                                                                                             cGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTc        >  1:180759/1‑100 (MQ=255)
                                                                                             cGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCACCTTTAACGCGTCCGCACCATCCGGATCGCAAAATCCCTGCACAATACCCTc        >  2:3168449/1‑100 (MQ=255)
                                                                                               gTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGAACCGGAACGAAAAATACCTGCACGATACCTTcgc      >  1:2042465/1‑100 (MQ=255)
                                                                                                 gCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCgg    <  1:2780262/100‑1 (MQ=255)
                                                                                                  cccGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATACGAATCGAAAAATACCTGCACGATACCTTCGCGGt   <  2:69197/100‑1 (MQ=255)
                                                                                                   ccGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTc  <  2:608897/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 < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑