Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,106,320 G→A *167* (TAG→TAA cpxP → inhibitor of the cpx response, periplasmic adaptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,106,3200GA98.7% 247.9 / NA 76*167* (TAG→TAAcpxPinhibitor of the cpx response, periplasmic adaptor protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (37/38);  minor base C (0/1);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTTAACGCCGGAGCAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAG  >  NC_000913/4106186‑4106450
                                                                                                                                      |                                                                                                                                   
cTGTTAACGCCGGAGCAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATaaa                                                                                                                                 <  2:1995615/139‑1 (MQ=255)
 tGTTAACGCCGGAGCAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAAc                                                                                                                                >  1:431069/1‑139 (MQ=255)
   ttAACGCCGGAGCAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAAccc                                                                                                                              >  2:1322912/1‑139 (MQ=255)
       cGCCGGAGCAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGtt                                                                                                                          <  1:611508/139‑1 (MQ=255)
           ggAGCAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTcc                                                                                                                      <  2:427810/139‑1 (MQ=255)
             agcagcAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCtt                                                                                                                    >  2:1258999/1‑139 (MQ=255)
             agcagcAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCtt                                                                                                                    <  2:431069/139‑1 (MQ=255)
             agcagcAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCtt                                                                                                                    >  1:412389/1‑139 (MQ=255)
                 gcaAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCa                                                                                                                <  1:1322912/139‑1 (MQ=255)
                   aaGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATa                                                                                                              <  1:1258999/139‑1 (MQ=255)
                   aaGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATa                                                                                                              <  2:828271/139‑1 (MQ=255)
                        gTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGacac                                                                                                         <  2:264369/139‑1 (MQ=255)
                                  agaAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTc                                                                                               <  2:412389/139‑1 (MQ=255)
                                   gaAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCt                                                                                              >  1:1648298/1‑139 (MQ=255)
                                        aTCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTcccc                                                                                         >  1:127319/1‑139 (MQ=255)
                                             caacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGacac                                                                                                         <  2:1637428/118‑1 (MQ=255)
                                             caacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGacac                                                                                                         >  1:1637428/1‑118 (MQ=255)
                                                  aaTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATcc                                                                                                    <  2:432886/118‑1 (MQ=255)
                                                  aaTGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATcc                                                                                                    >  1:432886/1‑118 (MQ=255)
                                                        gCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGt                                                                        >  2:1192960/1‑139 (MQ=255)
                                                             tGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGtttttt                                                                   <  2:127319/139‑1 (MQ=255)
                                                              gCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCcaca                                                                                     <  1:743352/121‑1 (MQ=255)
                                                              gCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCcaca                                                                                     >  2:743352/1‑121 (MQ=255)
                                                                 tgacgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTAt                                                               >  1:1063813/1‑139 (MQ=255)
                                                                   acgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTcc                                                                                                                      <  2:461025/83‑1 (MQ=255)
                                                                   acgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTcc                                                                                                                      >  1:461025/1‑83 (MQ=255)
                                                                      tgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCa                                                          >  1:15834/1‑139 (MQ=255)
                                                                         cgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAAtt                                                       >  1:2050081/1‑139 (MQ=255)
                                                                              tGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATg                                                                                   <  1:1504268/107‑1 (MQ=255)
                                                                              tGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATg                                                                                   >  2:1504268/1‑107 (MQ=255)
                                                                               ggCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCctg                                                 <  1:1192960/139‑1 (MQ=255)
                                                                                    aaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTgctgct                                            <  2:1648298/139‑1 (MQ=255)
                                                                                    aaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTgctgct                                            <  2:92391/139‑1 (MQ=255)
                                                                                           cATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCctg                                                 <  1:605846/127‑1 (MQ=255)
                                                                                           cATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCctg                                                 >  2:605846/1‑127 (MQ=255)
                                                                                                           gAGTAGTCGCAACTCACGTTTCCCGTTCTTCACCCTCTTTTTCTTGTCTTAGTCACATTTTCTTCATTTCCCCACTTCTTTTGTTGGTTTTTTTTATCCTCATTTTTTCTTCTGCTTAATGCATTCCAGATGATTTgtt                     <  1:834516/139‑3 (MQ=255)
                                                                                                                 tagCAACTCACGTTCCTAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAt                                    >  2:471352/1‑118 (MQ=255)
                                                                                                                 tagCAACTCACGTTCCTAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAt                                    <  1:471352/118‑1 (MQ=255)
                                                                                                                 tagCAACTCACGTTCCCAGTAATAAACCCTGTTTTcc                                                                                                                      >  2:504366/1‑37 (MQ=255)
                                                                                                                 tagCAACTCACGTTCCCAGTAATAAACCCTGTTTTcc                                                                                                                      <  1:504366/37‑1 (MQ=255)
                                                                                                                 tagCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATc                                                              <  1:1389790/92‑1 (MQ=255)
                                                                                                                 tagCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGtt               >  2:1819453/1‑139 (MQ=255)
                                                                                                                 tagCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCAGCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATc                                                              >  2:1389790/1‑92 (MQ=255)
                                                                                                                   gCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGacac                                                                                                         >  1:962750/1‑48 (MQ=255)
                                                                                                                   gCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGacac                                                                                                         <  2:962750/48‑1 (MQ=255)
                                                                                                                   gCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGt                                                                        <  2:524575/80‑1 (MQ=255)
                                                                                                                   gCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGt                                                                        >  1:524575/1‑80 (MQ=255)
                                                                                                                   gCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCctg                                                 <  1:1285920/103‑1 (MQ=255)
                                                                                                                   gCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCctg                                                 >  2:1285920/1‑103 (MQ=255)
                                                                                                                   gCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCa                                     >  1:1002016/1‑115 (MQ=255)
                                                                                                                   gCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCa                                     <  2:1002016/115‑1 (MQ=255)
                                                                                                                      aCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATa                                                                                                              >  1:1925007/1‑40 (MQ=255)
                                                                                                                      aCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATa                                                                                                              <  2:1925007/40‑1 (MQ=255)
                                                                                                                      aCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTcc                                                                                           <  1:236241/59‑1 (MQ=255)
                                                                                                                      aCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTcc                                                                                           >  2:236241/1‑59 (MQ=255)
                                                                                                                      aCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATAc           <  2:135123/139‑1 (MQ=255)
                                                                                                                      aCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATAc           >  1:135123/1‑139 (MQ=255)
                                                                                                                      aCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATc                                                              <  2:911205/87‑1 (MQ=255)
                                                                                                                      aCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATc                                                              >  1:911205/1‑87 (MQ=255)
                                                                                                                      aCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTc                                                           >  1:1860766/1‑90 (MQ=255)
                                                                                                                      aCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTc                                                           <  2:1860766/90‑1 (MQ=255)
                                                                                                                      aCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACt          >  1:1005789/1‑139 (MQ=255)
                                                                                                                       cTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCctg                                                 >  2:1384366/1‑99 (MQ=255)
                                                                                                                       cTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCctg                                                 <  1:1384366/99‑1 (MQ=255)
                                                                                                                        tCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTg                                  <  2:186731/113‑1 (MQ=255)
                                                                                                                        tCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTg                                  >  1:186731/1‑113 (MQ=255)
                                                                                                                          aCGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCg      <  2:1063813/139‑1 (MQ=255)
                                                                                                                           cGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGacac                                                                                                         <  2:1072196/40‑1 (MQ=255)
                                                                                                                           cGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGacac                                                                                                         >  1:1072196/1‑40 (MQ=255)
                                                                                                                           cGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTcc                                                                                           <  2:1938387/54‑1 (MQ=255)
                                                                                                                           cGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTcc                                                                                           >  1:1938387/1‑54 (MQ=255)
                                                                                                                           cGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCcaca                                                                                     <  2:1578780/60‑1 (MQ=255)
                                                                                                                           cGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCcaca                                                                                     >  1:1578780/1‑60 (MQ=255)
                                                                                                                           cGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGt     >  2:1809268/1‑139 (MQ=255)
                                                                                                                           cGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGt     <  2:15834/139‑1 (MQ=255)
                                                                                                                              tCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAg  >  1:515977/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                   
CTGTTAACGCCGGAGCAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAG  >  NC_000913/4106186‑4106450

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