Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,377,791 A→T 16.0% intergenic (‑69/+20) blc ← / ← ampC outer membrane lipoprotein cell division and growth lipocalin/penicillin‑binding protein; beta‑lactamase, intrinsically weak

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,377,7910AT16.0% 508.3 / 82.3 207intergenic (‑69/+20)blc/ampCouter membrane lipoprotein cell division and growth lipocalin/penicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (93/80);  new base T (10/23);  total (103/104)
Fisher's exact test for biased strand distribution p-value = 2.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-01

Only 100 of 214 total aligned reads displayed.

TCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACG  >  NC_000913/4377692‑4377889
                                                                                                   |                                                                                                  
 cGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGaaaa                                                                                                   <  1:699746/100‑1 (MQ=255)
 cGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGaaa                                                                                                    >  1:2213307/1‑99 (MQ=255)
  gCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAAt                                                                                                  >  1:598558/1‑100 (MQ=255)
   cTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAAtt                                                                                                 <  2:2032651/100‑1 (MQ=255)
    tGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTc                                                                                                <  1:1487296/100‑1 (MQ=255)
     gCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCg                                                                                               <  2:1012243/100‑1 (MQ=255)
     gCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCg                                                                                               >  2:2713703/1‑100 (MQ=255)
     gCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCg                                                                                               <  2:3524487/100‑1 (MQ=255)
       ggCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGa                                                                                             >  2:2849416/1‑100 (MQ=255)
        gCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAc                                                                                            >  2:2184453/1‑100 (MQ=255)
         cAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAcc                                                                                           <  2:2870544/100‑1 (MQ=255)
         cAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAcc                                                                                           >  2:1168066/1‑100 (MQ=255)
          aaCGAGAGGGGGCAGGCGCATAAATGTTTCCTTGCTGGTTTTTTTCTAAGTGTGGATGACAGCAAGGAAGAGCGGGGAAAGGGTCCGAAGTTTCGGAcct                                                                                          <  1:1030902/100‑2 (MQ=255)
          aaCGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcct                                                                                          <  2:662811/100‑2 (MQ=255)
          aaCGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAccc                                                                                          >  2:395486/1‑100 (MQ=255)
          aaCGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAccc                                                                                          >  2:2671297/1‑100 (MQ=255)
           aCGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcctt                                                                                         <  1:2311754/100‑3 (MQ=255)
            cGAGAGGGAGTAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttt                                                                                        <  1:2859483/100‑4 (MQ=255)
            cGAGAGGGAGCGGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGGAAAGGTCCGAATTTTCGGAccttt                                                                                        <  1:2522580/100‑4 (MQ=255)
            cGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGa                                                                                        >  1:1735195/1‑100 (MQ=255)
              agagGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttttc                                                                                      >  1:1612376/1‑95 (MQ=255)
                agGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttttctc                                                                                    <  1:3163439/100‑8 (MQ=255)
                agGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttttctc                                                                                    <  2:1612383/100‑8 (MQ=255)
                 gggAGCAGGCGCATAAATGTTTCCTTGCTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttttctc                                                                                    <  2:3789448/99‑8 (MQ=38)
                 gggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGaa                                                                                   >  2:974970/1‑100 (MQ=255)
                 gggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGaa                                                                                   <  1:3334351/100‑1 (MQ=255)
                 gggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGaa                                                                                   >  2:4262854/1‑100 (MQ=255)
                  ggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAt                                                                                  <  2:2564856/100‑1 (MQ=255)
                  ggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAt                                                                                  >  2:3566907/1‑100 (MQ=255)
                   gAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGACAAGCGGAGAAAAAGACCGAAAATTCGGACCCGATGGAAtt                                                                                 >  1:1340101/1‑100 (MQ=255)
                     gCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAtttt                                                                               <  2:1096099/100‑1 (MQ=255)
                      cAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTa                                                                              <  1:3429340/100‑1 (MQ=255)
                        ggCGCATAAATGTTTCCTTACTGGTTTTTTTTTTTTTTTTGATGACAGCAAGGAAAAGCTGAGTAAAGGTCGGAAAATTCGGACCCGATGGAATTTTACt                                                                            >  1:2460951/1‑100 (MQ=255)
                        ggCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGATAAGGTCCGAAAATTCGGACCCGATGGAATTTTACt                                                                            >  1:2406940/1‑100 (MQ=255)
                        ggCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTaat                                                                            >  1:3725101/1‑98 (MQ=255)
                        ggCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAAGTCCGAAAATTCGGACCCGATGGAATTTTACt                                                                            >  1:3205637/1‑100 (MQ=255)
                         gcgcATAAATGTTTCCTGACTGGTTTTTTTCTAAGTGTAGAAGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGATCCGATGGAATTTTACTg                                                                           <  1:1819275/100‑1 (MQ=255)
                          cgcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGt                                                                          <  2:2334073/100‑1 (MQ=255)
                           gcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTa                                                                         >  2:2807352/1‑100 (MQ=255)
                           gcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATCTCGGCGATCGCCGGCTCATTAAAAGTCAGAAAATTCGGACCCGATGGAATTTTACTGTa                                                                         >  1:1561038/1‑100 (MQ=255)
                            cATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTa                                                                         >  1:1757732/1‑99 (MQ=255)
                              tAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTagag                                                                      <  1:887027/100‑1 (MQ=255)
                               aaaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGc                                                                     <  1:2011156/100‑1 (MQ=255)
                               aaaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGc                                                                     <  2:4412737/100‑1 (MQ=255)
                                aaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCg                                                                    <  1:4333492/100‑1 (MQ=255)
                                     ttCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAg                                                               <  1:3727114/100‑1 (MQ=255)
                                        cTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAAt                                                            >  2:488029/1‑100 (MQ=255)
                                        cTTACTGGTTTTTTTCTAAGTGTAGATCACAGCTAGGAAAGAAGGAGTAAAGGTCAGAAAATTCGGCCTCGATGGGATCTTACAGTATAGCGTTaaaga                                                             >  1:3992965/1‑96 (MQ=255)
                                          ttcTGGTTTTTTTGTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCt                                                          >  2:2985232/3‑100 (MQ=255)
                                          tACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAACATCGGATAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAATAATCt                                                          >  1:3360808/1‑100 (MQ=255)
                                          tACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCt                                                          >  1:4261428/1‑100 (MQ=255)
                                            ttccgatctttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGc                                                        <  2:2453132/92‑1 (MQ=255)
                                            cTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGTAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGc                                                        >  1:3908782/1‑100 (MQ=255)
                                            cTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTg                                                         >  1:1658103/1‑99 (MQ=255)
                                            cTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGc                                                        <  2:822862/100‑1 (MQ=255)
                                            cTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGc                                                        >  2:1660216/1‑100 (MQ=255)
                                                 ttttttCTAAGTGTAGATGACAGCAAGGAACAGATGATTAAAGATAAGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTACCAgg                                                    >  1:4282387/1‑99 (MQ=255)
                                                       cTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                             <  2:967140/100‑1 (MQ=255)
                                                       cTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                             >  1:3761518/1‑100 (MQ=255)
                                                          aGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGAc                                          <  1:2631033/100‑1 (MQ=255)
                                                          aGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGAc                                          <  2:1406833/100‑1 (MQ=255)
                                                           gtgtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGAct                                         >  2:2771786/1‑100 (MQ=255)
                                                            tgtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGActc                                        <  2:943260/100‑1 (MQ=255)
                                                              tAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACtctc                                      <  2:468184/100‑1 (MQ=255)
                                                               aGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCg                                     <  2:193612/100‑1 (MQ=255)
                                                                gATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGc                                    <  1:280726/100‑1 (MQ=255)
                                                                gATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATATGCCAGGCGGCGTCGACTCTCg                                     >  2:4165600/1‑99 (MQ=255)
                                                                 aTGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGACTCTGCCAGGCCGCGTCGACTCTCGc                                    >  2:4338440/1‑99 (MQ=255)
                                                                    aCAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGa                                >  1:1221529/1‑100 (MQ=255)
                                                                     cAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAt                               >  1:420097/1‑100 (MQ=255)
                                                                     cAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAt                               >  1:3399601/1‑100 (MQ=255)
                                                                         aaGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTggg                           <  1:3595739/100‑1 (MQ=255)
                                                                          aGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGa                          <  1:407701/100‑1 (MQ=255)
                                                                          aGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTGACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGa                          <  2:4104541/100‑1 (MQ=255)
                                                                           ggAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGAt                         >  2:4361341/1‑100 (MQ=255)
                                                                           ggAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGAt                         >  2:1376925/1‑100 (MQ=255)
                                                                                aGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttt                    >  2:4232535/1‑100 (MQ=255)
                                                                                   ggAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttg                 <  1:1142193/100‑1 (MQ=255)
                                                                                    tccatcgggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGACAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgt                >  2:3241386/8‑100 (MQ=255)
                                                                                    gagaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgt                <  2:1757739/100‑1 (MQ=255)
                                                                                       atcgggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTTTCGCTGGATTGGGATAGTTtttgttt              >  2:4139196/5‑99 (MQ=255)
                                                                                       aaaaGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg             <  2:878140/100‑1 (MQ=255)
                                                                                       aaaaGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg             <  2:3398321/100‑1 (MQ=255)
                                                                                        aaaGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCTATAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgttt              >  1:832139/1‑98 (MQ=255)
                                                                                        aaaGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGc            >  2:2497447/1‑100 (MQ=255)
                                                                                         aaGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGcc           <  2:3598091/100‑1 (MQ=255)
                                                                                         aaGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGcc           <  1:3459446/100‑1 (MQ=255)
                                                                                          gggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGc            >  1:2666854/2‑98 (MQ=255)
                                                                                          gggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCa          >  1:3796008/2‑100 (MQ=255)
                                                                                           ggTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAg         >  2:3099405/1‑100 (MQ=255)
                                                                                           ggTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAg         <  2:1235498/100‑1 (MQ=255)
                                                                                             tCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCa       >  2:1556403/1‑100 (MQ=255)
                                                                                             tCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCa       >  2:906651/1‑100 (MQ=255)
                                                                                               cGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATc     >  1:647916/1‑100 (MQ=255)
                                                                                               cGAAAAGTCGGGCCCGAGGGAATTCTACTGTAGAGCGTTAAGGAGCTGCCAGGCGGCGTCCACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCa          <  1:2512210/95‑1 (MQ=255)
                                                                                                gAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCa    >  1:522467/1‑100 (MQ=255)
                                                                                                gAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCa    >  1:4364435/1‑100 (MQ=255)
                                                                                                gAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCa    >  2:2956081/1‑100 (MQ=255)
                                                                                                 aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCAc   <  2:2946865/100‑1 (MQ=255)
                                                                                                  aaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACg  <  1:2990274/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACG  >  NC_000913/4377692‑4377889

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