Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,310,176 A→G R81R (CGT→CGC truB ← tRNA pseudouridine 5S synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,310,1760AG99.6% 835.2 / ‑5.0 242R81R (CGT→CGCtruBtRNA pseudouridine 5S synthase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (127/114);  total (128/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 256 total aligned reads displayed.

GCCGCTGCCAGCTGCTCTGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGACGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCGGT  >  NC_000913/3310087‑3310265
                                                                                         |                                                                                         
gCCGCTGCCAGCTGCTCTGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGgcgc                                                                                        >  2:3584566/1‑93 (MQ=255)
gCCGCTGCCAGCTGCTCTGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGgcgc                                                                                        <  2:1165126/93‑1 (MQ=255)
     tGCCAGCTGCTCTGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAAt                                                                                   <  1:2741372/93‑1 (MQ=255)
      gCCAGCTGCTCTGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATg                                                                                  >  1:5562670/1‑93 (MQ=255)
       ccAGCTGCTCTGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGa                                                                                 <  2:842621/93‑1 (MQ=255)
       ccAGCTGCTCTGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGa                                                                                 <  2:3302378/93‑1 (MQ=255)
        cAGCTGCTCTGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGAc                                                                                >  1:1151383/1‑93 (MQ=255)
         aGCTGCTCTGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGaac                                                                               >  2:3827325/1‑91 (MQ=255)
         aGCTGCTCTGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGAcc                                                                               >  1:1587562/1‑93 (MQ=255)
         aGCTGCTCTGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGAcc                                                                               >  2:5172061/1‑93 (MQ=255)
          gctgctCTGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGAccc                                                                              >  2:5880717/1‑93 (MQ=255)
            tgctCTGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGa                                                                            <  2:257899/93‑1 (MQ=255)
            tgctCTGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGa                                                                            <  1:3320186/93‑1 (MQ=255)
              ctctGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATa                                                                          >  2:1123585/1‑93 (MQ=255)
               tctGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAg                                                                         <  1:3599290/93‑1 (MQ=255)
               tctGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAg                                                                         >  1:455686/1‑93 (MQ=255)
                ctGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGc                                                                        <  2:656283/93‑1 (MQ=255)
                 tGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCg                                                                       <  1:6304977/93‑1 (MQ=255)
                  gcgcTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGt                                                                      >  1:5602536/1‑93 (MQ=255)
                    gcTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGttt                                                                    >  1:3831698/1‑93 (MQ=255)
                      tAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGt                                                                  >  2:5668534/1‑93 (MQ=255)
                      tAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGt                                                                  <  1:3381558/93‑1 (MQ=255)
                         aGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCgg                                                               >  2:544575/1‑93 (MQ=255)
                         aGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAACGACCCGATAGCGTTTGTCgg                                                               >  1:350835/1‑93 (MQ=255)
                           gTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAg                                                             >  2:1599059/1‑93 (MQ=255)
                            ttACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGt                                                            >  2:3133372/1‑93 (MQ=255)
                             tACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTc                                                           >  1:4737639/1‑93 (MQ=255)
                              aCCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTcc                                                          >  2:142339/1‑93 (MQ=255)
                              aCCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTcc                                                          >  2:1785512/1‑93 (MQ=255)
                              aCCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTcc                                                          <  2:2772194/93‑1 (MQ=255)
                              aCCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTcc                                                          <  2:6184604/93‑1 (MQ=255)
                              aCCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGCTTGTCGGAGTcc                                                          >  2:2307556/1‑93 (MQ=255)
                              aCCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATCACCCGATAGCGTTTGTCGGAGTcc                                                          <  1:6228274/93‑1 (MQ=255)
                               ccGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCca                                                         <  1:934192/93‑1 (MQ=255)
                               ccGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCca                                                         >  2:4612061/1‑93 (MQ=255)
                               ccGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCca                                                         <  2:4591283/93‑1 (MQ=255)
                                cGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCcag                                                        >  2:5309861/1‑93 (MQ=255)
                                 ggACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCcagc                                                       <  2:4034570/93‑1 (MQ=255)
                                 ggACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCcagc                                                       >  2:3964274/1‑93 (MQ=255)
                                 ggACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGCCcagc                                                       >  2:3675497/1‑93 (MQ=255)
                                   aCGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCcagcag                                                     >  1:5811732/1‑93 (MQ=255)
                                   aCGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCcagcag                                                     <  2:4107264/93‑1 (MQ=255)
                                    cGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGa                                                    <  1:619654/93‑1 (MQ=255)
                                     gTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGAt                                                   >  2:2476857/1‑93 (MQ=255)
                                     gTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGAt                                                   >  1:4550562/1‑93 (MQ=255)
                                     gTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGAt                                                   >  1:3494311/1‑93 (MQ=255)
                                       tcttcAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAACGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATAc                                                 >  2:5618829/1‑93 (MQ=255)
                                        cttcAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACt                                                >  2:1070702/1‑93 (MQ=255)
                                             aCGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGaa                                           >  1:1853894/1‑93 (MQ=255)
                                               gATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGaaaa                                         >  2:1836966/1‑93 (MQ=255)
                                                aTCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAAc                                        <  1:3086608/93‑1 (MQ=255)
                                                aTCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAAc                                        >  1:2222593/1‑93 (MQ=255)
                                                 tCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACt                                       <  1:6054016/93‑1 (MQ=255)
                                                  cTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACtt                                      <  2:6027938/93‑1 (MQ=255)
                                                   tGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGACGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTtc                                     >  2:4979128/1‑93 (MQ=255)
                                                     tCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTtcgt                                   <  2:4291615/93‑1 (MQ=255)
                                                     tCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTtcgt                                   <  2:5851796/93‑1 (MQ=255)
                                                      ccGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTtcgtc                                  <  1:2499734/93‑1 (MQ=255)
                                                      ccGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTtcgtc                                  <  2:6104708/93‑1 (MQ=255)
                                                       cGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTtcgtcg                                 >  2:5097399/1‑93 (MQ=255)
                                                        gTCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGc                                >  2:4491961/1‑93 (MQ=255)
                                                         tCGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCt                               >  1:3975552/1‑93 (MQ=255)
                                                          cGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCtt                              <  2:2351139/93‑1 (MQ=255)
                                                          cGGCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCtt                              <  2:2983203/93‑1 (MQ=255)
                                                            gCATCAGAAGTATCGGTACGCTGTCCAAGGCGTGCAATGACCCGATAGCGTTTGTCGTAGTCCAGCAGATACTGGGAACACTTCGTCGCTTcc                            >  1:5392852/1‑93 (MQ=255)
                                                            gCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTcc                            >  1:1014596/1‑93 (MQ=255)
                                                            gCATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTcc                            >  1:2083345/1‑93 (MQ=255)
                                                             cATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTccc                           <  2:5332516/93‑1 (MQ=255)
                                                             cATCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTccc                           <  1:4874435/93‑1 (MQ=255)
                                                              aTCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGCGTCCAGCAGATCCTGGGCACACTTCGTCGCTTcacc                          >  1:2816857/1‑93 (MQ=255)
                                                              aTCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTcccc                          >  2:3431178/1‑93 (MQ=255)
                                                              aTCAGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTcccc                          <  1:6181865/93‑1 (MQ=255)
                                                                 ggAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAg                       >  2:3250914/2‑93 (MQ=255)
                                                                 aGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAg                       >  2:174722/1‑93 (MQ=255)
                                                                 aGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAg                       >  1:2904079/1‑93 (MQ=255)
                                                                 aGAAGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAg                       >  1:523631/1‑93 (MQ=255)
                                                                   aaGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGc                     >  2:1932439/1‑93 (MQ=255)
                                                                    aGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCa                    >  1:1637127/1‑93 (MQ=255)
                                                                    aGTATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCa                    >  2:6316129/1‑93 (MQ=255)
                                                                      tATCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCaaa                  <  1:1220324/93‑1 (MQ=255)
                                                                       aTCGGTACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAAt                 >  2:2716279/1‑93 (MQ=255)
                                                                            tACGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCa            >  2:4864704/1‑93 (MQ=255)
                                                                             aCGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCaa           <  1:1363362/93‑1 (MQ=255)
                                                                              cGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAAc          >  1:2114426/1‑93 (MQ=255)
                                                                              cGCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAAc          <  2:5742487/93‑1 (MQ=255)
                                                                               gCTGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACa         >  1:4512345/1‑93 (MQ=255)
                                                                                 tGTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATg       <  2:1655449/93‑1 (MQ=255)
                                                                                  ggccAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAACCTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGc      <  1:1742145/91‑1 (MQ=255)
                                                                                  gTCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGc      >  2:4368979/1‑93 (MQ=255)
                                                                                   tCCAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGcc     >  1:1928731/1‑93 (MQ=255)
                                                                                    ccAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCg    >  2:2083078/1‑93 (MQ=255)
                                                                                    ccAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCACCGAGGCAAATCGGCAACATGCCg    >  1:2855516/1‑93 (MQ=255)
                                                                                     cAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCgg   >  1:2664785/1‑93 (MQ=255)
                                                                                     cAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCgg   <  1:421235/93‑1 (MQ=255)
                                                                                     cAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCgg   >  2:6364627/1‑93 (MQ=255)
                                                                                     cAAGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCgg   <  1:5648464/93‑1 (MQ=255)
                                                                                      aaGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCGGt  >  2:5903775/1‑93 (MQ=255)
                                                                                      aaGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCGGt  <  1:3224338/93‑1 (MQ=255)
                                                                                      aaGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCGGt  >  1:4287710/1‑93 (MQ=255)
                                                                                      aaGGCGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCGGt  >  2:646171/1‑93 (MQ=255)
                                                                                         |                                                                                         
GCCGCTGCCAGCTGCTCTGCGCTAAAGGTTACCGGACGTTCTTCAACGATCTGTCCGTCGGCATCAGAAGTATCGGTACGCTGTCCAAGACGCGCAATGACCCGATAGCGTTTGTCGGAGTCCAGCAGATACTGGGAAAACTTCGTCGCTTCCCCGAGGCAAATCGGCAACATGCCGGT  >  NC_000913/3310087‑3310265

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