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.7% 1120.1 / NA 334A40V (GCG→GTG) ksgAS‑adenosylmethionine‑6‑N',N'‑adenosyl (rRNA) dimethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (121/212);  minor base C (0/1);  total (121/213)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 346 total aligned reads displayed.

TCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAG  >  NC_000913/52222‑52404
                                                                                          |                                                                                            
tCGCGGTCAAGTTCGATGGCCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCAcc                                                                                            <  1:3200712/93‑1 (MQ=255)
 cgcgGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCt                                                                                           <  2:4826518/93‑1 (MQ=255)
  gcgGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTg                                                                                          <  1:3045864/93‑1 (MQ=255)
   cgGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTgg                                                                                         >  2:9703684/1‑93 (MQ=255)
      tCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGccc                                                                                      >  1:5644830/1‑93 (MQ=255)
      tCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGccc                                                                                      >  1:4464690/1‑93 (MQ=255)
        aaGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCtt                                                                                    <  1:2246546/93‑1 (MQ=255)
        aaGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCtt                                                                                    <  1:4342653/93‑1 (MQ=255)
         aGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCttt                                                                                   >  1:6458083/1‑93 (MQ=255)
         aGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCttt                                                                                   <  2:7528961/93‑1 (MQ=255)
           ttCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTg                                                                                 <  2:1921982/93‑1 (MQ=255)
            tCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGc                                                                                >  1:5941535/1‑93 (MQ=255)
            tCGATGACCGTCAGCTGGCCCAGCCGTCCGCCGACCGGTTCGGTCAATGCCCCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGc                                                                                <  1:5207926/93‑1 (MQ=255)
             cGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCg                                                                               <  1:2666700/93‑1 (MQ=255)
              gctGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCga                                                                              <  2:6163561/91‑2 (MQ=255)
              gATGGCCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCgg                                                                              <  2:3286512/93‑1 (MQ=255)
              gATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCgg                                                                              >  1:6938342/1‑93 (MQ=255)
               aTGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCggg                                                                             >  1:9985176/1‑93 (MQ=255)
                tGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGt                                                                            >  1:8055326/1‑93 (MQ=255)
                 gACGGTCAGCTGGTCCAGACGTGCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGGCCATCACCTGGCCCTTTTGCGGGtt                                                                           <  2:951458/93‑1 (MQ=255)
                  gccGTCGGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGACGATTTCGACCATAACCTGGCCCTTTTGCGGGTTa                                                                          <  1:1793839/92‑1 (MQ=255)
                   ccGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTaa                                                                         <  1:10022298/93‑1 (MQ=255)
                   ccGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTaa                                                                         <  1:806571/93‑1 (MQ=255)
                   ccGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTaa                                                                         >  2:8824270/1‑93 (MQ=255)
                   ccGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTaa                                                                         <  1:2371743/93‑1 (MQ=255)
                   ccGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTaa                                                                         >  2:924300/1‑93 (MQ=255)
                   ccGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTaa                                                                         >  2:3549878/1‑93 (MQ=255)
                   ccGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTaa                                                                         <  1:5757547/93‑1 (MQ=255)
                   ccGGCAGCTGGCCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTaa                                                                         <  2:1552747/93‑1 (MQ=255)
                     gTCAGCTGGTCCAGACGTTCGCCGACAGGTTCGGTCAATGCCGCCAGACCGGGGCCGATGTCGCCCATCACCTGGCCCTTTTGCGGGTTAATg                                                                       <  1:350245/93‑1 (MQ=255)
                       cAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGc                                                                     <  1:908001/93‑1 (MQ=255)
                       cAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGc                                                                     <  2:1301458/93‑1 (MQ=255)
                       cAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGc                                                                     <  1:8560188/93‑1 (MQ=255)
                         gCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAg                                                                   >  2:6639532/1‑93 (MQ=255)
                         gCTGGGCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAg                                                                   <  1:6851674/93‑1 (MQ=255)
                           tGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGac                                                                 <  2:4985922/93‑1 (MQ=255)
                            ggTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGaca                                                                <  1:9592817/93‑1 (MQ=255)
                            ggTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGACCTTTTGCGGGTTAATGGCAGaaa                                                                >  2:4175863/1‑91 (MQ=255)
                             gTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGacac                                                               <  1:9354336/93‑1 (MQ=255)
                                 aGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATa                                                           <  2:6526693/93‑1 (MQ=255)
                                   aCGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGATCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACt                                                         >  2:4042693/1‑93 (MQ=255)
                                   aCGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACt                                                         <  1:4040203/93‑1 (MQ=255)
                                     gatCGCCGACCGGGTCGGTCAAGGCCGCGAGACAGGGGCCGAGTTCGACCATCACCTGGCCCGTTTGCGGGTTAATGGCAGACACAATACTGt                                                       <  1:9527442/91‑1 (MQ=255)
                                     gTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCTATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGt                                                       <  1:4851590/93‑1 (MQ=255)
                                      ttCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTc                                                      <  1:3663542/93‑1 (MQ=255)
                                        cGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCCCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGa                                                    <  1:2774369/93‑1 (MQ=255)
                                         gCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGAt                                                   >  2:6330794/1‑93 (MQ=255)
                                         cccGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGAt                                                   <  2:2254628/92‑1 (MQ=255)
                                           cGACCGGTTCGGTCACTGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTAAATGGCAGACACAATACTGTCGATCa                                                 <  2:7685251/93‑1 (MQ=255)
                                            gACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCAc                                                <  2:2996154/93‑1 (MQ=255)
                                             aCCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACg                                               <  2:730695/93‑1 (MQ=255)
                                                ggTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAAc                                            <  1:7757525/93‑1 (MQ=255)
                                                 gTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACt                                           >  1:370863/1‑93 (MQ=255)
                                                   tCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGa                                         <  2:6587008/93‑1 (MQ=255)
                                                   tCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGa                                         <  1:745758/93‑1 (MQ=255)
                                                   tCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACCGTCGATCACGAACTGa                                         >  1:8983942/1‑93 (MQ=255)
                                                    cGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGAt                                        >  2:6761584/1‑93 (MQ=255)
                                                    cGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGAt                                        <  2:4652244/93‑1 (MQ=255)
                                                    cGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGAt                                        <  2:5431257/93‑1 (MQ=255)
                                                     ggTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATc                                       <  2:3494044/93‑1 (MQ=255)
                                                       tCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGt                                     >  2:5171059/1‑93 (MQ=255)
                                                       tCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGt                                     <  1:3326602/93‑1 (MQ=255)
                                                        cAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGtt                                    >  2:3673389/1‑93 (MQ=255)
                                                         aaTGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTg                                   >  1:3801870/1‑93 (MQ=255)
                                                         aaTGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTg                                   <  2:4338960/93‑1 (MQ=255)
                                                          aTGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTga                                  >  2:7974298/1‑93 (MQ=255)
                                                          aTGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGGTga                                  >  2:71020/1‑93 (MQ=255)
                                                           tGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTgag                                 >  1:7237890/1‑93 (MQ=255)
                                                            gccgccAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAGCTGATCGTTgaga                                <  2:3927681/93‑1 (MQ=255)
                                                                  aGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGttt                          >  1:1309539/1‑93 (MQ=255)
                                                                   gACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAGCTGATCGTTGAGAAAGtttg                         >  2:8117843/1‑92 (MQ=255)
                                                                   gACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGtttt                         >  1:6290850/1‑93 (MQ=255)
                                                                     ccGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGc                       >  1:547607/1‑93 (MQ=255)
                                                                      cGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGcc                      <  1:9905309/93‑1 (MQ=255)
                                                                       ggggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGccc                     >  2:7408791/1‑93 (MQ=255)
                                                                       ggggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGccc                     >  1:6467094/1‑93 (MQ=255)
                                                                        gggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCg                    >  1:1919038/1‑93 (MQ=255)
                                                                        gggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCg                    >  2:5957661/1‑93 (MQ=255)
                                                                        gggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCg                    <  2:1096868/93‑1 (MQ=255)
                                                                         ggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGa                   <  2:2487974/93‑1 (MQ=255)
                                                                          gCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGaa                  <  1:3969892/93‑1 (MQ=255)
                                                                          gCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGaa                  <  1:6089492/93‑1 (MQ=255)
                                                                          gCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGaa                  >  2:9400112/1‑93 (MQ=255)
                                                                           ccGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAg                 <  2:8157888/93‑1 (MQ=255)
                                                                           ccGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAg                 <  2:5650302/93‑1 (MQ=255)
                                                                            cGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGc                <  1:1505087/93‑1 (MQ=255)
                                                                              aTTGCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGt              <  1:8330181/93‑1 (MQ=255)
                                                                               tttCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGtt             >  2:1775379/1‑93 (MQ=255)
                                                                               tttCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGtt             <  2:9116806/93‑1 (MQ=255)
                                                                                ttCGACCAGCACCTGGCCCTTGTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGttt            <  1:5334051/93‑1 (MQ=255)
                                                                                 tCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTa           <  2:2067265/93‑1 (MQ=255)
                                                                                 tCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTa           <  2:6297917/93‑1 (MQ=255)
                                                                                  cGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTAc          <  2:2600489/93‑1 (MQ=255)
                                                                                   gACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACg         <  1:387722/93‑1 (MQ=255)
                                                                                     ccATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACggg       >  1:2503071/1‑93 (MQ=255)
                                                                                      cATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGc      >  1:8788988/1‑93 (MQ=255)
                                                                                      cATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGc      <  1:1093392/93‑1 (MQ=255)
                                                                                       aTCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCt     >  1:3431032/1‑93 (MQ=255)
                                                                                       aTCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCt     >  2:7796513/1‑93 (MQ=255)
                                                                                          gCCTGGCCCTTTTGCGGGTTAATGGCAGACCCAACACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAg  <  2:5073317/93‑1 (MQ=255)
                                                                                          |                                                                                            
TCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAG  >  NC_000913/52222‑52404

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: