Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,961,966 A→G intergenic (‑122/+14) gppA ← / ← rhlB guanosine pentaphosphatase/putative ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,961,9660AG99.6% 872.9 / NA 264intergenic (‑122/+14)gppA/rhlBguanosine pentaphosphatase/putative ATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (124/139);  total (124/140)
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 268 total aligned reads displayed.

GCTGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTATGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGGCGC  >  NC_000913/3961876‑3962055
                                                                                          |                                                                                         
gCTGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGtgtg                                                                                         <  1:3099371/93‑1 (MQ=255)
  tGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGAc                                                                                       >  1:6580198/1‑93 (MQ=255)
  tGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGAc                                                                                       >  1:1173904/1‑93 (MQ=255)
  tGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGAc                                                                                       <  2:512113/93‑1 (MQ=255)
    cATTCTTTCACCTGAATCGCTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAt                                                                                     <  1:5485279/93‑1 (MQ=255)
    cATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAt                                                                                     <  2:5901019/93‑1 (MQ=255)
      ttCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAttt                                                                                   <  1:2059454/93‑1 (MQ=255)
      ttCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAttt                                                                                   >  1:1464062/1‑93 (MQ=255)
       tCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAtttt                                                                                  <  2:4165345/93‑1 (MQ=255)
        cTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAttttt                                                                                 >  1:2565129/1‑93 (MQ=255)
         tttCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAtttttt                                                                                >  2:2199435/1‑93 (MQ=255)
         tttCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAtttttt                                                                                <  1:932167/93‑1 (MQ=255)
          ttCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTttta                                                                               <  1:6535929/93‑1 (MQ=255)
           tCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttat                                                                              <  2:2353325/93‑1 (MQ=255)
           tCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttat                                                                              <  1:4047724/93‑1 (MQ=255)
           tCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttat                                                                              >  1:1823301/1‑93 (MQ=255)
            cACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttatt                                                                             <  1:4517484/93‑1 (MQ=255)
             aCCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttattt                                                                            <  1:20106/93‑1 (MQ=255)
             aCCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttattt                                                                            >  2:1075545/1‑93 (MQ=255)
               cTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTaa                                                                          >  2:5383285/1‑93 (MQ=255)
                tGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAAc                                                                         >  2:5629547/1‑93 (MQ=255)
                  aaTCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCt                                                                       >  2:3815503/1‑93 (MQ=255)
                  aaTCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCt                                                                       >  1:3607693/1‑93 (MQ=255)
                  aaTCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCt                                                                       <  1:6395626/93‑1 (MQ=255)
                     cACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGaa                                                                    <  1:5354143/93‑1 (MQ=255)
                     cACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGaa                                                                    <  2:4655146/93‑1 (MQ=255)
                     cACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGaa                                                                    <  1:231425/93‑1 (MQ=255)
                     cACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGaa                                                                    <  2:6490433/93‑1 (MQ=255)
                      aCTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAac                                                                   >  1:3530883/1‑93 (MQ=255)
                      aCTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAac                                                                   <  2:4528242/93‑1 (MQ=255)
                       cTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacg                                                                  >  2:1446588/1‑93 (MQ=255)
                        ttACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacga                                                                 >  2:3390493/1‑93 (MQ=255)
                        ttACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacga                                                                 >  1:6242790/1‑93 (MQ=255)
                           ctacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgacga                                                              >  2:3310239/1‑93 (MQ=255)
                             acAGTAAGCTCCTCGGGATTCATTCTCTTGCCGTCTTGGTACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCGGAacggcgacg                                                            >  1:6731611/1‑93 (MQ=255)
                             acAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgacgacg                                                            <  1:1008288/93‑1 (MQ=255)
                              cAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAAcgacgacga                                                           >  2:1981267/1‑93 (MQ=255)
                               aGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGAt                                                          <  1:3695540/93‑1 (MQ=255)
                               aGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGAt                                                          <  1:6487226/93‑1 (MQ=255)
                                gTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGAtt                                                         >  1:5022738/1‑93 (MQ=255)
                                gTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGAtt                                                         <  2:4062891/93‑1 (MQ=255)
                                 tAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTa                                                        <  2:968620/93‑1 (MQ=255)
                                 tAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGGGACATTTTTTATTTAACCTGAACGACGACGATTa                                                        <  1:2719358/93‑1 (MQ=255)
                                  taGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAGTGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTAc                                                       <  2:6734670/92‑1 (MQ=255)
                                  aaNCTCATCGGGATTCATTCTNTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTAc                                                       <  2:2874159/93‑1 (MQ=255)
                                   aGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTAcg                                                      >  2:6719127/1‑93 (MQ=255)
                                    gCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGACCGACGACGATTAcgc                                                     <  1:6022891/93‑1 (MQ=255)
                                     cTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTAcgcg                                                    >  1:2243005/1‑93 (MQ=255)
                                     cTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTAcgcg                                                    <  2:3455312/93‑1 (MQ=255)
                                      tCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCgg                                                   <  1:6421412/93‑1 (MQ=255)
                                       cATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGa                                                  <  2:4677960/93‑1 (MQ=255)
                                       cATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGa                                                  <  2:495900/93‑1 (MQ=255)
                                         tCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAgc                                                >  2:1040049/1‑93 (MQ=255)
                                             gATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCa                                            >  1:3357374/1‑93 (MQ=255)
                                                tcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTa                                         <  2:4419500/93‑1 (MQ=255)
                                                 cattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTac                                        <  1:2372968/93‑1 (MQ=255)
                                                  attcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacg                                       >  1:4301461/1‑93 (MQ=255)
                                                   ttcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacga                                      <  2:3861033/93‑1 (MQ=255)
                                                   ttcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacga                                      >  2:3110824/1‑93 (MQ=255)
                                                    tctcTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacgac                                     <  2:4262886/93‑1 (MQ=255)
                                                    tctcTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacgac                                     <  2:780582/93‑1 (MQ=255)
                                                    tctcTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacgac                                     <  2:6509070/93‑1 (MQ=255)
                                                      tctTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGAcgc                                   >  1:3022532/1‑93 (MQ=255)
                                                          gCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGAc                               >  2:1713609/1‑93 (MQ=255)
                                                          gCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGAc                               <  1:723728/93‑1 (MQ=255)
                                                          gCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGAc                               <  1:4160565/93‑1 (MQ=255)
                                                            cATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCa                             >  1:2115584/1‑93 (MQ=255)
                                                              tCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCAtt                           >  1:2029757/1‑93 (MQ=255)
                                                              tCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCAtt                           >  1:1660528/1‑93 (MQ=255)
                                                              tCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCAtt                           >  2:6288360/1‑93 (MQ=255)
                                                                ttGATACAGTTTGAATGATTTTGAGTATGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGc                         <  2:6079036/93‑1 (MQ=255)
                                                                  gATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCt                       >  2:6194324/1‑93 (MQ=255)
                                                                      cAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTgc                   <  1:5462717/93‑1 (MQ=255)
                                                                       aGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTgcg                  <  2:1515174/93‑1 (MQ=255)
                                                                       aGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTgcg                  >  2:2803365/1‑93 (MQ=255)
                                                                       aGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTgcg                  >  2:4341845/1‑93 (MQ=255)
                                                                        gTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTgcgc                 >  1:1048202/1‑93 (MQ=255)
                                                                        gTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTgcgc                 >  1:5953740/1‑93 (MQ=255)
                                                                        gTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTgcgc                 >  2:1167342/1‑93 (MQ=255)
                                                                         tttGAATGATTTTGCGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGcgcg                >  2:5190816/1‑93 (MQ=255)
                                                                         tttGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGcgcg                >  1:5683068/1‑93 (MQ=255)
                                                                         tttGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGcgcg                >  2:1372565/1‑93 (MQ=255)
                                                                            gAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGgc             >  2:4471665/1‑93 (MQ=255)
                                                                            gAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGgc             >  1:4888917/1‑93 (MQ=255)
                                                                               tGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGcgcg          >  1:663794/1‑93 (MQ=255)
                                                                               tGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGcgcg          <  1:4612654/93‑1 (MQ=255)
                                                                                  ttttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGa       <  1:4594381/93‑1 (MQ=255)
                                                                                  ttttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGa       <  2:1329908/93‑1 (MQ=255)
                                                                                  ttttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGa       <  1:2793587/93‑1 (MQ=255)
                                                                                  ttttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGa       >  1:2145247/1‑93 (MQ=255)
                                                                                  tttgGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGa       <  2:2934293/93‑1 (MQ=255)
                                                                                   tttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAg      <  2:1844666/93‑1 (MQ=255)
                                                                                   tttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAg      <  2:255222/93‑1 (MQ=255)
                                                                                   tttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAg      <  2:2779567/93‑1 (MQ=255)
                                                                                    ttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAgg     <  2:6172351/93‑1 (MQ=255)
                                                                                     tGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGgc    >  1:3799098/1‑93 (MQ=255)
                                                                                      gAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGgcg   >  2:5565943/1‑93 (MQ=255)
                                                                                      gAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGgcg   <  2:2225402/93‑1 (MQ=255)
                                                                                      gAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGgcg   >  2:3692331/1‑93 (MQ=255)
                                                                                       aGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGgcgc  >  1:5132090/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCTGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTATGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGGCGC  >  NC_000913/3961876‑3962055

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