Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 426,471 G→A 31.3% 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,4710GA31.3% 120.0 / 67.2 89R112H (CGT→CAT) tgttRNA‑guanine transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base A (14/14);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

TTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTTGATGA  >  NC_000913/426372‑426567
                                                                                                   |                                                                                                
ttATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCg                                                                                                  <  1:1821685/100‑1 (MQ=255)
 tATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGt                                                                                                 <  1:891450/100‑1 (MQ=255)
   tGCTGTGGAAGGGGCCGATCATGACCGACTCCGAGGAATTTCAGGTCTTCAGCCTTCGCGATATTCGTAAAATAACTGAATAGGGCGTGCACTTCCGTa                                                                                                >  2:773215/1‑99 (MQ=255)
   tGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATaa                                                                                               >  1:1834806/1‑100 (MQ=255)
    gCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAAc                                                                                              <  1:947491/100‑1 (MQ=255)
    gCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAAc                                                                                              >  1:1291112/1‑100 (MQ=255)
     cAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAAcc                                                                                             >  1:1662930/1‑100 (MQ=255)
     cAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAAc                                                                                              >  1:203612/1‑99 (MQ=255)
       gTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCg                                                                                           >  2:481103/1‑100 (MQ=255)
        tGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATCTTCGTAAAATCACCGAACAGGGCGTGCATTTCCGTAAcccca                                                                                          <  1:1517262/100‑3 (MQ=255)
         ggAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGAt                                                                                         <  1:659224/100‑1 (MQ=255)
         ggAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGAt                                                                                         >  2:303217/1‑100 (MQ=255)
          gAAGGGGCCGGTCCTCACCGGCTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATTTTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAAccc                                                                                            <  2:459166/96‑1 (MQ=255)
          gAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATc                                                                                        >  1:1767261/1‑100 (MQ=255)
             ggggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAAc                                                                                     >  2:372056/1‑100 (MQ=255)
             ggggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAAc                                                                                     <  1:854993/100‑1 (MQ=255)
               ggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAAc                                                                                     >  1:813247/1‑98 (MQ=255)
               ggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACgg                                                                                   <  2:2166065/100‑1 (MQ=255)
               ggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACgg                                                                                   >  2:116083/1‑100 (MQ=255)
               ggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACg                                                                                    >  1:1710063/1‑99 (MQ=255)
                 ccGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCg                                                                                 <  2:3802/100‑1 (MQ=255)
                    aTCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATc                                                                              >  2:2058904/1‑100 (MQ=255)
                     tctttaCCGACTCCGGCGGTTTTCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATcc                                                                             <  1:223977/95‑1 (MQ=255)
                        tCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGAt                                                                          <  2:2082572/100‑1 (MQ=255)
                        tCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTCCCATAACCCGATCAACGGCGATCCGAt                                                                          >  2:623892/1‑100 (MQ=255)
                          aCCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGAttt                                                                        <  2:1405618/100‑1 (MQ=255)
                             gACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTcc                                                                     >  1:1800447/1‑100 (MQ=255)
                             gACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTcc                                                                     <  1:2205999/100‑1 (MQ=255)
                             gACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTcc                                                                     <  2:1202880/100‑1 (MQ=255)
                               cTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTc                                                                   >  1:644873/1‑100 (MQ=255)
                               cTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTc                                                                   <  1:2213950/100‑1 (MQ=255)
                                 ccGGCGGCTTCTAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGa                                                                 <  2:1835937/100‑1 (MQ=255)
                                 ccGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGa                                                                 >  1:1410333/1‑100 (MQ=255)
                                 ccGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGa                                                                 >  2:69305/1‑100 (MQ=255)
                                   gtcggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATc                                                               >  1:2079970/3‑100 (MQ=255)
                                   ggcggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATc                                                               <  1:2114074/100‑1 (MQ=255)
                                   ggcggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATc                                                               >  1:1463086/1‑100 (MQ=255)
                                   ggcggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGAt                                                                <  1:253624/99‑1 (MQ=255)
                                    gcggcTTCCAGGTGTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATcc                                                              <  2:1556390/100‑1 (MQ=255)
                                      ggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTg                                                            >  2:1226989/1‑100 (MQ=255)
                                       gcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGa                                                           >  2:897918/1‑100 (MQ=255)
                                        cTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaa                                                          >  1:1129720/1‑100 (MQ=255)
                                          tCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaaaa                                                        >  2:219961/1‑100 (MQ=255)
                                          tCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaaaa                                                        <  1:1453437/100‑1 (MQ=255)
                                            cAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAAt                                                      <  1:261161/100‑1 (MQ=255)
                                             aGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATc                                                     >  2:490226/1‑100 (MQ=255)
                                             aGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATc                                                     <  2:1077570/100‑1 (MQ=255)
                                              ggTCTTCAGCCTTGGCGATATTCGTAAAATCTCCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAAt                                                      <  1:734811/98‑1 (MQ=255)
                                              ggTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCa                                                    <  1:935502/100‑1 (MQ=255)
                                              ggTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCa                                                    >  1:1633351/1‑100 (MQ=255)
                                              ggTCTTCAGACTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCa                                                    <  2:2058029/100‑1 (MQ=255)
                                              ggTCGTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATc                                                     <  1:1054834/99‑1 (MQ=255)
                                                  ttCAGCCTTGGCGATATTCGTAAAATCACCGAACAGCGCGTGCACTTCCGTAACCCGATCAACGTCGATCCGATTTTCCTCGATCCTGAAAAATCAATgg                                                <  2:136293/100‑1 (MQ=255)
                                                   tCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGga                                               >  1:1007626/1‑100 (MQ=255)
                                                      gCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGAt                                            >  1:1479836/1‑100 (MQ=255)
                                                      gCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGAt                                            >  2:745422/1‑100 (MQ=255)
                                                       ccTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGgaga                                             >  1:624110/1‑98 (MQ=255)
                                                        cTTGGCGATATTCGTAAAATCACCGAACTGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTc                                          <  2:2079877/100‑1 (MQ=255)
                                                         ttGGCTATATTCGTATAATCACCTAACAGGGCGTGCTCTTCCGTAACCCGATCAACGGCGTTCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCa                                         <  2:17905/100‑1 (MQ=255)
                                                             cGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTAc                                     >  1:855047/1‑100 (MQ=255)
                                                                  ttCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCt                                >  2:1607116/1‑100 (MQ=255)
                                                                  ttCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCt                                >  1:2036365/1‑100 (MQ=255)
                                                                   tCGAAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTAAAAAAAAAAAAAAAATCAATGGAGATTCAGTACGATCtt                               <  2:1633340/100‑1 (MQ=255)
                                                                    cGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTg                              <  1:984518/100‑1 (MQ=255)
                                                                     gTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTgg                             >  1:997090/1‑100 (MQ=255)
                                                                     gTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTgg                             <  1:1091657/100‑1 (MQ=255)
                                                                      tAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATTGAGATTCAGTACGATCTTgg                             <  2:112146/99‑1 (MQ=255)
                                                                          aTCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCgg                        <  1:1457845/100‑1 (MQ=255)
                                                                          aTCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCg                         <  1:243325/99‑1 (MQ=255)
                                                                          aTCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGTAGATTCAGTACGATCTTGGTTCgn                        <  2:2077517/100‑2 (MQ=255)
                                                                                gAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCg                  <  2:813247/100‑1 (MQ=255)
                                                                                gAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCg                  >  2:23430/1‑100 (MQ=255)
                                                                                gAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCg                  <  2:1742311/100‑1 (MQ=255)
                                                                                 aaCAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGt                 <  1:2175232/100‑1 (MQ=255)
                                                                                 aaCAAGGCGAGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGAAGATTCAGTACGATCTTGGTTCGGATATCGt                 <  2:1669055/100‑1 (MQ=255)
                                                                                     gggCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATg             >  2:394870/1‑100 (MQ=255)
                                                                                          tGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCtt        >  1:503354/1‑100 (MQ=255)
                                                                                          tGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCtt        >  1:1238790/1‑100 (MQ=255)
                                                                                            cACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtg      <  1:1437781/100‑1 (MQ=255)
                                                                                             aCTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtga     <  2:324168/100‑1 (MQ=255)
                                                                                             aCTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtga     <  2:721626/100‑1 (MQ=255)
                                                                                              cTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgat    <  1:372055/100‑1 (MQ=255)
                                                                                              cTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgat    >  1:1864373/1‑100 (MQ=255)
                                                                                              cTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgat    <  2:1662919/100‑1 (MQ=255)
                                                                                               ttCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgatg   >  1:1447250/1‑100 (MQ=255)
                                                                                               ttCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgatg   <  1:1095753/100‑1 (MQ=255)
                                                                                               ttCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgatg   <  1:481930/100‑1 (MQ=255)
                                                                                               ttCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgatg   <  2:1487037/100‑1 (MQ=255)
                                                                                               ttCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgatg   >  1:2016572/1‑100 (MQ=255)
                                                                                                tCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgatga  >  2:1976317/1‑100 (MQ=255)
                                                                                                   |                                                                                                
TTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTTGATGA  >  NC_000913/426372‑426567

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: