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,3120GA98.7% 1006.3 / ‑4.9 303A40V (GCG→GTG) ksgAS‑adenosylmethionine‑6‑N',N'‑adenosyl (rRNA) dimethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (118/178);  total (119/184)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

Only 100 of 308 total aligned reads displayed.

TCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCT  >  NC_000913/52222‑52401
                                                                                          |                                                                                         
tCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCAcc                                                                                         <  2:3885199/93‑1 (MQ=255)
tCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCAcc                                                                                         <  2:3700924/93‑1 (MQ=255)
tCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCAcc                                                                                         <  1:298316/93‑1 (MQ=255)
 cgcgGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCt                                                                                        >  2:1290454/1‑93 (MQ=255)
 cgcgGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCt                                                                                        <  2:3570663/93‑1 (MQ=255)
      tCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGccc                                                                                   >  1:4184743/1‑93 (MQ=255)
      tCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGccc                                                                                   >  2:2022481/1‑93 (MQ=255)
           ttCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTg                                                                              >  2:485333/1‑93 (MQ=255)
            tCGATGACCGTCAGCTGGTCCAGACGTTCGCCGATCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGc                                                                             >  2:4954793/1‑93 (MQ=255)
            tCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGc                                                                             <  2:4762388/93‑1 (MQ=255)
               aTGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCggg                                                                          <  2:3053877/93‑1 (MQ=255)
               aTGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCggg                                                                          >  2:3499981/1‑93 (MQ=255)
               aTGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCggg                                                                          <  2:7301500/93‑1 (MQ=255)
               aTGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGACAGACCGGGGCCGATTTCGACCATCACCTGGCCCTGTTGCggg                                                                          >  1:7202389/1‑93 (MQ=255)
                 gACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGtt                                                                        <  2:5795322/93‑1 (MQ=255)
                  aCCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTa                                                                       >  1:190709/1‑93 (MQ=255)
                   ccGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTaa                                                                      >  1:3354803/1‑93 (MQ=255)
                   ccGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTaa                                                                      <  2:2495829/93‑1 (MQ=255)
                   ccGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTaa                                                                      <  1:4080099/93‑1 (MQ=255)
                    cGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAAt                                                                     <  1:7764069/93‑1 (MQ=255)
                    cGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAAt                                                                     >  1:4322560/1‑93 (MQ=255)
                     gTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATg                                                                    >  1:620785/1‑93 (MQ=255)
                      tCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATgg                                                                   >  2:7081776/1‑93 (MQ=255)
                      gcAGCGGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATgg                                                                   <  1:1611198/92‑1 (MQ=255)
                       cAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGc                                                                  <  2:5636824/93‑1 (MQ=255)
                        aGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCa                                                                 <  1:7451554/93‑1 (MQ=255)
                        aGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCa                                                                 >  2:1070425/1‑93 (MQ=255)
                        aGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCa                                                                 <  1:4984778/93‑1 (MQ=255)
                        aGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCa                                                                 <  2:5345168/93‑1 (MQ=255)
                        aGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCa                                                                 >  1:935680/1‑93 (MQ=255)
                         gCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAg                                                                >  2:3199072/1‑93 (MQ=255)
                           tGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGac                                                              >  2:2215584/1‑93 (MQ=255)
                            ggTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGaca                                                             <  2:1723891/93‑1 (MQ=255)
                             gTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGacac                                                            <  1:5760800/93‑1 (MQ=255)
                               ccAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACaa                                                          >  1:7547644/1‑93 (MQ=255)
                               ccAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACaa                                                          <  2:7722211/93‑1 (MQ=255)
                                cAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAAt                                                         <  2:6719350/93‑1 (MQ=255)
                                cAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAAt                                                         <  1:2315187/93‑1 (MQ=255)
                                 aggcGTTCGCCGACCGGTTCGTTCAATGCCGCCAGACCGGGGCCAATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATa                                                        <  1:59795/90‑1 (MQ=255)
                                  gACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATAc                                                       <  1:3749831/93‑1 (MQ=255)
                                  gACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATAc                                                       <  2:5112249/93‑1 (MQ=255)
                                  gACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATAc                                                       >  2:7812046/1‑93 (MQ=255)
                                     gTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGt                                                    >  1:4438176/1‑93 (MQ=255)
                                      ttCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTc                                                   <  2:3978038/93‑1 (MQ=255)
                                        cGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGa                                                 <  2:6605010/93‑1 (MQ=255)
                                          ccGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATc                                               <  2:288591/93‑1 (MQ=255)
                                          ccGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATc                                               <  2:7775999/93‑1 (MQ=255)
                                          ccGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATc                                               <  1:1340161/93‑1 (MQ=255)
                                           cGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCa                                              <  1:5434862/93‑1 (MQ=255)
                                           cGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCa                                              <  1:1774007/93‑1 (MQ=255)
                                            gACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCAc                                             <  1:4725302/93‑1 (MQ=255)
                                              ccGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGa                                           <  2:5692271/93‑1 (MQ=255)
                                              ccGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAAAGGCAGACACAATACTGTCGATCACGa                                           >  1:6756026/1‑93 (MQ=255)
                                               cGGTTCGGTCAATGCCGCCAGACAGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGaa                                          <  2:4310545/93‑1 (MQ=255)
                                                ggTTCGTTCAATGCCGCCAGACCGGGGCTGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATTCTGTCGATCACGAAc                                         <  2:746921/93‑1 (MQ=255)
                                                ggTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGCAc                                         >  1:1047305/1‑93 (MQ=255)
                                                ggTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAAc                                         <  2:4143582/93‑1 (MQ=255)
                                                ggTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAAc                                         <  2:4407056/93‑1 (MQ=255)
                                                ggTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAAc                                         >  1:1456350/1‑93 (MQ=255)
                                                 gTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCTCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACt                                        <  1:4854974/93‑1 (MQ=255)
                                                     ggTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATc                                    <  2:1443407/93‑1 (MQ=255)
                                                     ggTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATc                                    >  2:6989538/1‑93 (MQ=255)
                                                      gTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCg                                   <  1:7147434/93‑1 (MQ=255)
                                                      gTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCg                                   >  2:6488920/1‑93 (MQ=255)
                                                        cAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGtt                                 >  1:859080/1‑93 (MQ=255)
                                                        cAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGtt                                 <  1:3059050/93‑1 (MQ=255)
                                                         aaTGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTg                                <  2:5454085/93‑1 (MQ=255)
                                                            gccgccAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTgaga                             <  1:4611693/93‑1 (MQ=255)
                                                              cgccAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGaaa                           <  1:7873901/93‑1 (MQ=255)
                                                              cgccAGACAGGGACAGATTTCGACCATCACCTGGCCCTTTTGAGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGaaa                           <  1:7003336/93‑1 (MQ=255)
                                                               gccAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGGTAACGAACTGATCGTTGAGAAAg                          >  2:5376867/1‑93 (MQ=255)
                                                                  aGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGttt                       >  1:1127130/1‑93 (MQ=255)
                                                                   gACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGtttt                      <  2:2684505/93‑1 (MQ=255)
                                                                   gACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGtttt                      <  2:6077970/93‑1 (MQ=255)
                                                                      cGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGCCACAATACTGTCGATGCCGACCTCGTCGTTCACAACGTTTTGcc                   >  1:4867835/1‑93 (MQ=255)
                                                                       ggggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGccc                  <  2:1872697/93‑1 (MQ=255)
                                                                          gCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGaa               <  2:805220/93‑1 (MQ=255)
                                                                          gCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGaa               <  1:3553793/93‑1 (MQ=255)
                                                                          gCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGaa               >  1:3862331/1‑93 (MQ=255)
                                                                          gCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGaa               <  1:4492155/93‑1 (MQ=255)
                                                                          gCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGaa               <  2:1510579/93‑1 (MQ=255)
                                                                           ccGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAg              >  2:454175/1‑93 (MQ=255)
                                                                           ccGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAg              <  2:7657692/93‑1 (MQ=255)
                                                                            cGATTTCGACCATCACCTGGCGCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGc             <  2:4066917/93‑1 (MQ=255)
                                                                            cGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGc             <  1:2828168/93‑1 (MQ=255)
                                                                            cGATTTCGACCATCACCTGGCCATTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGc             >  1:5493971/1‑93 (MQ=255)
                                                                             gATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCg            >  1:6356034/1‑93 (MQ=255)
                                                                               tttCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGtt          <  1:5139868/93‑1 (MQ=255)
                                                                               tttCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGtt          <  2:4678398/93‑1 (MQ=255)
                                                                                 tCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTa        <  1:6515876/93‑1 (MQ=255)
                                                                                 tCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTa        >  1:3457544/1‑93 (MQ=255)
                                                                                  cGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTAc       <  1:5278327/93‑1 (MQ=255)
                                                                                  cGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTAc       <  1:2598082/93‑1 (MQ=255)
                                                                                  cGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTAc       <  2:1620134/93‑1 (MQ=255)
                                                                                   gACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACg      >  1:2419304/1‑93 (MQ=255)
                                                                                    gccATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACgg     <  1:1378333/92‑1 (MQ=255)
                                                                                    aCCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACgg     >  2:7530343/1‑93 (MQ=255)
                                                                                    aCCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACgg     >  2:3385556/1‑93 (MQ=255)
                                                                                      cATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGc   <  2:2331319/93‑1 (MQ=255)
                                                                                       aTCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCt  <  2:7230255/93‑1 (MQ=255)
                                                                                          |                                                                                         
TCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCT  >  NC_000913/52222‑52401

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