Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,296:1 +C 16.1% coding (227/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,2951.C16.1% 418.2 / 82.1 124coding (228/1773 nt)aspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base . (52/52);  new base C (8/12);  total (60/64)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.94e-04

Only 100 of 128 total aligned reads displayed.

GGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATA  >  NC_000913/1950197‑1950394
                                                                                                   |                                                                                                   
  cAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGt                                                                                                  >  1:3766023/1‑100 (MQ=255)
   aGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTa                                                                                                 <  2:1357182/100‑1 (MQ=255)
   aGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGAGGTTAATATTTTTTTCGTCACGCGCA‑CGTa                                                                                                 >  1:154184/1‑100 (MQ=255)
   aGAACATCTGCGCGGTTCATGAGAGTCAGCAAGGACGCAAGCACTTCGCTTTCGGCGGTCGCCATCTCGCGGTCAATATTTTTTTCGACCGGCGCatat                                                                                                   >  2:1280732/1‑96 (MQ=255)
     aaCATCTGCGCGGTTGCTGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCACGGTTAATATTTTTTTCGTCACGAGCA‑CGTAc                                                                                                >  1:236713/1‑99 (MQ=255)
     aaCATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACg                                                                                               <  2:2622729/100‑1 (MQ=255)
       cATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCAGTCGCCATATCGCAGTTAATATTTTTTTCGTCACGCGCA‑CGTACg                                                                                               >  1:432854/1‑98 (MQ=255)
        aTCTGCGCGGTTGATGATAGTTAGCGATGACGTTATTATTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGt                                                                                             <  1:1096280/99‑1 (MQ=255)
           tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGccc                                                                                         >  2:3824012/1‑100 (MQ=255)
           tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGccc                                                                                         >  1:491365/1‑100 (MQ=255)
           tGCGCGGTTGATGATAGTCAGCGAGGACGCAAGCACCTCGATTTCGCCGATCGCCATATCGCGGTTAATATGTTTGTCGACAGGCGCA‑CGTACGGTg                                                                                            >  2:3050177/1‑97 (MQ=255)
             cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTTTTTATGTGCA‑AGTACGGTGCACGt                                                                                       >  2:893231/1‑100 (MQ=255)
             cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGt                                                                                       <  2:465288/100‑1 (MQ=255)
                ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTTTTCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                    <  1:1030199/100‑1 (MQ=37)
                ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                    <  2:1201011/100‑1 (MQ=38)
                 gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGTCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcg                                                                                   <  1:3200376/100‑2 (MQ=37)
                 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)
                   tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTg                                                                                 <  1:3573579/100‑1 (MQ=255)
                   tgatgaTAGACATCGAGGATGCCAGCATTTTGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGTGCACCGTAC‑GTGCGCGTgtcgaa                                                                                 <  2:3571268/100‑6 (MQ=255)
                     atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGTCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaaaa                                                                               <  1:1179590/100‑6 (MQ=25)
                     atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGaa                                                                               >  1:36905/1‑100 (MQ=255)
                       gatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTTGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATg                                                                             <  1:2442554/100‑1 (MQ=255)
                         tAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCa                                                                           >  2:3570097/1‑100 (MQ=255)
                           gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGAGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTTCCCGTGACCTGAATGCAg                                                                          >  2:3508311/1‑99 (MQ=255)
                             cAGCGTGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                       <  1:3663895/100‑1 (MQ=255)
                             cAGCGAGGACGCCAGCACTTCGATTTCTCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                       <  1:636010/100‑1 (MQ=255)
                             cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGTTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                       <  1:1594674/100‑1 (MQ=255)
                              aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACt                                                                      <  1:160845/100‑1 (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)
                                      cGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGc                                                              >  1:1337159/1‑100 (MQ=255)
                                      cGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGc                                                              >  1:3576122/1‑100 (MQ=255)
                                       gCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGACCGTAACCTAAATGCAGAACTCATTACGCa                                                             >  1:3195188/1‑100 (MQ=255)
                                        ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                            >  2:2140644/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)
                                               ttCGATTTCGCCGGTCGCCATATCGCGGTTCATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTTagag                                                     <  2:598529/100‑1 (MQ=255)
                                               ttCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                     >  1:482757/1‑100 (MQ=255)
                                                 cGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGc                                                   >  1:1796755/1‑100 (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:1994681/1‑100 (MQ=255)
                                                        gCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCttt                                            >  1:13455/1‑100 (MQ=255)
                                                         ccGGTCGCCATATCGCGGTTAATATTTTTTTTTTTATTCGCA‑AGTACGGAGCCCACGACCTGAATGCAGAAATCTGTACGTAGTTCAGATACCAGCtt                                             >  2:1919062/1‑98 (MQ=255)
                                                           ggTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTa                                           <  2:314311/98‑1 (MQ=255)
                                                            gTCGCCATATCGCGGTTAATATTTTTTTCGTTTTGCGTA‑CGTACTGTGCTCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                        <  1:144851/100‑1 (MQ=255)
                                                            gTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CTTACGGTGTCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                        <  1:2113084/100‑1 (MQ=255)
                                                            gTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                        <  1:2631018/100‑1 (MQ=255)
                                                            gTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACGCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                        <  1:964535/100‑1 (MQ=255)
                                                                ccATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTc                                    >  2:3798042/1‑100 (MQ=255)
                                                                 cATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTcc                                   >  1:2724517/1‑100 (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)
                                                                             aaTATTTTTTTTGTCACGCGCA‑CGTACGTTGCCCGTGACATGAATACAGAACTCATTACGCAGTTCAGATGCGAACTTTAACGCGTCGTCGCGATccag                        >  2:58414/1‑97 (MQ=255)
                                                                               tatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACTTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGTGTTCGCACGATCCGGATc                     <  2:1379872/100‑1 (MQ=255)
                                                                               tatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                     <  2:2695714/100‑1 (MQ=255)
                                                                               tatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAAAGGCAAGCTTTAACGCGTCCGAACGATCCGGATc                     >  1:1748274/1‑100 (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)
                                                                                     tttCGTCACGCGCA‑CGTACTGTGCCCGTGACCTGAATGCAGAATTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaaa               <  1:3698019/100‑1 (MQ=255)
                                                                                     ttacGTCACTCGCA‑AGTTAGGTTATCGTGATCTGAATGCTGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaa                <  1:2210351/96‑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:2398707/2‑100 (MQ=38)
                                                                                        ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc            >  1:3298204/2‑100 (MQ=38)
                                                                                        ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATAGAAAAATAc            >  1:742235/2‑100 (MQ=37)
                                                                                         gTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc           >  1:3363193/1‑100 (MQ=38)
                                                                                         gTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCAGATCGAAAAATAcc           >  1:1938298/1‑100 (MQ=255)
                                                                                          tCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATaaga          >  2:856121/1‑97 (MQ=38)
                                                                                          tCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt          <  1:2140663/100‑1 (MQ=255)
                                                                                          tCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt          >  1:3383599/1‑100 (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)
                                                                                             cgcgCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCa       >  2:3903000/1‑100 (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)
                                                                                                   |                                                                                                   
GGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATA  >  NC_000913/1950197‑1950394

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: ‑