Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 201,032 G→A 6.2% G21D (GGC→GAC)  lpxD → UDP‑3‑O‑(3‑hydroxymyristoyl)‑glucosamine N‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913201,0320GA 6.2% 181.9 / 2.7 65G21D (GGC→GAC) lpxDUDP‑3‑O‑(3‑hydroxymyristoyl)‑glucosamine N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (28/32);  new base A (2/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

GGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTC  >  NC_000913/200886‑201174
                                                                                                                                                  |                                                                                                                                              
ggTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCg                                                                                                                                              >  1:529682/1‑149 (MQ=255)
    gttgATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGatat                                                                                                                                          >  1:549153/1‑149 (MQ=255)
        ttGCAAACGCCGGTGCTTACAACAGGAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGtc                                                                                                                                      <  2:241697/148‑1 (MQ=255)
        aTGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCTCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGtc                                                                                                                                      <  1:378412/149‑1 (MQ=255)
        aTGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGtc                                                                                                                                      <  2:774678/149‑1 (MQ=255)
          gCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGTCTAAATAAGTAATGCCTTCAAATCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGtcat                                                                                                                                    <  2:124888/149‑1 (MQ=255)
            aaaCGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGACGATATCGccatca                                                                                                                                  <  1:833802/149‑1 (MQ=255)
              aCGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCAcc                                                                                                                                >  2:629210/1‑149 (MQ=255)
                 ccGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGACGATATCGCCATCACCgtg                                                                                                                             <  2:245665/149‑3 (MQ=255)
                   gTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGt                                                                                                                           <  1:600917/149‑1 (MQ=255)
                    ttGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGtt                                                                                                                          <  1:696980/149‑1 (MQ=255)
                    gtGCTTACACCAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCACTTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGtt                                                                                                                          <  2:273107/148‑1 (MQ=255)
                      gCTTACAACACCAGCGATGTAAAAGACATCACTGCCGACGTACGGAAACAGGTTAAATAAGTAATGCCTTCAAATCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGc                                                                                                                        <  1:461817/149‑1 (MQ=255)
                        ttACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGt                                                                                                                      <  1:110504/149‑1 (MQ=255)
                         tACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTc                                                                                                                     >  1:837100/1‑149 (MQ=255)
                           aaacaGCAGATATGTAAAAGACATCACTGCCGACGTAATGAAACATGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGTACTACACGGTGATGGCGATAGCGTCATCACCGGCGTTGAGTCCa                                                                                                                   <  1:613229/148‑1 (MQ=255)
                            aacaGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTGGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCAt                                                                                                                  <  2:637726/149‑1 (MQ=255)
                                gcagcGAGGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCaa                                                                                                              <  2:750427/149‑1 (MQ=255)
                                      aTGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGca                                                                                                        <  2:833/149‑1 (MQ=255)
                                       tGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGcac                                                                                                       >  2:915154/1‑149 (MQ=255)
                                          aaaaGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCATCAGTTGGATGCAGAACTACACGGTGATGGCGATCTCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACaaa                                                                                                    >  1:393685/1‑149 (MQ=255)
                                          aaaaGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACaaa                                                                                                    <  1:745030/149‑1 (MQ=255)
                                           aaaGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAAc                                                                                                   >  2:137289/1‑149 (MQ=255)
                                             aGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAg                                                                                                 >  2:488764/1‑149 (MQ=255)
                                                  tCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTcac                                                                                            <  1:488764/149‑1 (MQ=255)
                                                   cACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTcaca                                                                                           <  1:481467/149‑1 (MQ=255)
                                                      tGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTa                                                                                        >  1:510328/1‑149 (MQ=255)
                                                      tGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATGCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGGCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTa                                                                                        >  2:698037/1‑149 (MQ=255)
                                                         cGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGt                                                                                     <  1:880956/149‑1 (MQ=255)
                                                          gACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATCGCGATATCGTCATCACCGGAGTTGCGTACATGCAATCTGCACAAACAGGTCACATTACGtt                                                                                    >  2:37592/1‑149 (MQ=255)
                                                                 tGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGtt                                                                             >  1:603510/1‑149 (MQ=255)
                                                                 tGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGAAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCAGCACAAACAGGTCACATTACGTTCATGGtt                                                                             >  1:194150/1‑149 (MQ=255)
                                                                  gAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTa                                                                            >  1:753866/1‑149 (MQ=255)
                                                                   aaaCAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGTCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTaa                                                                           >  2:819879/1‑149 (MQ=255)
                                                                   aaaCAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTaa                                                                           >  1:676705/1‑149 (MQ=255)
                                                                      cAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTCCACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACAGTACGTTCATGGTTAAccc                                                                        >  2:498053/1‑149 (MQ=255)
                                                                        gtttAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGGCACATTACGTTCATGGTTAACCCaa                                                                      <  2:262082/147‑1 (MQ=255)
                                                                          ttAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCaaaa                                                                    <  1:291696/149‑1 (MQ=255)
                                                                           tAAATAAGTAATTCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGTCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCAGATTACGTTCATGGTTAACCCAAAAt                                                                   <  2:911198/149‑1 (MQ=255)
                                                                                  gTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGa                                                            >  1:15270/1‑149 (MQ=255)
                                                                                   tAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAg                                                           <  1:49348/149‑1 (MQ=255)
                                                                                   tAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCCAAATCGTCATCACCGGCGTTGCTTCAAAGAAATCTGAAAAAACAGCTCAAAATACTTTCATGCTTAACCCAAACTAcctagag                                                           >  1:668832/1‑144 (MQ=255)
                                                                                        ccTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCAttt                                                      <  1:972781/149‑1 (MQ=255)
                                                                                           tCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAgg                                                   <  2:529682/149‑1 (MQ=255)
                                                                                              aTTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACAAGGGGATGGCGATATCGTCATCACCGGCGGTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCtt                                                <  2:846158/149‑1 (MQ=255)
                                                                                                  gACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGGCACATTACGTGCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGc                                            >  1:719043/1‑149 (MQ=255)
                                                                                                    cTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCa                                          <  1:20921/149‑1 (MQ=255)
                                                                                                       gCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGc                                       <  2:231704/149‑1 (MQ=255)
                                                                                                             ttAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCcgc                                 <  2:549153/149‑1 (MQ=255)
                                                                                                                 ggcagcaGTTGTCTGCAGAACTACACGGGGATGGCGATATCGTCATCACCGGCGTTGCGTCGATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACACAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGtt                             <  2:393685/148‑1 (MQ=255)
                                                                                                                 ggcagcaGTTGGATGCAGAACTACACGGTGAGGGCGATAGCGTCATCACCGGCGTTGCGTCCATGCAAGCTGCACTAACAGGTCACATTACGTTCATGGTTCACCCAAACTACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGtt                             <  2:890430/148‑1 (MQ=255)
                                                                                                                 ctcagcagTTGTATTCAGAAATACACGCTGAAGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTATCCCAAAATACCGTGAGCTTTTAGGCTTGTGCCAGGCGTCCGCGGtt                             <  2:477030/147‑1 (MQ=255)
                                                                                                                 cgcAGCAGTTGGATGCAGAACTACACGGTGATGGCGATCTCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGtt                             <  1:879522/149‑1 (MQ=255)
                                                                                                                  gcagcaGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTACACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTg                            >  1:796491/1‑149 (MQ=255)
                                                                                                                       agtacgATACAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATCCAATCAGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTGGTGCCAGGCGACCGCGGTAGTCATg                       >  2:249807/6‑149 (MQ=255)
                                                                                                                               gCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGa               >  1:281497/1‑149 (MQ=255)
                                                                                                                                 aGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCCCAAACAGGTCACATTACGTTCCTGGTTAACCCAAAATACCGTGCGCTTTGAGGGTGGGGCCCGGCGGCGCGGGTGGGCAGGCCCCAGGAct             >  2:624672/1‑148 (MQ=255)
                                                                                                                                 aGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGaag             >  1:768426/1‑147 (MQ=255)
                                                                                                                                 aGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCAAGGACg             >  1:484813/1‑149 (MQ=255)
                                                                                                                                        gccGGTGATGACGATATCGCCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTcc      >  2:5882/3‑149 (MQ=255)
                                                                                                                                        cacGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTTATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTcc      >  1:22212/1‑149 (MQ=255)
                                                                                                                                          cGGTGATGACGATATCGCCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCtt    >  2:15805/1‑149 (MQ=255)
                                                                                                                                            gTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCTGGACGATCTTCCTTTc  >  2:823277/1‑149 (MQ=255)
                                                                                                                                            gTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTc  <  1:915154/149‑1 (MQ=255)
                                                                                                                                            gTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTc  <  1:137289/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTC  >  NC_000913/200886‑201174

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑