Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,151,969 T→A 17.5% 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,9690TA17.5% 299.3 / 49.9 127R168S (AGA→AGTyegKser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base T (65/38);  new base A (12/10);  total (78/49)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

Only 100 of 131 total aligned reads displayed.

TGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACC  >  NC_000913/2151871‑2152068
                                                                                                  |                                                                                                   
tGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTc                                                                                                    <  2:1794218/100‑1 (MQ=255)
  aggagAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTg                                                                                                  <  2:1144000/100‑1 (MQ=255)
  aggagAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTg                                                                                                  <  2:3220202/100‑1 (MQ=255)
  aggagAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCAGCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTg                                                                                                  <  2:1453257/100‑1 (MQ=255)
   ggagAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTgg                                                                                                 <  1:879894/100‑1 (MQ=255)
   ggagAAGTATCCAGCATATTTAACGCCAGTCGCTCGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTgg                                                                                                 <  2:1921072/100‑1 (MQ=255)
        tctATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAAct                                                                                            <  2:1063037/98‑2 (MQ=255)
         ttATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACgg                                                                                           <  1:1478466/100‑1 (MQ=255)
           aTCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActttn                                                                                         <  2:268551/100‑5 (MQ=255)
           aTCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActttc                                                                                         <  1:158873/100‑5 (MQ=255)
           aTCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActaga                                                                                         >  1:1063045/1‑96 (MQ=255)
            tCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActttca                                                                                        <  1:2638398/100‑6 (MQ=255)
            tCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActttca                                                                                        <  1:2752999/100‑6 (MQ=255)
            tCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActttca                                                                                        <  1:2542989/100‑6 (MQ=255)
            tCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActttca                                                                                        <  2:420079/100‑6 (MQ=255)
            tCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTCCTGGTGAAAGTTTCCAGTCTGGAAACGGCAt                                                                                        <  2:163620/100‑1 (MQ=255)
              cAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCTAGTGCTGGTGAAAGTTTCCAGACTGGAAActttcacc                                                                                      <  2:1094541/100‑8 (MQ=255)
              cAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCt                                                                                      >  1:2481340/1‑100 (MQ=255)
              cAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCt                                                                                      >  2:1420358/1‑100 (MQ=255)
               aGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCtt                                                                                     >  1:3040733/1‑100 (MQ=255)
                gCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActttcaccag                                                                                    <  1:2731967/100‑10 (MQ=39)
                gCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCACAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtc                                                                                    >  1:1802223/1‑100 (MQ=255)
                acATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCACTGCTGGTGAAAGTTTCCAATCTGGAAACGGCATCTtc                                                                                    >  2:2824208/2‑100 (MQ=255)
                   tatTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtcgtc                                                                                 <  2:1034872/100‑1 (MQ=255)
                   tatTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtcgtc                                                                                 >  1:806083/1‑100 (MQ=255)
                   tatTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGGGAAAGTTTCCAGTCTGGAAACGGCATCTtcgtc                                                                                 >  2:3336021/1‑100 (MQ=255)
                   tatTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGGGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtcgtc                                                                                 >  2:942757/1‑100 (MQ=255)
                     tttAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCgg                                                                               >  1:2594878/1‑100 (MQ=255)
                     tttAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCgg                                                                               <  1:2666725/100‑1 (MQ=255)
                     tttAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGt                                                                              <  1:396422/100‑1 (MQ=255)
                       tAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTa                                                                             <  2:2676921/100‑1 (MQ=255)
                       tAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTa                                                                             >  1:1243852/1‑100 (MQ=255)
                           gCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGa                                                                         >  1:3577664/1‑100 (MQ=255)
                            ccAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGaa                                                                        >  2:1030292/1‑100 (MQ=255)
                              aGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGt                                                                      <  1:1267138/100‑1 (MQ=255)
                                    tGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCa                                                                >  2:3470843/1‑100 (MQ=255)
                                    tGGATACCATCCGTAAATGCTGCAACTTAGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCa                                                                >  1:940551/1‑100 (MQ=255)
                                     ggATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCAt                                                               >  1:1956266/1‑100 (MQ=255)
                                     ggATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCAt                                                               >  1:922451/1‑100 (MQ=255)
                                        tACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGtt                                                            >  1:821193/1‑100 (MQ=255)
                                        tACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGtt                                                            >  1:3445946/1‑100 (MQ=255)
                                        tACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGtt                                                            >  1:2818356/1‑100 (MQ=255)
                                          ccATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTgg                                                          >  1:439707/1‑100 (MQ=255)
                                            cttccgatctgctgcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCa                                                        <  2:2061533/91‑1 (MQ=255)
                                            aTCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCa                                                        >  1:1680754/1‑100 (MQ=255)
                                            aTCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCa                                                        >  2:634854/1‑100 (MQ=255)
                                             tCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCat                                                       >  1:1892016/1‑100 (MQ=255)
                                               cGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCatat                                                     <  1:591742/100‑1 (MQ=255)
                                                gTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATAtt                                                    >  1:1373404/1‑100 (MQ=255)
                                                  aaaTGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCa                                                  <  2:1483313/100‑1 (MQ=255)
                                                      gctgcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCaagatcggaag                                              >  1:2061545/1‑90 (MQ=255)
                                                      gctgcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCg                                              >  1:285911/1‑100 (MQ=255)
                                                       ctgcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCAGATTCACCCGc                                             >  1:113318/1‑100 (MQ=255)
                                                         gcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTTCAGTCTGGAAACGGCCTCTTCGTCGGTAATGCAGTGAGTCATGTTGGCATATTCACCCGCCa                                           >  2:1316974/1‑100 (MQ=255)
                                                         gcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCa                                           >  2:2296943/1‑100 (MQ=255)
                                                         gcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCa                                           >  1:497757/1‑100 (MQ=255)
                                                          cAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCAt                                          >  1:1292788/1‑100 (MQ=255)
                                                           aaCTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATc                                         <  1:2870714/100‑1 (MQ=255)
                                                            aCTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCg                                        >  1:2409129/1‑100 (MQ=255)
                                                            aCTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGGCTGGAAACGGCATCTTCGGCGGGAAGGAAGTGCGTCATGTTGGCATATTCACCCGCCATCg                                        >  1:1983501/1‑100 (MQ=255)
                                                              tttGTGCGCACGCCCAGTGCTGGTGAAAGTTTCAAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCgg                                       <  1:505269/99‑1 (MQ=255)
                                                              tgtgtgCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCgg                                       <  1:1295513/97‑1 (MQ=255)
                                                               ttGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCg                                     <  2:2463583/100‑1 (MQ=255)
                                                                    cgcACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTgaga                                <  2:2317169/100‑1 (MQ=255)
                                                                    cgatctcccAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTgaga                                <  2:1434363/94‑1 (MQ=255)
                                                                     gcACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCTTCATGTTGGCATATTCACCCGCCATCGGCGTGagag                               <  1:2574961/100‑1 (MQ=255)
                                                                      cACGCCCAGTGCTGGTGAAAGTTTCCAGTCTCGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAgg                              >  2:2309717/1‑100 (MQ=255)
                                                                       aCGCCCAGTGGTGGTGAAAGTTTCCAGTCTGGAACCGGCATCTTCGTCGGTAATGAAGGGCGTCATGTTGGCAGGTTCGCCCGCCATCGGCGTGAGAGGc                             <  2:871370/100‑1 (MQ=255)
                                                                        cGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCa                            >  2:2369127/1‑100 (MQ=255)
                                                                          cccAGTGCTGGTGAAAGTTTCCAGTCTGGAAGCGGCATCTTCGTCGGTAATGGAGTGCGTCATGTTGGAATATTCACCCGCCATCGGGGTGGGAGGAAAt                          >  2:376558/1‑100 (MQ=255)
                                                                          cccAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTgagaagatag                          >  1:1434371/1‑94 (MQ=255)
                                                                           ccAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAAtt                         >  1:916617/1‑100 (MQ=255)
                                                                            cAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTg                        <  2:2150426/100‑1 (MQ=255)
                                                                             aGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCCATTGc                       >  2:1665414/1‑100 (MQ=255)
                                                                             aGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGc                       <  1:1100258/100‑1 (MQ=255)
                                                                             aGTGCTGGTGAAAGCTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGc                       <  2:2481324/100‑1 (MQ=255)
                                                                               tGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAg                     >  1:2157897/1‑100 (MQ=255)
                                                                                 cTGGTGAAAGTTTCCAGTCTGGAAACGGTATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGAc                   <  1:314460/100‑1 (MQ=255)
                                                                                   ggTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCg                 >  1:558049/1‑100 (MQ=255)
                                                                                    gTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGt                >  1:1878852/1‑100 (MQ=255)
                                                                                    gTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGt                >  1:2215937/1‑100 (MQ=255)
                                                                                     tGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTg               >  1:2716638/1‑100 (MQ=255)
                                                                                     atgccgTTTCCAGACTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTg               >  1:578610/6‑100 (MQ=255)
                                                                                     atgccgTTTCCAGACTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTg               >  1:3059263/6‑100 (MQ=255)
                                                                                      tgccgTTTTCAGACTGGAACCGGCATCTTCGTCGGTAATGAATTGCGTCGTTTTGGCGTATTCGCCCGCCATCGGCGTGAGAGGCAATTGCAGAcc                  >  2:1626592/5‑96 (MQ=255)
                                                                                       gccgTTTCCAGACTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGTGTGAGAGGCAATTGCAGACCGTGAc             >  2:3110353/4‑100 (MQ=255)
                                                                                        ccgTTTCCAGACTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGAcc            >  1:3596821/3‑100 (MQ=255)
                                                                                        ccgTTTCCAGACTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGAcc            >  1:139850/3‑100 (MQ=255)
                                                                                        ccgTTTCCAGACTGGAAACGCCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGAcc            >  2:2218376/3‑100 (MQ=255)
                                                                                        aaGTTGCCGGTCTGGAAACGGCATCTTCGTCGGTAATGAGGTGCGTGATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGAc             <  2:269523/99‑1 (MQ=255)
                                                                                          gTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCa          <  2:2691804/100‑1 (MQ=255)
                                                                                          gTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCa          >  1:1834864/1‑100 (MQ=255)
                                                                                           tttCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCaa         >  1:2088653/1‑100 (MQ=255)
                                                                                            ttCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGCGAGAGGCAATTGCAGACCGTGACCCaaa        >  2:2567680/1‑100 (MQ=255)
                                                                                               cAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCt     >  2:2316422/1‑100 (MQ=255)
                                                                                                aGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTa    <  2:806077/100‑1 (MQ=255)
                                                                                                aGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGACGTGACAGGCAATCGCAGACCGTTACCCaaaac      >  2:3331073/1‑96 (MQ=255)
                                                                                                  tCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTacc  <  2:1124237/100‑1 (MQ=255)
                                                                                                  tCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTacc  <  2:1662683/100‑1 (MQ=255)
                                                                                                  tCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTacc  <  2:1663097/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACC  >  NC_000913/2151871‑2152068

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