Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 426,471 G→A 69.6% 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.6% 125.0 / 70.4 92R112H (CGT→CAT) tgttRNA‑guanine transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (19/9);  new base A (40/24);  total (59/33)
Fisher's exact test for biased strand distribution p-value = 8.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTTGATGAG  >  NC_000913/426374‑426568
                                                                                                 |                                                                                                 
aTGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATa                                                                                                 <  1:2043050/100‑1 (MQ=255)
 tGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATaa                                                                                                >  2:310023/1‑100 (MQ=255)
   cAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAAcc                                                                                              >  2:1020312/1‑100 (MQ=255)
   cAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAAcc                                                                                              >  2:256909/1‑100 (MQ=255)
   cAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAAcc                                                                                              >  2:2128157/1‑100 (MQ=255)
     gTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCg                                                                                            >  2:1321186/1‑100 (MQ=255)
        gAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATc                                                                                         <  2:1294747/100‑1 (MQ=255)
        gAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATc                                                                                         <  1:1815482/100‑1 (MQ=255)
          aGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCaa                                                                                       >  2:2394913/1‑100 (MQ=255)
          aGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCaa                                                                                       <  1:461470/100‑1 (MQ=255)
          aGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCaa                                                                                       >  1:2105989/1‑100 (MQ=255)
              gCCGATCCTCACCGACCCCGGCGGCTTCCAGGTCTTCAGCGTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACggn                                                                                   <  1:1235610/100‑2 (MQ=255)
                  aTCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATc                                                                               <  2:1456368/100‑1 (MQ=255)
                  aTCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATc                                                                               <  2:620452/100‑1 (MQ=255)
                   tCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATcc                                                                              <  2:885322/100‑1 (MQ=255)
                    ccTCACCGACTACGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGACGTGCACTTCCATAACCCGATCAACGGCGATCCg                                                                             >  2:1396280/1‑100 (MQ=255)
                     cTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGa                                                                            >  1:2174384/1‑100 (MQ=255)
                     cTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGa                                                                            >  1:2453586/1‑100 (MQ=255)
                      tCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGAt                                                                           >  1:661938/1‑100 (MQ=255)
                      tCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGACGTGCACTTCCATAACCCGATCAACGGCGATCCGa                                                                            >  1:809290/1‑99 (MQ=255)
                       cACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGAtt                                                                          >  2:2135960/1‑100 (MQ=255)
                        aCCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGAttt                                                                         >  1:2169693/1‑100 (MQ=255)
                         ccGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGAtttt                                                                        >  2:1252542/1‑100 (MQ=255)
                          cGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTACCCCGATCAACGGCGATCCGATTTTc                                                                       >  2:2440142/1‑100 (MQ=255)
                          cGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGAtttt                                                                        >  1:1327052/1‑99 (MQ=255)
                           gACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTcc                                                                      >  2:458512/1‑100 (MQ=255)
                           gACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTcc                                                                      >  1:2195686/1‑100 (MQ=255)
                           gACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTcc                                                                      >  2:863947/1‑100 (MQ=255)
                             cTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTc                                                                    <  1:2128183/100‑1 (MQ=255)
                             cTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCATAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTc                                                                    <  1:2384454/100‑1 (MQ=255)
                               ccGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGa                                                                  >  1:742909/1‑100 (MQ=255)
                               ccGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGa                                                                  >  2:1531245/1‑100 (MQ=255)
                               ccGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGa                                                                  >  2:1714556/1‑100 (MQ=255)
                                cggcggCTTCCAGGTCTTCAGCCTTGGTGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGAt                                                                 <  1:322664/100‑1 (MQ=255)
                                 ggcggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATc                                                                >  1:1644974/1‑100 (MQ=255)
                                      cTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaa                                                           >  2:2339840/1‑100 (MQ=255)
                                      cTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaa                                                           >  2:1042791/1‑100 (MQ=255)
                                       ttCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaaa                                                          >  1:2489040/1‑100 (MQ=255)
                                       ttCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaaa                                                          >  1:1603650/1‑100 (MQ=255)
                                         ccAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaaaaa                                                        >  2:386524/1‑100 (MQ=255)
                                          cAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAAt                                                       >  1:1654686/1‑100 (MQ=255)
                                           aGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATc                                                      <  1:1403251/100‑1 (MQ=255)
                                           aGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATc                                                      >  2:900023/1‑100 (MQ=255)
                                            ggTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATAAACGGCGATACGATTTTCCTCGATCCTGAAAAATCa                                                     >  1:986545/1‑100 (MQ=255)
                                             gTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCaa                                                    >  1:2133754/1‑100 (MQ=255)
                                              tCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATc                                                      >  2:1758650/1‑97 (MQ=255)
                                                ttCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATgg                                                 >  1:1244537/1‑100 (MQ=255)
                                                 tCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGga                                                >  2:1519606/1‑100 (MQ=255)
                                                  cAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGTTTTTCCTCGATCCTGAAAAATCAATGgag                                               >  2:1587939/1‑100 (MQ=255)
                                                    gCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGAt                                             >  1:523002/1‑100 (MQ=255)
                                                    gCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGAt                                             <  1:458516/100‑1 (MQ=255)
                                                      cTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTc                                           <  2:581964/100‑1 (MQ=255)
                                                      cTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTc                                           >  2:1402500/1‑100 (MQ=255)
                                                      cTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTc                                           >  1:36143/1‑100 (MQ=255)
                                                      cTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTc                                           >  1:1884409/1‑100 (MQ=255)
                                                      cTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTc                                           <  1:863956/100‑1 (MQ=255)
                                                          gCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTa                                       <  2:2195660/100‑1 (MQ=255)
                                                           cGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTAc                                      <  2:934924/100‑1 (MQ=255)
                                                           cGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTAc                                      >  2:790711/1‑100 (MQ=255)
                                                            gATATTCGTAAAATCACCGAACAGGGCGTGCCCTTCCGTAACCCGATCAACGACGATCAGATTTTCATCAATCCTGAAAAATCAATGGAGATTCAGt                                        >  1:852202/1‑97 (MQ=255)
                                                            gATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACg                                     <  2:1244525/100‑1 (MQ=255)
                                                               atTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGa                                    >  1:2580492/1‑98 (MQ=255)
                                                                 tCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCtt                                >  2:211228/1‑100 (MQ=255)
                                                                 tCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCtt                                >  1:2474918/1‑100 (MQ=255)
                                                                  cGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTg                               <  2:1602579/100‑1 (MQ=255)
                                                                   gTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTgg                              >  2:1707898/1‑100 (MQ=255)
                                                                   gTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTgg                              >  2:174925/1‑100 (MQ=255)
                                                                    tAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGt                             >  2:2074107/1‑100 (MQ=255)
                                                                     aaaaTCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGtt                            <  1:1714572/100‑1 (MQ=255)
                                                                          cACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGat                       <  1:2339870/100‑1 (MQ=255)
                                                                          cACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGat                       >  1:1430808/1‑100 (MQ=255)
                                                                            ccGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGatat                     <  2:48741/100‑1 (MQ=255)
                                                                             cGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATc                    >  2:491274/1‑100 (MQ=255)
                                                                              gAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCg                   >  2:1132541/1‑100 (MQ=255)
                                                                               aaCAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGt                  <  2:95622/100‑1 (MQ=255)
                                                                               aaCAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGt                  >  1:1419999/1‑100 (MQ=255)
                                                                               aaCAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGt                  <  1:2255970/100‑1 (MQ=255)
                                                                                 cAGGGCGTGCACTTCCATAACCCGATCAATGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCa                <  2:2635954/100‑1 (MQ=255)
                                                                                 cAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCa                >  1:700940/1‑100 (MQ=255)
                                                                                     gCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGAt            >  2:1183316/1‑100 (MQ=255)
                                                                                     gCGTGCACTTCCATAACCCGATCAACGGCGATTCGATTTTCCTTGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGAt            <  1:734611/100‑1 (MQ=255)
                                                                                      cGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATc           >  2:361785/1‑100 (MQ=255)
                                                                                      cGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATc           <  1:386526/100‑1 (MQ=255)
                                                                                      cGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATc           <  1:286039/100‑1 (MQ=255)
                                                                                         gCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCttt        <  1:1531260/100‑1 (MQ=255)
                                                                                         gCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCttt        <  1:293696/100‑1 (MQ=255)
                                                                                           aCTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtga      <  1:1110017/100‑1 (MQ=255)
                                                                                            cTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTAAAGATCTTtg       >  2:1982158/1‑98 (MQ=255)
                                                                                            cTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgat     >  2:1513278/1‑100 (MQ=255)
                                                                                            cTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGAAAATGAAAAATCAATGGAGATTCAGTACGATCTTGGATCGGATATCGTCATGTTCTTtg       >  2:308238/1‑98 (MQ=255)
                                                                                             ttCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgatg    <  1:704982/100‑1 (MQ=255)
                                                                                               ccGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTTGATGAg  <  1:46745/100‑1 (MQ=255)
                                                                                                 |                                                                                                 
ATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTTGATGAG  >  NC_000913/426374‑426568

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