Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,805,475	T→C	100%	N159S (AAT→AGT)	waaP ←	kinase that phosphorylates core heptose of lipopolysaccharide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,805,475	0	T	C	99.0%	335.6 / ‑5.8	99	N159S (AAT→AGT)	waaP	kinase that phosphorylates core heptose of lipopolysaccharide
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (44/53);  total (45/54)
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

TGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCAATC  >  NC_000913/3805376‑3805569                                                                                                    |                                                                                               ttagCCAGATGCAGATCGATAACTGAAATTTTTAATGCATCTTGCCGGGCAGTAGATGGCAAATGAGGCAAAAAGTGGCAAATGTAACAGTCGCGGTGAt                                                                                                <  2:2200841/97‑1 (MQ=255)    gCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTa                                                                                             >  1:1055149/1‑100 (MQ=255)    gCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTa                                                                                             <  1:1517308/100‑1 (MQ=255)    gCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTa                                                                                             >  1:1505477/1‑100 (MQ=255)     cccGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTAt                                                                                            >  1:446927/1‑100 (MQ=255)      ccGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATc                                                                                           >  1:755072/1‑100 (MQ=255)      ccGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATc                                                                                           >  2:173828/1‑100 (MQ=255)        gATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATccc                                                                                         <  1:1431372/100‑1 (MQ=255)          tGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTg                                                                                       <  1:2141664/100‑1 (MQ=255)            cAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCGAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCa                                                                                     >  2:1576570/1‑100 (MQ=255)            cAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCa                                                                                     >  1:723574/1‑100 (MQ=255)            cAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCa                                                                                     >  2:1259930/1‑100 (MQ=255)              gATCGATAACTGAAATTTTTAATTCATCTTCCCGGTCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAAAAGTCGCGGTGACGTATCCCTGGAGt                                                                                   >  2:1474114/1‑100 (MQ=255)              gATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGt                                                                                   <  2:2454345/100‑1 (MQ=255)               aTCGATAATTGAATTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAACAAAAAGTGGCAAATGTCACACTCGCGCTGATTTATCCCTGCat                                                                                    <  2:884163/98‑2 (MQ=255)               aTCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTa                                                                                  >  2:1352924/1‑100 (MQ=255)                 cGATAACTGAAATTTTTAATTCATCTTCCCGGCAAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATg                                                                                >  1:895448/1‑100 (MQ=255)                  gATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGCCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGc                                                                               >  2:1946538/1‑100 (MQ=255)                     aaCTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCAtt                                                                            <  2:1849953/100‑1 (MQ=255)                      aCTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCAttt                                                                           <  2:37564/100‑1 (MQ=255)                       ntGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCAtttt                                                                          >  2:1715874/2‑100 (MQ=255)                       cTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCAtttt                                                                          >  1:698452/1‑100 (MQ=255)                       cTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCAtttt                                                                          >  2:1745957/1‑100 (MQ=255)                        tGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTa                                                                         <  1:266779/100‑1 (MQ=255)                         gAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTAc                                                                        >  2:936909/1‑100 (MQ=255)                          aaaTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACg                                                                       <  2:1246631/100‑1 (MQ=255)                           aaTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGc                                                                      >  2:1453782/1‑100 (MQ=255)                            aTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCa                                                                     <  1:861350/100‑1 (MQ=255)                             tttttAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCAc                                                                    >  1:893801/1‑100 (MQ=255)                              ttttAATACACCTTCCCGGGCAGTAAATGTCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCAc                                                                    <  2:2430901/99‑1 (MQ=255)                              gtacgacccATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCAcc                                                                   <  2:814093/92‑1 (MQ=255)                                   aTTCATCTACCCGGCCCGTAAATGGCAAATGAAGCAAAAAGAGGCAACTGTAACAGTCGCGGAGACTTATCCCTGCAGTATGCATTTTACTCACCATAGt                                                              >  2:1109130/1‑100 (MQ=255)                                      cATGTTCCCGGCCAGAAAATGGCAAATGAAGCAAAAAGGAGCAAATGTAACAGTCGCGGCGACTTATCCCTGCCGTATGCATTTTACGCACCATAGtt                                                             >  2:330351/1‑98 (MQ=255)                                      cATCTTCCCGGCCAGTAAATGGCAAATTAAGTAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTg                                                            <  2:149842/99‑1 (MQ=255)                                      cATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGc                                                           <  2:2517292/100‑1 (MQ=255)                                      cATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGc                                                           <  2:2039653/100‑1 (MQ=255)                                      cATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGc                                                           <  1:561733/100‑1 (MQ=255)                                       aTCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACTATAGTTGCt                                                          >  2:1488409/1‑100 (MQ=255)                                       aTCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCt                                                          <  2:953578/100‑1 (MQ=255)                                        tCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTa                                                         <  1:1441391/100‑1 (MQ=255)                                         cTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTac                                                        <  1:2263720/100‑1 (MQ=255)                                          ttCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTaca                                                       >  1:1832438/1‑100 (MQ=255)                                            cccGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACg                                                     <  1:561239/100‑1 (MQ=255)                                            cccGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACg                                                     >  1:1275663/1‑100 (MQ=255)                                             ccGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGt                                                    >  1:2400965/1‑100 (MQ=255)                                               ggCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGc                                                  <  1:245751/100‑1 (MQ=255)                                               ggCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGc                                                  <  1:2253944/100‑1 (MQ=255)                                               ggCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGc                                                  <  2:1575635/100‑1 (MQ=255)                                                 ncAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGa                                                >  2:547450/2‑100 (MQ=255)                                                  cAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGAt                                               <  1:1240284/100‑1 (MQ=255)                                                  cAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGAt                                               <  1:1248263/100‑1 (MQ=255)                                                   aGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATc                                              <  1:161038/100‑1 (MQ=255)                                                   aGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATc                                              >  2:1483077/1‑100 (MQ=255)                                                    gTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCg                                                 >  2:1592837/1‑96 (MQ=255)                                                    gTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCa                                             <  1:1004591/100‑1 (MQ=255)                                                     tAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAg                                            >  1:1744803/1‑100 (MQ=255)                                                     tAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAg                                            <  1:1657592/100‑1 (MQ=255)                                                      aaaTGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGc                                           <  1:1634778/100‑1 (MQ=255)                                                      aaaTAGCAAATGAAGCAAAAAGTGGCAAATGTAAGAGACGCGGTGACTTATCCCTGGATTATGCAGTTTACGCACCATAGTTGCTACACGTGCGATCAGc                                           <  2:358072/100‑1 (MQ=255)                                                       aaTGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCa                                          >  1:52895/1‑100 (MQ=255)                                                       aaTGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCa                                          <  1:1714940/100‑1 (MQ=255)                                                        aTGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCAt                                         <  1:936917/100‑1 (MQ=255)                                                         tGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATa                                        >  1:1094199/1‑100 (MQ=255)                                                            cAAATGAGGCAAAAAGTGGCAAATGTAACAGTCGCGTTGATGAATCACTCCAGAATTCATTTTATGCACCATAGATGCTACACGTGCGATCAGCATatg                                      >  2:1876066/1‑97 (MQ=255)                                                            cAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGc                                     >  1:818177/1‑100 (MQ=255)                                                            cAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGc                                     >  1:2140886/1‑100 (MQ=255)                                                             aaaTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCt                                    <  1:1083346/100‑1 (MQ=255)                                                             aaaTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCt                                    >  1:615464/1‑100 (MQ=255)                                                             aaaTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCt                                    <  2:633318/100‑1 (MQ=255)                                                               aTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTa                                  >  2:409664/1‑100 (MQ=255)                                                               aTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTa                                  <  2:1310939/100‑1 (MQ=255)                                                                tGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTaa                                 >  2:2210809/1‑100 (MQ=255)                                                                 gAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAac                                <  1:870225/100‑1 (MQ=255)                                                                   aGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAacac                              >  2:717495/1‑100 (MQ=255)                                                                   aGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAacac                              <  1:513889/100‑1 (MQ=255)                                                                    gCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACg                             <  1:1364135/100‑1 (MQ=255)                                                                    gCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACg                             >  2:2112925/1‑100 (MQ=255)                                                                      aaaaaGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGta                           >  1:1906665/1‑100 (MQ=255)                                                                      aaaaaGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGta                           <  2:2471209/100‑1 (MQ=255)                                                                      aaaaaGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGta                           <  2:1266250/100‑1 (MQ=255)                                                                       aaaaGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGtat                          >  1:400524/1‑100 (MQ=255)                                                                        aaaGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGtata                         <  1:288834/100‑1 (MQ=255)                                                                               cAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCg                  >  1:1781848/1‑100 (MQ=255)                                                                                aaaTGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCgg                 <  1:1910204/100‑1 (MQ=255)                                                                                 aaTGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCggg                <  1:381884/100‑1 (MQ=255)                                                                                   tGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGtt              >  1:1174799/1‑100 (MQ=255)                                                                                   tGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGtt              <  1:2048058/100‑1 (MQ=255)                                                                                      aaCAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGAc           <  1:1702977/100‑1 (MQ=255)                                                                                      aaCAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGAc           <  2:1887679/100‑1 (MQ=255)                                                                                        cAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTg         <  2:1361384/100‑1 (MQ=255)                                                                                         aGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGc        >  1:402942/1‑100 (MQ=255)                                                                                         aGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCAACATAGTTGCTACACGTGGGATCAGCATACGCTTAACAAGTATATCAGGCGGGTTGACTGc        >  2:2169108/1‑100 (MQ=255)                                                                                          gTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACt          <  2:1663323/97‑1 (MQ=255)                                                                                          gTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGcc       <  1:1389184/100‑1 (MQ=255)                                                                                           tCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGccc      <  1:2506867/100‑1 (MQ=255)                                                                                           tCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGccc      <  2:1505466/100‑1 (MQ=255)                                                                                            cgcgGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGcccn     <  2:1275655/100‑2 (MQ=255)                                                                                             gcgGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCa     >  1:1936546/1‑99 (MQ=255)                                                                                               ggTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCAATc  <  2:446923/100‑1 (MQ=255)                                                                                                 tGACTTAGCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCAATc  <  2:504052/98‑1 (MQ=255)                                                                                                    |                                                                                               TGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCAATC  >  NC_000913/3805376‑3805569

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