Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 426,471 G→A 69.5% R112H (CGT→CAT)  tgt → tRNA‑guanine transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913426,4710GA69.5% 212.0 / 124.8 160R112H (CGT→CAT) tgttRNA‑guanine transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (30/19);  new base A (58/53);  total (88/72)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 161 total aligned reads displayed.

TGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTTGATGAGTG  >  NC_000913/426375‑426570
                                                                                                |                                                                                                   
tGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTaa                                                                                                  >  1:2192315/1‑100 (MQ=255)
tGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTaa                                                                                                  <  2:2038011/100‑1 (MQ=255)
tGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATaa                                                                                                  >  2:2427325/1‑100 (MQ=255)
 gCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAAc                                                                                                 <  2:1816334/100‑1 (MQ=255)
  cAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGCTATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAAcc                                                                                                <  2:4378344/100‑1 (MQ=255)
   aGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAAccc                                                                                               <  2:1237359/100‑1 (MQ=255)
   aGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCAAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAAccc                                                                                               <  1:1468972/100‑1 (MQ=255)
      ggAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAACTCACCGAACAGGGCGTGCACTTCCATAACCCGAt                                                                                            >  2:2581137/1‑100 (MQ=255)
      ggAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGAt                                                                                            >  1:1865339/1‑100 (MQ=255)
       gAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATc                                                                                           >  1:3052886/1‑100 (MQ=255)
        aaGGGGCCTATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCa                                                                                          <  1:190683/100‑1 (MQ=255)
         aGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCaa                                                                                         >  1:2242948/1‑100 (MQ=255)
          ggggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAAc                                                                                        >  2:1617026/1‑100 (MQ=255)
            ggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACgg                                                                                      >  2:2728994/1‑100 (MQ=255)
            ggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACgg                                                                                      >  2:1092339/1‑100 (MQ=255)
             gCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGc                                                                                     <  2:3403264/100‑1 (MQ=255)
             gCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGc                                                                                     >  2:3015102/1‑100 (MQ=255)
               cGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGa                                                                                   <  2:1569361/100‑1 (MQ=255)
               cGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGa                                                                                   >  1:2521366/1‑100 (MQ=255)
                 aTCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATc                                                                                 <  1:1138405/100‑1 (MQ=255)
                   ccTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCg                                                                               >  1:2479450/1‑100 (MQ=255)
                    cTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGa                                                                              >  2:3586697/1‑100 (MQ=255)
                      cACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGAtt                                                                            >  2:2995960/1‑100 (MQ=255)
                       aCCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGAttt                                                                           >  2:1299876/1‑100 (MQ=255)
                       aCCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGAttt                                                                           >  2:3675595/1‑100 (MQ=255)
                        ccGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGAtttt                                                                          >  1:306043/1‑100 (MQ=255)
                        ccGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGAttt                                                                           <  2:1284106/99‑1 (MQ=255)
                          gACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTcc                                                                        <  1:4354073/100‑1 (MQ=255)
                          gACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTcc                                                                        <  2:3836662/100‑1 (MQ=255)
                          gACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTcc                                                                        >  1:3716529/1‑100 (MQ=255)
                            cTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTc                                                                      <  2:737354/100‑1 (MQ=255)
                            cTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTc                                                                      <  2:3893636/100‑1 (MQ=255)
                            cTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTc                                                                      <  1:2244895/100‑1 (MQ=255)
                            cTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTc                                                                      <  1:3305021/100‑1 (MQ=255)
                              ccGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGa                                                                    >  2:1282565/1‑100 (MQ=255)
                              ccGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGa                                                                    >  2:3467785/1‑100 (MQ=255)
                                ggcggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATc                                                                  <  2:401707/100‑1 (MQ=255)
                                ggcggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATc                                                                  >  1:4242178/1‑100 (MQ=255)
                                ggcggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATc                                                                  <  2:1543173/100‑1 (MQ=255)
                                 gcggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATcc                                                                 >  1:2997324/1‑100 (MQ=255)
                                 gcggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATcc                                                                 <  1:3647183/100‑1 (MQ=255)
                                 gcggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATcc                                                                 <  2:3995577/100‑1 (MQ=255)
                                   ggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTg                                                               <  2:565105/100‑1 (MQ=255)
                                   ggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTg                                                               >  2:2907909/1‑100 (MQ=255)
                                    acTTCCAGGTCTTCAGCCGTGGCGTGATTCGTAAAATCACCGCGCAGGGCGTTCACTACCGTAAGCCGACGAACGGCGATCCGATTATCCTCGATCCTGa                                                              >  1:3758741/2‑100 (MQ=255)
                                     cTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaa                                                             >  2:610206/1‑100 (MQ=255)
                                         cAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAAt                                                         >  2:1312131/1‑100 (MQ=255)
                                         cAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAAt                                                         <  2:2507775/100‑1 (MQ=255)
                                         cAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAAt                                                         <  2:1937465/100‑1 (MQ=255)
                                            gTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCaa                                                      <  1:806365/100‑1 (MQ=255)
                                             tCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAAt                                                     <  1:598605/100‑1 (MQ=255)
                                              cTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATg                                                    <  1:2835582/100‑1 (MQ=255)
                                                tCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGga                                                  >  2:3306995/1‑100 (MQ=255)
                                                tCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCAATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAACATCAATGga                                                  >  2:812479/1‑100 (MQ=255)
                                                 cAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGgag                                                 >  2:3118803/1‑100 (MQ=255)
                                                   gCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGAt                                               <  1:3141291/100‑1 (MQ=255)
                                                   gCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGAt                                               >  2:63807/1‑100 (MQ=255)
                                                   gCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGAt                                               >  1:3947273/1‑100 (MQ=255)
                                                    ccTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGgaga                                                >  2:1114816/1‑98 (MQ=255)
                                                     cTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTc                                             >  1:3238524/1‑100 (MQ=255)
                                                        ggCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGt                                          <  2:1613678/100‑1 (MQ=255)
                                                         gCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTa                                         >  2:1951515/1‑100 (MQ=255)
                                                         gCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTa                                         <  1:13838/100‑1 (MQ=255)
                                                             tatTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGAt                                     >  2:2728979/1‑100 (MQ=255)
                                                                tCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCtt                                  >  1:2199795/1‑100 (MQ=255)
                                                                 cGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTg                                 >  1:206922/1‑100 (MQ=255)
                                                                  gTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTgg                                >  2:1063859/1‑100 (MQ=255)
                                                                  gTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTgg                                <  1:2941727/100‑1 (MQ=255)
                                                                      aaTCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCg                            <  2:2479462/100‑1 (MQ=255)
                                                                      aaTCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCg                            >  1:4218361/1‑100 (MQ=255)
                                                                      aaTCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCg                            <  1:4093999/100‑1 (MQ=255)
                                                                       aTCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCgg                           <  2:2908466/100‑1 (MQ=255)
                                                                       aTCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCgg                           >  2:3111927/1‑100 (MQ=255)
                                                                        tCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGa                          <  1:1090022/100‑1 (MQ=255)
                                                                         cACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGat                         <  2:1456470/100‑1 (MQ=255)
                                                                         cACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGat                         >  1:2020372/1‑100 (MQ=255)
                                                                         cACCGAACAGGGCATGAACTTCCATAACCCGATCAACGGCTATCCGAAAAAAAAAAAAAAAAAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGat                         <  1:4069834/100‑1 (MQ=255)
                                                                          aCCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGata                        >  1:3558968/1‑100 (MQ=255)
                                                                           ccGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGatat                       >  2:3136697/1‑100 (MQ=255)
                                                                            cGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATc                      >  2:4083950/1‑100 (MQ=255)
                                                                             gAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAACTCAATGGAGATTCAGTACGATCTTGGTTCGGATATCg                     >  2:876603/1‑100 (MQ=255)
                                                                             gAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCg                     >  1:4403735/1‑100 (MQ=255)
                                                                              aaCAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGt                    <  1:1628031/100‑1 (MQ=255)
                                                                                cAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAACATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCa                  >  1:2172451/1‑100 (MQ=255)
                                                                                cAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCa                  <  2:3053525/100‑1 (MQ=255)
                                                                                  ggggGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCAt                 <  1:1575055/99‑1 (MQ=255)
                                                                                  gggCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATg                >  2:4296687/1‑100 (MQ=255)
                                                                                  gggCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATg                >  2:781236/1‑100 (MQ=255)
                                                                                       tGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCtt           >  2:2564800/1‑100 (MQ=255)
                                                                                       tGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCtt           >  2:629364/1‑100 (MQ=255)
                                                                                        gCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCttt          <  2:3253057/100‑1 (MQ=255)
                                                                                         cACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtg         >  2:202413/1‑100 (MQ=255)
                                                                                          aCTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtga        >  1:3336421/1‑100 (MQ=255)
                                                                                           cTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgat       >  2:1265456/1‑100 (MQ=255)
                                                                                           cTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtg         <  1:1824120/98‑1 (MQ=255)
                                                                                           cTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgat       <  2:3947303/100‑1 (MQ=255)
                                                                                            ttCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgatg      <  1:1766784/100‑1 (MQ=255)
                                                                                            ttCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgatg      >  2:2797/1‑100 (MQ=255)
                                                                                            ttCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgatg      <  2:3712873/100‑1 (MQ=255)
                                                                                             tCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgatga     >  2:269134/1‑100 (MQ=255)
                                                                                                |                                                                                                   
TGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTTGATGAGTG  >  NC_000913/426375‑426570

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