Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 416,945 A→G T194A (ACG→GCG)  phoB → positive response regulator for pho regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913416,9450AG97.6% 385.9 / ‑2.4 125T194A (ACG→GCG) phoBpositive response regulator for pho regulon
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base G (66/56);  total (66/59)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 133 total aligned reads displayed.

TGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGATA  >  NC_000913/416853‑417033
                                                                                            |                                                                                        
tgcACTTTTTCAGGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCa                                                                                          <  2:4407806/93‑1 (MQ=255)
  cACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACcgcgcg                                                                                        >  2:3078965/1‑93 (MQ=255)
  cACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTAAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACcgcgcg                                                                                        >  1:506346/1‑93 (MQ=255)
   aCTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCgg                                                                                       >  1:1884111/1‑93 (MQ=255)
    cTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGt                                                                                      <  2:3832356/93‑1 (MQ=255)
     ttttttATGACGCATCCTGAGCGCGTGTAGAGCCGCGAGAAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTc                                                                                     <  1:3559074/93‑1 (MQ=255)
      tttttATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCg                                                                                    <  1:4575774/93‑1 (MQ=255)
      tttttATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCg                                                                                    <  1:968870/93‑1 (MQ=255)
       ttttATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGa                                                                                   <  1:2519093/93‑1 (MQ=255)
       ttttATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGa                                                                                   <  1:633723/93‑1 (MQ=255)
        tttATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGAt                                                                                  <  1:3004433/93‑1 (MQ=255)
        tttATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGAt                                                                                  <  1:4691329/93‑1 (MQ=255)
        tttATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGAt                                                                                  >  1:1915966/1‑93 (MQ=255)
         ttATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATg                                                                                 <  2:91355/93‑1 (MQ=255)
         ttATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATg                                                                                 >  1:496320/1‑93 (MQ=255)
          tATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGt                                                                                >  2:4092905/1‑93 (MQ=255)
          tATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGt                                                                                <  2:2886968/93‑1 (MQ=255)
           aTGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTc                                                                               >  2:2629252/1‑93 (MQ=255)
               cGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCcaca                                                                           <  1:1808760/93‑1 (MQ=255)
                gCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCCCGCGGTCGATGTCCACAt                                                                          >  2:766807/1‑93 (MQ=255)
                 cATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACAtt                                                                         >  2:4156486/1‑93 (MQ=255)
                 cATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACAtt                                                                         >  1:2426412/1‑93 (MQ=255)
                 cATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACAtt                                                                         <  1:2859418/93‑1 (MQ=255)
                 cATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACAtt                                                                         <  2:4399153/93‑1 (MQ=255)
                  aTCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATtc                                                                        >  2:1909575/1‑93 (MQ=255)
                   tCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATtcg                                                                       >  1:2156844/1‑93 (MQ=255)
                   tCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATtcg                                                                       >  1:1638816/1‑93 (MQ=255)
                    ccTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATtcgt                                                                      <  1:3335489/93‑1 (MQ=255)
                     cTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATtcgtc                                                                     >  1:2892443/1‑93 (MQ=255)
                      tGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATtcgtcg                                                                    >  1:128086/1‑93 (MQ=255)
                      tGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATtcgtcg                                                                    >  1:2129521/1‑93 (MQ=255)
                       gAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGc                                                                   >  1:549682/1‑93 (MQ=255)
                        aGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGcc                                                                  >  1:3826940/1‑93 (MQ=255)
                        aGCGCGTGTAAAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAAGGTGTATGTGGCAGAACGCGCGGTCGATGTACACATTCGTCAcc                                                                  >  1:3826939/1‑93 (MQ=255)
                         gcgcGTGTACAGCCCCGAGCAGCTGTTAAACCACGTCTGGGGGACTTACGTGTATGTGGAAGACCGCGCGGTCGATGTCCCCATTCGTCGcac                                                                 >  2:4703709/1‑91 (MQ=255)
                            cgTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCg                                                              >  2:3252839/1‑93 (MQ=255)
                             gtgtACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGt                                                             >  1:459929/1‑93 (MQ=255)
                                  cAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGc                                                        <  2:3542409/93‑1 (MQ=255)
                                     ccGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCGTGCGTAAAGCACt                                                     >  1:3379218/1‑93 (MQ=255)
                                     ccGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACt                                                     <  1:2804881/93‑1 (MQ=255)
                                     ccGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACt                                                     >  1:3379217/1‑93 (MQ=255)
                                     ccGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACt                                                     <  1:229114/93‑1 (MQ=255)
                                      cgcgAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTg                                                    <  2:4573877/93‑1 (MQ=255)
                                      cgcgAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTg                                                    <  2:451580/93‑1 (MQ=255)
                                       gcgAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTgg                                                   >  1:4513501/1‑93 (MQ=255)
                                       gcgAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTgg                                                   >  2:3448925/1‑93 (MQ=255)
                                         gAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAg                                                 >  1:3527603/1‑93 (MQ=255)
                                         gAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAg                                                 >  1:194778/1‑93 (MQ=255)
                                          agcagcTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGc                                                >  2:3510391/1‑93 (MQ=255)
                                           gcagcTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGcc                                               <  1:1210051/93‑1 (MQ=255)
                                           gcagcTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGcc                                               <  2:269673/93‑1 (MQ=255)
                                            cagcTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGccc                                              <  1:1522353/93‑1 (MQ=255)
                                             agcTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCcg                                             <  2:3466210/93‑1 (MQ=255)
                                              gcTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCcgg                                            >  2:83637/1‑93 (MQ=255)
                                              gcTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCcgg                                            >  2:1658824/1‑93 (MQ=255)
                                               gtGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCcggc                                           <  1:3527148/92‑1 (MQ=255)
                                                tGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCcggcg                                          >  2:607818/1‑93 (MQ=255)
                                                tGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCcggcg                                          >  1:3931435/1‑93 (MQ=255)
                                                    aaaCCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCa                                      <  2:2830111/93‑1 (MQ=255)
                                                     aaCCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCAt                                     >  2:1756036/1‑93 (MQ=255)
                                                        cACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGAc                                  <  1:3050180/93‑1 (MQ=255)
                                                            tCTGGTGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCa                              <  1:1115594/93‑1 (MQ=255)
                                                              tGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATg                            <  1:790168/93‑1 (MQ=255)
                                                              tGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATg                            >  2:3025330/1‑93 (MQ=255)
                                                              tGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATg                            <  1:3388998/93‑1 (MQ=255)
                                                              tGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATg                            >  1:1564674/1‑93 (MQ=255)
                                                              tGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTAGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATg                            >  2:918455/1‑93 (MQ=255)
                                                               ggggAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATgg                           >  1:1994496/1‑93 (MQ=255)
                                                                gggAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGt                          <  2:1354498/93‑1 (MQ=255)
                                                                  gAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGc                        <  1:154551/93‑1 (MQ=255)
                                                                   aaCTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCa                       >  2:2267330/1‑93 (MQ=255)
                                                                   aaCTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCa                       >  1:1245140/1‑93 (MQ=255)
                                                                    aCTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAg                      <  1:811658/93‑1 (MQ=255)
                                                                     cTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGa                     <  2:2287194/93‑1 (MQ=255)
                                                                      tAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGAc                    >  2:749173/1‑93 (MQ=255)
                                                                       aaCGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGagc                   >  1:1040536/1‑91 (MQ=255)
                                                                       aaCGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGAcc                   >  1:1040384/1‑93 (MQ=255)
                                                                       aaCGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGAcc                   <  1:3199214/93‑1 (MQ=255)
                                                                        aCGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCg                  >  2:2425640/1‑93 (MQ=255)
                                                                        aCGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCg                  <  1:2258046/93‑1 (MQ=255)
                                                                         cGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGt                 <  1:3300574/93‑1 (MQ=255)
                                                                           tgtATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTgc               <  1:962279/93‑1 (MQ=255)
                                                                            gtATGTGGAAGACCGGGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTgcg              >  2:4057782/1‑93 (MQ=255)
                                                                             tATGTGGAAGACCGCGCGGTCGATGTCCACATTAGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTgcgc             <  1:918455/93‑1 (MQ=255)
                                                                              aTGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGCAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGcgcg            <  2:2142408/93‑1 (MQ=255)
                                                                               tgtgGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCgg           >  1:1074029/1‑93 (MQ=255)
                                                                                gtgGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGt          >  2:133245/1‑93 (MQ=255)
                                                                                gtgGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGt          <  1:2967352/93‑1 (MQ=255)
                                                                                 tgGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTa         >  2:4353575/1‑93 (MQ=255)
                                                                                 tgGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTa         >  1:3013616/1‑93 (MQ=255)
                                                                                  ggAAGGCCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTAc        <  2:3588973/93‑1 (MQ=255)
                                                                                  ggAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTAc        >  2:3820172/1‑93 (MQ=255)
                                                                                  ggAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTAc        >  1:2275762/1‑93 (MQ=255)
                                                                                  agAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTAc        <  1:591389/92‑1 (MQ=255)
                                                                                   gAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACa       <  1:770600/93‑1 (MQ=255)
                                                                                    aaGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAg      <  2:1556788/93‑1 (MQ=255)
                                                                                     aGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAgg     <  2:1537922/93‑1 (MQ=255)
                                                                                     aGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAgg     >  2:1932675/1‑93 (MQ=255)
                                                                                        ccGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGata  >  2:3725006/1‑93 (MQ=255)
                                                                                        ccGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGata  >  1:1764963/1‑93 (MQ=255)
                                                                                            |                                                                                        
TGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGATA  >  NC_000913/416853‑417033

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