Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,151,969 T→A 5.2% R168S (AGA→AGT yegK ← ser/thr phosphatase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,151,9690TA 5.2% 410.7 / 10.6 115R168S (AGA→AGTyegKser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base T (57/52);  new base A (3/3);  total (60/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-02

Only 100 of 118 total aligned reads displayed.

CCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGAT  >  NC_000913/2151821‑2152116
                                                                                                                                                    |                                                                                                                                                   
ccTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGt                                                                                                                                                     <  2:1463360/149‑1 (MQ=255)
  tGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACt                                                                                                                                                   >  2:286010/1‑149 (MQ=255)
   gCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTg                                                                                                                                                  <  4:41266/149‑1 (MQ=255)
          tgcCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActt                                                                                                                                           <  1:286010/149‑3 (MQ=255)
             cgccAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCAt                                                                                                                                        <  2:988262/149‑1 (MQ=255)
                  gTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtcgt                                                                                                                                  >  2:420929/1‑149 (MQ=255)
                    ccATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtcgtc                                                                                                                                 <  1:3521/149‑1 (MQ=255)
                      aTTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCgg                                                                                                                               >  4:344474/1‑149 (MQ=255)
                      aTTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCgg                                                                                                                               >  2:182701/1‑149 (MQ=255)
                      aTTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGGGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCgg                                                                                                                               >  1:200137/1‑149 (MQ=255)
                       ttGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGt                                                                                                                              <  1:182701/149‑1 (MQ=255)
                       ttGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGt                                                                                                                              >  1:1255508/1‑149 (MQ=255)
                         gAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTaa                                                                                                                            <  2:200137/149‑1 (MQ=255)
                            aaaCGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGaa                                                                                                                        <  1:420929/149‑1 (MQ=255)
                              aCGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATc                                                                                                                                       >  4:38662/1‑133 (MQ=255)
                              aCGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATc                                                                                                                                       <  3:38662/133‑1 (MQ=255)
                                 gggTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGc                                                                                                                    <  3:1287151/149‑1 (MQ=255)
                                    tAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCt                                                                                                                                                   <  1:595959/115‑1 (MQ=255)
                                        gaaagGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGt                                                                                                             >  3:886604/1‑149 (MQ=255)
                                        gaaagGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGt                                                                                                             >  4:526834/1‑149 (MQ=255)
                                           agGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTgg                                                                                                          >  3:1371542/1‑149 (MQ=255)
                                             gCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGt                                                                                                                      <  1:984296/135‑1 (MQ=255)
                                             gCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGt                                                                                                                      >  2:984296/1‑135 (MQ=255)
                                               acaTGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCata                                                                                                      >  4:640907/1‑149 (MQ=255)
                                                    aggagAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTg                                                                                                           <  1:1271077/139‑1 (MQ=255)
                                                    aggagAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTg                                                                                                           >  2:1271077/1‑139 (MQ=255)
                                                      gagaATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTgg                                                                                                          <  2:1204730/138‑1 (MQ=255)
                                                      gagaATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTgg                                                                                                          >  1:1204730/1‑138 (MQ=255)
                                                      gagaATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCAccc                                                                                               >  1:528196/1‑149 (MQ=255)
                                                         aaTTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGcc                                                                                            >  1:1367297/1‑149 (MQ=255)
                                                           ttATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGt                                                                                                                      >  3:362794/1‑121 (MQ=255)
                                                           ttATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGt                                                                                                                      <  4:362794/121‑1 (MQ=255)
                                                               ccAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGc                                                                                      >  4:328734/1‑149 (MQ=255)
                                                                cAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCg                                                                                     >  4:544758/1‑149 (MQ=255)
                                                                cAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAAct                                                                                                                                            <  1:331821/94‑2 (MQ=255)
                                                                cAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAAct                                                                                                                                            >  2:331821/1‑93 (MQ=255)
                                                                   cATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTga                                                                                  <  3:526834/149‑1 (MQ=255)
                                                                     tatTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTgaga                                                                                >  4:1431780/1‑149 (MQ=255)
                                                                       tttAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAgg                                                                              >  4:1103173/1‑149 (MQ=255)
                                                                          aaCGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCaa                                                                           <  3:344474/149‑1 (MQ=255)
                                                                              ccAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtcgtc                                                                                                                                 <  3:770756/91‑1 (MQ=255)
                                                                              ccAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtcgtc                                                                                                                                 >  4:770756/1‑91 (MQ=255)
                                                                              ccAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTga                                                                                  >  1:720355/1‑138 (MQ=255)
                                                                              ccAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTga                                                                                  <  2:720355/138‑1 (MQ=255)
                                                                              ccAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGc                                                                       >  4:825928/1‑149 (MQ=255)
                                                                               cAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGcg                                                                      <  2:528196/149‑2 (MQ=255)
                                                                               cAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCa                                                                      <  2:1367297/149‑1 (MQ=255)
                                                                                 gTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGa                                                                    <  2:1302983/149‑1 (MQ=255)
                                                                                 gTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGa                                                                    <  3:1103173/149‑1 (MQ=255)
                                                                                  tCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTgg                                                                                                          >  3:97959/1‑110 (MQ=255)
                                                                                     cTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGt                                                                <  3:1431780/149‑1 (MQ=255)
                                                                                       ggATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGa                                                              <  2:1255508/149‑1 (MQ=255)
                                                                                         aTACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtc                                                                                                                                    >  2:845437/1‑77 (MQ=255)
                                                                                         aTACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtc                                                                                                                                    <  1:256254/77‑1 (MQ=255)
                                                                                         aTACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtc                                                                                                                                    <  1:845437/77‑1 (MQ=255)
                                                                                         aTACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtc                                                                                                                                    >  2:256254/1‑77 (MQ=255)
                                                                                              aTCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTgaga                                                                                >  1:1459018/1‑124 (MQ=255)
                                                                                              aTCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTgaga                                                                                <  2:1459018/124‑1 (MQ=255)
                                                                                                 cGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTa                                                    <  4:1371542/149‑1 (MQ=255)
                                                                                                  gTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTac                                                   <  3:544758/149‑1 (MQ=255)
                                                                                                   tAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTga                                                                                  >  1:345360/1‑117 (MQ=255)
                                                                                                   tAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTga                                                                                  <  2:345360/117‑1 (MQ=255)
                                                                                                     aaTGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTaccac                                                >  1:799211/1‑149 (MQ=255)
                                                                                                     aaTGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTaccac                                                >  3:414542/1‑149 (MQ=255)
                                                                                                        gctgcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACcaccac                                             <  4:414542/149‑1 (MQ=255)
                                                                                                         ctgcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTgaga                                                                                <  3:1435883/113‑1 (MQ=255)
                                                                                                         ctgcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTgaga                                                                                >  4:1435883/1‑113 (MQ=255)
                                                                                                          tgcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCAccc                                                                                               >  1:1116573/1‑97 (MQ=255)
                                                                                                          tgcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCAcgc                                           >  1:731701/1‑149 (MQ=255)
                                                                                                          tgcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCAcgc                                           >  4:324066/1‑149 (MQ=255)
                                                                                                               cTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGTCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCa                                      >  4:871225/1‑149 (MQ=255)
                                                                                                                tttGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCAt                                     >  2:1360182/1‑149 (MQ=255)
                                                                                                                 ttGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATc                                    <  3:640907/149‑1 (MQ=255)
                                                                                                                 ttGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATc                                    <  3:324066/149‑1 (MQ=255)
                                                                                                                           gCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGc                          <  3:58209/149‑1 (MQ=255)
                                                                                                                           gCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGc                          >  3:967133/1‑149 (MQ=255)
                                                                                                                             ccAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTga                                                                                  >  1:1108253/1‑91 (MQ=255)
                                                                                                                             ccAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTga                                                                                  <  2:1108253/91‑1 (MQ=255)
                                                                                                                             ccAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCAt                        >  1:607236/1‑149 (MQ=255)
                                                                                                                               aGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATaa                      <  2:752278/149‑1 (MQ=255)
                                                                                                                                gTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAAt                     <  2:731701/149‑1 (MQ=255)
                                                                                                                                 tGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATc                    <  2:607236/149‑1 (MQ=255)
                                                                                                                                      gTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTacc                                                  >  3:1209274/1‑114 (MQ=255)
                                                                                                                                      gTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGt               >  3:669212/1‑149 (MQ=255)
                                                                                                                                      gTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGt               <  4:428157/149‑1 (MQ=255)
                                                                                                                                       tGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTga                                                                                  <  3:383605/81‑1 (MQ=255)
                                                                                                                                        gAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGc             >  3:1458704/1‑149 (MQ=255)
                                                                                                                                        gAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGc             >  2:671704/1‑149 (MQ=255)
                                                                                                                                        gAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGc             >  4:1273120/1‑149 (MQ=255)
                                                                                                                                          ccgTTTCCAGACTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTac                                                   >  3:215770/3‑109 (MQ=255)
                                                                                                                                          ccgTTTCCAGACTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTac                                                   <  4:215770/107‑1 (MQ=255)
                                                                                                                                           aGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGt          <  4:669212/149‑1 (MQ=255)
                                                                                                                                            gTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTac                                                   <  2:353355/107‑1 (MQ=255)
                                                                                                                                            gTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTac                                                   >  1:353355/1‑107 (MQ=255)
                                                                                                                                             tttCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAATTCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCa        >  4:127715/1‑149 (MQ=255)
                                                                                                                                             tttCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCa        <  3:871225/149‑1 (MQ=255)
                                                                                                                                                   gTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGc             >  4:1357428/1‑138 (MQ=255)
                                                                                                                                                   gTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGAt  <  1:1360182/149‑1 (MQ=255)
                                                                                                                                                    tCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTac                                                   >  2:1081015/1‑99 (MQ=255)
                                                                                                                                                    tCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTac                                                   <  1:1081015/99‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGAT  >  NC_000913/2151821‑2152116

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