Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,377,792 A→T 18.4% intergenic (‑70/+19) 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,7920AT18.4% 313.7 / 64.5 136intergenic (‑70/+19)blc/ampCouter membrane lipoprotein cell division and growth lipocalin/penicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (61/50);  new base T (9/16);  total (70/66)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-01

Only 100 of 140 total aligned reads displayed.

CGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCA  >  NC_000913/4377693‑4377887
                                                                                                   |                                                                                               
cGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGaaaa                                                                                                 >  2:1786629/1‑100 (MQ=255)
 gCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAAt                                                                                                >  1:1486016/1‑100 (MQ=255)
   tGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTc                                                                                              >  2:2479199/1‑100 (MQ=255)
   tGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTc                                                                                              <  1:148588/100‑1 (MQ=255)
     cGGCAACGAGAGGGAGCATGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCgg                                                                                            >  1:1875745/1‑100 (MQ=255)
     cGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCgg                                                                                            <  1:2419168/100‑1 (MQ=255)
      ggCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGa                                                                                           >  2:1266111/1‑100 (MQ=255)
       gCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAc                                                                                          <  2:737671/100‑1 (MQ=255)
        cAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcc                                                                                         <  2:518617/100‑1 (MQ=255)
          aCGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcctt                                                                                       <  2:455848/100‑3 (MQ=255)
          aCGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcctt                                                                                       <  1:360638/100‑3 (MQ=255)
          aCGAAAGGGAGCAGGGGCATAAATGATTCCTTACTGGTTTTTGTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcctt                                                                                       <  1:1024237/100‑3 (MQ=255)
           cGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTTTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAcc                                                                                         >  1:1575314/1‑97 (MQ=255)
           cGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGTAAGGAAAAGATGAGTAAAGGTCCGAAAATTCGGAcc                                                                                         >  1:1836040/1‑97 (MQ=255)
           cGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttt                                                                                      <  2:2648230/100‑4 (MQ=255)
              gagGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttttct                                                                                   <  2:688820/100‑7 (MQ=255)
               agGGCGCAGGCGCATAAATGTTTGCTTACTTGTTTTTTTCGAAGTGTAGATGACTGCCAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGTccttttct                                                                                   <  1:1133270/99‑7 (MQ=255)
               agGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTTTTTTTTTATATGACAACAAGGAAAAACGGAGTAAAAATGCGAAAAAGCGGACCCGAttg                                                                                   >  1:2919508/1‑97 (MQ=255)
               agGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTACGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttttct                                                                                   <  1:1862952/99‑7 (MQ=255)
               agGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttttctc                                                                                  <  2:2155726/100‑8 (MQ=255)
                gggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttttctcc                                                                                 <  2:2128857/100‑9 (MQ=38)
                gggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGaa                                                                                 >  2:2004854/1‑100 (MQ=255)
                 ggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttttctccg                                                                                <  2:377831/100‑10 (MQ=38)
                 ggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttttctccg                                                                                <  1:10786/100‑10 (MQ=38)
                   aGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAttt                                                                              >  2:1631233/1‑100 (MQ=255)
                      aGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAc                                                                           <  2:2664774/100‑1 (MQ=255)
                      aGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAc                                                                           <  2:2304880/100‑1 (MQ=255)
                      aGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAc                                                                           <  2:2737905/100‑1 (MQ=255)
                       ggCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACt                                                                          >  2:32609/1‑100 (MQ=255)
                       ggCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACt                                                                          >  2:1783467/1‑100 (MQ=255)
                       ggCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGATCCGAAAATTCGGACCCGATGGAATTTTACt                                                                          >  1:2465709/1‑100 (MQ=255)
                        gcgcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTg                                                                         <  2:2311022/100‑1 (MQ=255)
                        gcgcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTg                                                                         <  2:245052/100‑1 (MQ=255)
                           cATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTa                                                                       >  1:2150534/1‑99 (MQ=255)
                              aaaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGc                                                                   >  1:145854/1‑100 (MQ=255)
                              aaaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGc                                                                   <  2:1197291/100‑1 (MQ=255)
                              aaaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGc                                                                   <  1:544114/100‑1 (MQ=255)
                               aaTGTTTCCTTACTGGTTTTTTTTTTTTTGTATATGACAGCAAGGAAAGTCGGAGAAATTGTTTTTAAATTGGGACCCGATGGAATTTTACTGTa                                                                       >  1:1435101/1‑95 (MQ=255)
                               aaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGcn                                                                  <  2:809267/100‑2 (MQ=255)
                               aaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCg                                                                  >  2:1656954/1‑100 (MQ=255)
                                aTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGt                                                                 <  2:1587454/100‑1 (MQ=255)
                                aTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGt                                                                 >  2:1621377/1‑100 (MQ=255)
                                  gTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTa                                                               >  2:1766192/1‑100 (MQ=255)
                                    ttCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCTAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAg                                                             <  2:2050563/100‑1 (MQ=255)
                                     tCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGTGAAAAGGTCCGAAAATTCGGACCCTATGGATTTTTACTGTATAGCGTTaa                                                              >  1:537997/1‑98 (MQ=255)
                                          aCTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTg                                                       <  2:1857974/100‑1 (MQ=255)
                                          aCTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTg                                                       <  1:2167181/100‑1 (MQ=255)
                                           cTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGc                                                      >  1:833504/1‑100 (MQ=255)
                                            tGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAATGTCAGAAAAGTCGGACCCCGTGAGATGTTACTGTGGAGCGTTGAGAAGCtac                                                      >  1:2398576/1‑97 (MQ=255)
                                            tGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGcc                                                     >  1:2131154/1‑100 (MQ=255)
                                             ggTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCa                                                    >  2:215132/1‑100 (MQ=255)
                                             ggTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCa                                                    >  1:216673/1‑100 (MQ=255)
                                               tttttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAgg                                                  >  2:2840896/1‑100 (MQ=255)
                                                ttttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggc                                                 >  2:529330/1‑100 (MQ=255)
                                                 tttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggcg                                                <  2:2012494/100‑1 (MQ=255)
                                                  ttttCTAAGTGTAGATGACAGCGAGGGCCGTCTGATTAAAAAACAGAAAATTCGGACACGATGGAATTTTACTGTAGAGCGTTAAAAATTTGCCAgac                                                 >  1:17103/1‑96 (MQ=255)
                                                  ttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggcgg                                               >  2:300943/1‑100 (MQ=255)
                                                  ttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggcgg                                               <  1:1266106/100‑1 (MQ=255)
                                                  ttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAGGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggcgg                                               <  1:2880927/100‑1 (MQ=255)
                                                     tCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGt                                            >  2:432328/1‑100 (MQ=255)
                                                       tAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAGTCTGACAGGCGGCGTCg                                          >  2:2418213/1‑100 (MQ=255)
                                                       tAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCg                                          <  2:1704446/100‑1 (MQ=255)
                                                          gtgtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGAct                                       >  2:218929/1‑100 (MQ=255)
                                                          gtgtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGAct                                       <  1:525337/100‑1 (MQ=255)
                                                           tgtAGATGACAGCAAAGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGTGCGTTAAGAATCTGCCAGGCGGCGTCGAct                                       <  1:486711/99‑1 (MQ=255)
                                                               gATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGc                                  >  2:612637/1‑100 (MQ=255)
                                                                aTGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCt                                 <  2:2261689/100‑1 (MQ=255)
                                                                aTGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCt                                 >  1:640564/1‑100 (MQ=255)
                                                                aTGACAGCAAGGAAAAGCGGAGAAAAGATCCGAAAATTCGGACCCGATGGAATTTTACCGTAGACCGTTAAGAATCTGCCAGGCGGCGTCGACTGTCACt                                 >  1:414921/1‑100 (MQ=255)
                                                                 tGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTg                                <  1:2350651/100‑1 (MQ=255)
                                                                  gACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTgg                               <  2:2171991/100‑1 (MQ=255)
                                                                    cAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGa                              >  1:1331089/1‑99 (MQ=255)
                                                                     aGCAAGTAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTg                                <  1:2701871/96‑1 (MQ=255)
                                                                      gCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTg                           >  2:1336047/1‑100 (MQ=255)
                                                                      gCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTg                           <  2:1637205/100‑1 (MQ=255)
                                                                       cAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGTCTCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTgg                          <  2:2801404/100‑1 (MQ=255)
                                                                       cAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGTAATTTTACGGTAGAGCTTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTgg                          <  2:1619443/100‑1 (MQ=255)
                                                                         acgAAAAGCGTAGAAAATGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCTTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTgg                          <  1:11624/96‑1 (MQ=255)
                                                                          ggAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGAt                       <  1:1241830/100‑1 (MQ=255)
                                                                          ggAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGAt                       >  1:2341750/1‑100 (MQ=255)
                                                                          ggAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGAt                       >  2:1619171/1‑100 (MQ=255)
                                                                                gCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGtttt                 >  1:2269362/1‑100 (MQ=255)
                                                                                  ggAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttg               >  2:820919/1‑100 (MQ=255)
                                                                                  ggAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttg               >  1:976691/1‑100 (MQ=255)
                                                                                  ggAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttg               <  1:846600/100‑1 (MQ=255)
                                                                                   gagaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgt              >  2:1507648/1‑100 (MQ=255)
                                                                                   gagaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgt              >  2:1353623/1‑100 (MQ=255)
                                                                                     catcgggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATCGTTtttggtt            >  1:528197/6‑100 (MQ=255)
                                                                                     catcgggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgttt            >  1:2014575/6‑100 (MQ=255)
                                                                                     catcgggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgttt            >  2:1590970/6‑100 (MQ=255)
                                                                                      atcgggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg           >  1:2289366/5‑100 (MQ=255)
                                                                                       tcgggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGc          >  1:708518/4‑100 (MQ=255)
                                                                                       aaaGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgttt            >  1:987993/1‑98 (MQ=255)
                                                                                          ggTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAg       >  2:12244/1‑100 (MQ=255)
                                                                                             ncgcaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCCGCAt    >  1:89/5‑100 (MQ=255)
                                                                                             ccGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCAt    >  1:1892356/1‑100 (MQ=255)
                                                                                             ccGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCAt    >  1:1070080/1‑100 (MQ=255)
                                                                                              cGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATc   >  1:350582/1‑100 (MQ=255)
                                                                                               gAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTGTTGTTTGCCAGCAtaa  >  1:2713382/1‑98 (MQ=255)
                                                                                                aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATc   <  1:2793832/98‑1 (MQ=255)
                                                                                                   |                                                                                               
CGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCA  >  NC_000913/4377693‑4377887

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: