Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 52,312 G→A A40V (GCG→GTG)  ksgA ← S‑adenosylmethionine‑6‑N',N'‑adenosyl (rRNA) dimethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091352,3120GA99.2% 909.7 / ‑5.6 269A40V (GCG→GTG) ksgAS‑adenosylmethionine‑6‑N',N'‑adenosyl (rRNA) dimethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (89/176);  total (90/179)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

Only 100 of 273 total aligned reads displayed.

TCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCT  >  NC_000913/52222‑52401
                                                                                          |                                                                                         
tCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCAcc                                                                                         <  2:6400036/93‑1 (MQ=255)
 cgcgGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCt                                                                                        <  2:2982848/93‑1 (MQ=255)
  ggggTAAATTTCGATGCCCGTCAGCTGGTCCAGACGCTCGCCGACGGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTg                                                                                       <  2:5330954/91‑1 (MQ=255)
  gcgGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTg                                                                                       <  1:6455897/93‑1 (MQ=255)
  gcgGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTg                                                                                       <  1:3303563/93‑1 (MQ=255)
      tCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGccc                                                                                   <  1:993066/93‑1 (MQ=255)
      tCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGccc                                                                                   <  1:2914246/93‑1 (MQ=255)
      tCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGccc                                                                                   <  2:503587/93‑1 (MQ=255)
      tCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGccc                                                                                   <  1:2282105/93‑1 (MQ=255)
         aGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCttt                                                                                >  2:6505513/1‑93 (MQ=255)
         aGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCttt                                                                                <  1:4708537/93‑1 (MQ=255)
         aGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCttt                                                                                <  1:4331035/93‑1 (MQ=255)
         aGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCttt                                                                                >  2:3257133/1‑93 (MQ=255)
           ttCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTg                                                                              <  2:5489222/93‑1 (MQ=255)
            tCGATGCCCGTCAGCTGGTCCAGGCGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGc                                                                             <  1:5217541/93‑1 (MQ=255)
            tCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGc                                                                             <  1:6582211/93‑1 (MQ=255)
            tCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGc                                                                             <  2:4655043/93‑1 (MQ=255)
            tCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGc                                                                             <  1:306488/93‑1 (MQ=255)
             cGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCg                                                                            <  2:1990134/93‑1 (MQ=255)
              gagtaCCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCgg                                                                           <  2:4936265/89‑1 (MQ=255)
               aTGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCggg                                                                          <  2:3455422/93‑1 (MQ=255)
               aTGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCggg                                                                          <  1:1622541/93‑1 (MQ=255)
                tGNCCGNCAGCTGGTCCAGACGNTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGt                                                                         <  2:5553164/93‑1 (MQ=255)
                tGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGt                                                                         >  1:6734516/1‑93 (MQ=255)
                  aCCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTa                                                                       >  2:1297169/1‑93 (MQ=255)
                     gTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATg                                                                    <  1:500988/93‑1 (MQ=255)
                     gTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATg                                                                    <  1:2304193/93‑1 (MQ=255)
                     ctCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATg                                                                    <  2:5144802/92‑1 (MQ=255)
                      tCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATgg                                                                   <  1:6369555/93‑1 (MQ=255)
                         gCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAg                                                                >  1:1875775/1‑93 (MQ=255)
                         gCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAg                                                                >  2:1856135/1‑93 (MQ=255)
                         gCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAg                                                                >  2:4754148/1‑93 (MQ=255)
                         gCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAg                                                                >  1:5120617/1‑93 (MQ=255)
                          cTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCTGGTTAATGGCAGa                                                               <  1:1177485/93‑1 (MQ=255)
                              tCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGAcaca                                                           <  1:1328543/93‑1 (MQ=255)
                               ccAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACaa                                                          <  1:1125538/93‑1 (MQ=255)
                                  gACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATAc                                                       >  2:4032401/1‑93 (MQ=255)
                                  gACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATAc                                                       <  2:2123851/93‑1 (MQ=255)
                                  gACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATAc                                                       >  1:1037167/1‑93 (MQ=255)
                                  gACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTCTTGCGGGTTAATGGCAGACACAATAc                                                       >  2:6722767/1‑93 (MQ=255)
                                   aCGTCCNCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCCACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACt                                                      <  2:1429249/93‑1 (MQ=255)
                                     gTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGt                                                    <  2:3870254/93‑1 (MQ=255)
                                      ttCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTc                                                   <  1:5563128/93‑1 (MQ=255)
                                      ttCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTc                                                   <  2:2892652/93‑1 (MQ=255)
                                       tCGCCGACCGGTTCGGTCAATGCTGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCg                                                  <  1:4413376/93‑1 (MQ=255)
                                        cGCCGACCGGTTCGGTCAATGCCGCCGGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGa                                                 <  1:5309372/93‑1 (MQ=255)
                                        cGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGa                                                 <  1:3189836/93‑1 (MQ=255)
                                         gCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGAt                                                <  2:2639852/93‑1 (MQ=255)
                                          ccGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATc                                               <  1:5731031/93‑1 (MQ=255)
                                           cGACGGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGCCACAATACTGTGGCTCa                                              <  1:296000/93‑1 (MQ=255)
                                           cGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCa                                              <  2:2508013/93‑1 (MQ=255)
                                            gACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCAc                                             >  1:6035645/1‑93 (MQ=255)
                                            gACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGGTAATGGCAGACACAATACTGTCGATCAc                                             >  2:4532929/1‑93 (MQ=255)
                                            caCCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCAGCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCAc                                             <  2:4322365/92‑1 (MQ=255)
                                             gccGTTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACg                                            <  1:3195933/92‑1 (MQ=255)
                                             aCCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACg                                            >  2:6625948/1‑93 (MQ=255)
                                             aCCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACg                                            >  2:1466433/1‑93 (MQ=255)
                                                ggTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAAc                                         >  2:3428485/1‑93 (MQ=255)
                                                 gTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGCACt                                        >  2:5104005/1‑93 (MQ=255)
                                                     ggTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATc                                    <  1:2469235/93‑1 (MQ=255)
                                                     ggTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATc                                    <  2:3146510/93‑1 (MQ=255)
                                                      gTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCg                                   >  2:4120605/1‑93 (MQ=255)
                                                       tCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGt                                  <  1:2792691/93‑1 (MQ=255)
                                                        cAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGtt                                 <  2:6604130/93‑1 (MQ=255)
                                                        cAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGCCGATCACGAACTGATCGtt                                 <  1:6522199/93‑1 (MQ=255)
                                                         aaTGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTg                                >  1:2694776/1‑93 (MQ=255)
                                                         aaTGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTg                                >  2:3756615/1‑93 (MQ=255)
                                                          aTGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGTTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTga                               <  1:1529075/93‑1 (MQ=255)
                                                             ccgccAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGCACTGATCGTTGAGaa                            >  1:1234649/1‑93 (MQ=255)
                                                              cgccAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGaaa                           <  1:6348596/93‑1 (MQ=255)
                                                               gccAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAg                          <  1:1938327/93‑1 (MQ=255)
                                                               cccAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAg                          <  2:2073535/92‑1 (MQ=255)
                                                                 cAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGtt                        <  1:620309/93‑1 (MQ=255)
                                                                 cAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGtt                        <  1:4804182/93‑1 (MQ=255)
                                                                  aGACCGTGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGttt                       <  2:6113706/93‑1 (MQ=255)
                                                                  aGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGttt                       >  2:3090628/1‑93 (MQ=255)
                                                                   gACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGtttt                      >  1:2491704/1‑93 (MQ=255)
                                                                   gACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGtttt                      >  1:2269771/1‑93 (MQ=255)
                                                                     ccGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGc                    >  2:3594660/1‑93 (MQ=255)
                                                                     ccGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGc                    >  1:6442967/1‑93 (MQ=255)
                                                                      cGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGcc                   <  2:4161015/93‑1 (MQ=255)
                                                                       ggggCGGATTTCGCCCACATCCTGCCCCTTTTGCGGAGTAATGGCAGACATAATACTGTCTATCACGAACTGATCGTTGAGAAAGTTTTAccc                  <  1:3955080/93‑1 (MQ=255)
                                                                       ggggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGccc                  >  1:5407080/1‑93 (MQ=255)
                                                                        gggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCg                 >  1:2855087/1‑93 (MQ=255)
                                                                        gggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCg                 <  1:4122993/93‑1 (MQ=255)
                                                                        gggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCg                 >  1:2090959/1‑93 (MQ=255)
                                                                        gggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCg                 <  1:3479778/93‑1 (MQ=255)
                                                                          gCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGaa               >  2:6319400/1‑93 (MQ=255)
                                                                            cGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGc             <  2:113380/93‑1 (MQ=255)
                                                                             gATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCg            >  2:5966412/1‑93 (MQ=255)
                                                                              aTTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGt           <  2:1561704/93‑1 (MQ=255)
                                                                              aTTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGt           >  2:1707718/1‑93 (MQ=255)
                                                                                ttCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGttt         >  2:4230513/1‑93 (MQ=255)
                                                                                 tCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTa        <  2:2719404/93‑1 (MQ=255)
                                                                                  cGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTAc       <  1:790829/93‑1 (MQ=255)
                                                                                   gACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACg      >  2:1240645/1‑93 (MQ=255)
                                                                                    aCCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACgg     <  2:819246/93‑1 (MQ=255)
                                                                                      cATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGc   >  1:358391/1‑93 (MQ=255)
                                                                                       aTCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCt  <  2:5804357/93‑1 (MQ=255)
                                                                                       aTCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCt  >  1:1827989/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCT  >  NC_000913/52222‑52401

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