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,9450AG99.6% 788.0 / ‑5.1 253T194A (ACG→GCG) phoBpositive response regulator for pho regulon
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (137/115);  total (138/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

Only 100 of 261 total aligned reads displayed.

ACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGATAT  >  NC_000913/416856‑417034
                                                                                         |                                                                                         
aCTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCgg                                                                                        >  1:1994519/1‑93 (MQ=255)
aCTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCgg                                                                                        <  2:4160921/93‑1 (MQ=255)
  ttttttATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTc                                                                                      >  1:7237023/1‑93 (MQ=255)
   tttttATGACGCATCCAGAGCGAGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGACCTAACGTGTATGTGGAAGACCGCGCGGTCg                                                                                     <  2:7253394/93‑1 (MQ=255)
       tATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGt                                                                                 <  1:812308/93‑1 (MQ=255)
          gACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTccc                                                                              >  1:3547397/1‑92 (MQ=255)
            cGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCcaca                                                                            <  1:3149101/93‑1 (MQ=255)
               aTCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATtc                                                                         <  1:3793877/93‑1 (MQ=255)
               aTCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATtc                                                                         >  1:2039494/1‑93 (MQ=255)
                tCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGCACTAACGTGTATGTGGAAGACCGGGCGGTCGATGTCCACATtcg                                                                        >  2:2377206/1‑93 (MQ=255)
                  cTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATtcgtc                                                                      <  2:1427771/93‑1 (MQ=255)
                  cTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATtcgtc                                                                      >  2:2861511/1‑93 (MQ=255)
                  cTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATtcgtc                                                                      >  1:7691100/1‑93 (MQ=255)
                  cTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATtcgtc                                                                      <  2:126200/93‑1 (MQ=255)
                   tGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATtcgtcc                                                                     >  2:6312255/1‑92 (MQ=255)
                     aGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGcc                                                                   <  1:1768336/93‑1 (MQ=255)
                      gcgcGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCgccgccc                                                                  >  2:6124565/1‑92 (MQ=255)
                      gcgcGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCt                                                                  >  1:4769580/1‑93 (MQ=255)
                        gcgTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGc                                                                >  1:4604029/1‑93 (MQ=255)
                        gcgTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGc                                                                <  2:1044208/93‑1 (MQ=255)
                        gcgTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGc                                                                >  2:337900/1‑93 (MQ=255)
                         cgTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGGGCGGTCGATGTCCCCATTCGTCGCCTGCg                                                               >  1:1962064/1‑93 (MQ=255)
                          gtgtACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGGACTAACGTGTATGTGGAAGAGCGCGCGGTCGATGTCCACATTCGTCGGCTGCGt                                                              >  2:644009/1‑93 (MQ=255)
                           tgtACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTa                                                             >  1:5322559/1‑93 (MQ=255)
                           tgtACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTa                                                             <  2:3180545/93‑1 (MQ=255)
                             tACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTaaa                                                           <  1:4180577/93‑1 (MQ=255)
                              tcAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAg                                                          <  1:208160/92‑1 (MQ=255)
                               cAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGc                                                         >  2:1241320/1‑93 (MQ=255)
                               cAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGc                                                         >  2:5824970/1‑93 (MQ=255)
                               cAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGc                                                         >  2:6596775/1‑93 (MQ=255)
                                 gCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCAc                                                       >  1:7936191/1‑93 (MQ=255)
                                 gCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCAc                                                       <  2:5029067/93‑1 (MQ=255)
                                 gCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCAc                                                       >  1:6863904/1‑93 (MQ=255)
                                 gCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCAc                                                       <  2:7432324/93‑1 (MQ=255)
                                  ccGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACt                                                      <  2:3315148/93‑1 (MQ=255)
                                   cgcgAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTg                                                     <  1:2293862/93‑1 (MQ=255)
                                      gAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAg                                                  <  2:5628773/93‑1 (MQ=255)
                                       agcagcTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGc                                                 >  2:1922638/1‑93 (MQ=255)
                                       agcagcTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGc                                                 <  2:4660497/93‑1 (MQ=255)
                                       agcagcTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGc                                                 >  2:4324804/1‑93 (MQ=255)
                                       agcagcTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGc                                                 <  1:3090429/93‑1 (MQ=255)
                                       agcagcTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGc                                                 <  2:7480187/93‑1 (MQ=255)
                                        gcagcTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGcc                                                >  1:4537129/1‑93 (MQ=255)
                                        gcagcTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGcc                                                >  1:4276052/1‑93 (MQ=255)
                                            cTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCcggc                                            <  2:7796262/93‑1 (MQ=255)
                                            cTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCcggc                                            >  2:395234/1‑93 (MQ=255)
                                             tGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCcggcg                                           >  1:2803044/1‑93 (MQ=255)
                                              gTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCcggcgg                                          >  1:5576554/1‑93 (MQ=255)
                                              gTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCcggcgg                                          <  2:1759802/93‑1 (MQ=255)
                                               ttAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCggg                                         >  2:7087928/1‑93 (MQ=255)
                                                 aaaCCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCa                                       <  2:624553/93‑1 (MQ=255)
                                                 aaaCCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCa                                       <  2:5446698/93‑1 (MQ=255)
                                                 aaaCCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCa                                       >  1:6129066/1‑93 (MQ=255)
                                                  aaCCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCAt                                      <  1:4335906/93‑1 (MQ=255)
                                                      aCGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGAcc                                  >  1:3304030/1‑93 (MQ=255)
                                                      aCGTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGAcc                                  >  2:4838625/1‑93 (MQ=255)
                                                        gTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGc                                >  1:1137504/1‑93 (MQ=255)
                                                        gTCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGc                                >  1:5427106/1‑93 (MQ=255)
                                                         tCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCa                               >  2:4499111/1‑93 (MQ=255)
                                                         tCTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCa                               >  2:6915780/1‑93 (MQ=255)
                                                          cTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCAt                              <  1:5851703/93‑1 (MQ=255)
                                                          cTGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCAt                              <  2:6950318/93‑1 (MQ=255)
                                                           tGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATg                             <  2:1701041/93‑1 (MQ=255)
                                                           tGGGGAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATg                             <  1:950258/93‑1 (MQ=255)
                                                            ggggAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACTGCATgg                            <  1:3822187/93‑1 (MQ=255)
                                                            ggggAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATgg                            <  2:4086460/93‑1 (MQ=255)
                                                             gggAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATggg                           >  2:5327953/1‑92 (MQ=255)
                                                              ggAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTg                          >  1:2622991/1‑93 (MQ=255)
                                                               gAACTAACGTGTATGTGGAAGACTGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGc                         <  2:4298830/93‑1 (MQ=255)
                                                               gAACTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGGGc                         >  1:1155686/1‑93 (MQ=255)
                                                                aaCTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCa                        >  1:5668180/1‑93 (MQ=255)
                                                                aaCTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCa                        <  1:600022/93‑1 (MQ=255)
                                                                 aCTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAg                       >  2:6560235/1‑93 (MQ=255)
                                                                  cTAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGa                      <  1:6326254/93‑1 (MQ=255)
                                                                   tAACGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGAc                     <  1:343558/93‑1 (MQ=255)
                                                                     aCGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCg                   <  1:3981562/93‑1 (MQ=255)
                                                                      cGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGt                  >  1:5631441/1‑93 (MQ=255)
                                                                      cGTGTATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGt                  <  2:2306188/93‑1 (MQ=255)
                                                                       gtgtATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTg                 <  2:2092423/93‑1 (MQ=255)
                                                                        tgtATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTgc                <  2:3612057/93‑1 (MQ=255)
                                                                        tgtATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTgc                >  1:912018/1‑93 (MQ=255)
                                                                        tgtATGTGGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTgc                >  2:28020/1‑93 (MQ=255)
                                                                            tgtgGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCgg            >  2:7569677/1‑93 (MQ=255)
                                                                             gtgGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGt           >  1:7412002/1‑93 (MQ=255)
                                                                             gtgGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGt           <  2:4649930/93‑1 (MQ=255)
                                                                              tgGAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTa          >  1:1338987/1‑93 (MQ=255)
                                                                               ggAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTAc         <  1:6014677/93‑1 (MQ=255)
                                                                               ggAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTAc         >  1:7465115/1‑93 (MQ=255)
                                                                               ggAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTAc         >  1:2889449/1‑93 (MQ=255)
                                                                               ggAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTAc         <  1:2730174/93‑1 (MQ=255)
                                                                                gAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACa        <  2:5528315/93‑1 (MQ=255)
                                                                                gAAGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACa        <  2:3766569/93‑1 (MQ=255)
                                                                                 aaGACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAg       <  1:5895126/93‑1 (MQ=255)
                                                                                   gACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGa     >  1:7063541/1‑93 (MQ=255)
                                                                                   gACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGa     >  2:7726218/1‑93 (MQ=255)
                                                                                   gACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGa     >  1:1599008/1‑93 (MQ=255)
                                                                                   gACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGa     <  2:7237022/93‑1 (MQ=255)
                                                                                   gACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGa     >  2:2476233/1‑93 (MQ=255)
                                                                                   gACCGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGa     <  1:6443367/93‑1 (MQ=255)
                                                                                     ccGCGCGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGata   <  2:793424/93‑1 (MQ=255)
                                                                                         |                                                                                         
ACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGATAT  >  NC_000913/416856‑417034

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: