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,3120GA96.4% 438.4 / 1.6 139A40V (GCG→GTG) ksgAS‑adenosylmethionine‑6‑N',N'‑adenosyl (rRNA) dimethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base A (52/80);  total (53/86)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.12e-01

Only 100 of 147 total aligned reads displayed.

TCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCT  >  NC_000913/52222‑52401
                                                                                          |                                                                                         
tCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCAcc                                                                                         <  2:2376462/93‑1 (MQ=255)
  gcgGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTg                                                                                       <  2:1418013/93‑1 (MQ=255)
  gcgGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTg                                                                                       <  1:160110/93‑1 (MQ=255)
   cgGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTgg                                                                                      >  1:3215471/1‑93 (MQ=255)
     gTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGcc                                                                                    >  1:2194877/1‑93 (MQ=255)
      tCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGccc                                                                                   <  1:3439230/93‑1 (MQ=255)
       cAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCt                                                                                  >  2:321268/1‑93 (MQ=255)
           ttCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTg                                                                              >  1:988621/1‑93 (MQ=255)
             cGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCg                                                                            <  1:1942847/93‑1 (MQ=255)
              gATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCgg                                                                           <  1:3303388/93‑1 (MQ=255)
               aTGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCggg                                                                          <  2:3638671/93‑1 (MQ=255)
               aTGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCggg                                                                          <  2:1199326/93‑1 (MQ=255)
               aTGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCggg                                                                          <  2:2692082/93‑1 (MQ=255)
               aTGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCggg                                                                          <  1:1234449/93‑1 (MQ=255)
                tGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGt                                                                         <  2:2304513/93‑1 (MQ=255)
                  aCCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTa                                                                       >  2:916446/1‑93 (MQ=255)
                   ccGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTaa                                                                      >  1:2454203/1‑93 (MQ=255)
                   ccGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTaa                                                                      <  1:3132350/93‑1 (MQ=255)
                    cGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAAt                                                                     <  2:3241070/93‑1 (MQ=255)
                    cGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAAt                                                                     <  1:1679110/93‑1 (MQ=255)
                    cGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAAt                                                                     <  2:760645/93‑1 (MQ=255)
                     gTCGGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATg                                                                    <  2:1939253/93‑1 (MQ=255)
                      tCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATgg                                                                   <  2:3103052/93‑1 (MQ=255)
                       cAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCTGGTTAATGGc                                                                  >  1:1774446/1‑93 (MQ=255)
                       cAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGc                                                                  <  2:448510/93‑1 (MQ=255)
                       cAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGc                                                                  <  2:3451483/93‑1 (MQ=255)
                        aGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCa                                                                 >  1:4477208/1‑93 (MQ=255)
                        aGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCCCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGcg                                                                 >  1:4061158/1‑92 (MQ=255)
                          cTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGa                                                               >  2:4197755/1‑93 (MQ=255)
                            ggTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGaca                                                             <  1:1250458/93‑1 (MQ=255)
                            ggTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGaca                                                             >  1:3123905/1‑93 (MQ=255)
                            ggTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGaca                                                             <  2:1929542/93‑1 (MQ=255)
                              tCCAGACGTTCGCCGACCGGGTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGCCCCTTTTGCGAGTTAATGGCAGAcaca                                                           >  1:24358/1‑93 (MQ=255)
                              gccAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGAcaca                                                           <  1:2846259/92‑1 (MQ=255)
                               ccAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACaa                                                          <  1:4316544/93‑1 (MQ=255)
                               ccAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACaa                                                          <  1:1246571/93‑1 (MQ=255)
                                  gACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATAc                                                       >  1:262201/1‑93 (MQ=255)
                                   aCGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACt                                                      <  2:4676575/93‑1 (MQ=255)
                                    cGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTg                                                     <  2:351389/93‑1 (MQ=255)
                                     gTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGt                                                    <  1:38553/93‑1 (MQ=255)
                                     gTTCGCCGACCGGTTCGGTCAATGACGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGt                                                    <  1:2441753/93‑1 (MQ=255)
                                      ttCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCCCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTc                                                   <  2:2496547/93‑1 (MQ=255)
                                      ttCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTc                                                   <  2:820554/93‑1 (MQ=255)
                                      ttCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTc                                                   <  2:1746839/93‑1 (MQ=255)
                                      ttCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTc                                                   <  1:2031066/93‑1 (MQ=255)
                                        cGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGa                                                 <  1:1158231/93‑1 (MQ=255)
                                         gCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGAt                                                <  2:1707331/93‑1 (MQ=255)
                                         gCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGAt                                                <  1:2412695/93‑1 (MQ=255)
                                         gCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGAt                                                <  1:2367698/93‑1 (MQ=255)
                                         cccGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGAt                                                <  2:514398/92‑1 (MQ=255)
                                          ccGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATc                                               >  1:3740372/1‑93 (MQ=255)
                                           cGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCa                                              >  1:881747/1‑93 (MQ=255)
                                            gACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCAc                                             >  2:1822792/1‑93 (MQ=255)
                                                ggTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAAc                                         >  1:3157101/1‑93 (MQ=255)
                                                ggTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAAc                                         >  2:1076665/1‑93 (MQ=255)
                                                  ttCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTg                                       >  1:1020874/1‑93 (MQ=255)
                                                   tCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGa                                      >  2:2977808/1‑93 (MQ=255)
                                                   tCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGa                                      <  2:933355/93‑1 (MQ=255)
                                                    cGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGAt                                     <  1:4103933/93‑1 (MQ=255)
                                                    cGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGAt                                     >  2:3954083/1‑93 (MQ=255)
                                                     ggTCAATGCCGCAAGACCGGGGCCGATTTCGACCATCACCTGGCCTTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATc                                    <  2:3093147/93‑1 (MQ=255)
                                                      gTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCg                                   >  1:2504387/1‑93 (MQ=255)
                                                        cAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGtt                                 <  2:3608418/93‑1 (MQ=255)
                                                         aaTGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTg                                >  1:2406254/1‑93 (MQ=255)
                                                          aTGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTga                               <  2:1519630/93‑1 (MQ=255)
                                                            gccgccAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTgaga                             <  1:4272823/93‑1 (MQ=255)
                                                               gccAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAg                          <  2:4605803/93‑1 (MQ=255)
                                                                ccAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGt                         <  2:1498390/93‑1 (MQ=255)
                                                                ccAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGt                         <  1:1275539/93‑1 (MQ=255)
                                                                  aGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGttt                       <  1:2025864/93‑1 (MQ=255)
                                                                  aGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGttt                       <  1:2020689/93‑1 (MQ=255)
                                                                  aGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGttt                       <  1:4563213/93‑1 (MQ=255)
                                                                   gACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGtttt                      >  2:2759755/1‑93 (MQ=255)
                                                                     ccGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGc                    >  1:2232194/1‑93 (MQ=255)
                                                                     ccGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGc                    <  1:2174144/93‑1 (MQ=255)
                                                                      cGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGcc                   >  2:4379166/1‑93 (MQ=255)
                                                                      cGGGGCCGATTTCGACCATCACCTGGCCATGTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGcc                   <  2:447804/93‑1 (MQ=255)
                                                                       ggggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGccc                  >  1:3164393/1‑93 (MQ=255)
                                                                        gggCCGATTTCGACCATCACCTGGCCGTTTTGCGGTTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCg                 <  2:2112107/93‑1 (MQ=255)
                                                                        gggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGTGAAAGTTTTGCCCg                 >  2:1033163/1‑93 (MQ=255)
                                                                        gggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCg                 >  2:1826093/1‑93 (MQ=255)
                                                                         ggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGGAATGATCGTTGAGAAAGTTTTGCCCGa                >  1:2998947/1‑93 (MQ=255)
                                                                         ggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGa                >  2:2358522/1‑93 (MQ=255)
                                                                         ggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGa                >  1:1840403/1‑93 (MQ=255)
                                                                          gCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGaa               >  1:1567874/1‑93 (MQ=255)
                                                                          gCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGaa               >  2:557495/1‑93 (MQ=255)
                                                                            tgATTTCGACGATCGCTTGGCCATTTTGCGGGTTAATGGCAGACACAATCCTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGc             <  2:2776800/92‑1 (MQ=255)
                                                                            cGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGc             <  2:1066892/93‑1 (MQ=255)
                                                                            cGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGc             >  2:1454780/1‑93 (MQ=255)
                                                                            cGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGc             <  1:2770960/93‑1 (MQ=255)
                                                                             gATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCg            >  2:3980752/1‑93 (MQ=255)
                                                                             gATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCg            >  1:1919787/1‑93 (MQ=255)
                                                                               tttCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGtt          <  1:909516/93‑1 (MQ=255)
                                                                               tttCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGtt          <  1:33134/93‑1 (MQ=255)
                                                                                 tCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGATCTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTa        <  2:237602/93‑1 (MQ=255)
                                                                                 tCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTa        <  2:2538900/93‑1 (MQ=255)
                                                                                  cGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTAc       <  1:1049562/93‑1 (MQ=255)
                                                                                     ccATCACCTGGCCCTTTTGCGGGTTAATGGCAGacac                                                            <  2:3686056/37‑1 (MQ=255)
                                                                                      cATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGc   >  1:3410163/1‑93 (MQ=255)
                                                                                       aTCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCt  <  1:1657462/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 < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: