Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,587,006	A→C	5.7%	T426P (ACG→CCG) ‡	narQ →	sensory histidine kinase NarQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,587,006	0	A	C	5.7%	261.0 / 7.8	88	T426P (ACG→CCG) ‡	narQ	sensory histidine kinase NarQ
Reads supporting (aligned to +/- strand):  ref base A (47/36);  new base C (3/2);  total (50/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

TTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGC  >  NC_000913/2586861‑2587155                                                                                                                                                   |                                                                                                                                                      ttCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCt                                                                                                                                                   >  2:1074532/1‑150 (MQ=255)   cttacttaCGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTgc                                                                                                                                                 <  1:1380055/150‑1 (MQ=255)      acttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTgcagc                                                                                                                                              <  2:246338/150‑1 (MQ=255)      acttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTgcagc                                                                                                                                              >  2:1305943/1‑150 (MQ=255)               tCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGAtc                                                                                                                                     >  1:614841/1‑150 (MQ=255)                ccAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGCTGCAGCGTCAGGCGGaaag                                                                                                                                    <  2:1157618/150‑4 (MQ=255)                ccAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGAtct                                                                                                                                    >  1:911488/1‑150 (MQ=255)                ccAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTATCCGCCTGCTGCAGCGTCAGGCGGaaag                                                                                                                                    <  1:1520168/150‑4 (MQ=255)                 cAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGAtctc                                                                                                                                   >  1:1270722/1‑150 (MQ=255)                  aGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTcc                                                                                                                                  >  2:528878/1‑150 (MQ=255)                  aGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTAGCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTcc                                                                                                                                  <  2:1281802/150‑1 (MQ=255)                    ttGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTcccc                                                                                                                                <  2:1271157/150‑1 (MQ=255)                     tGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCt                                                                                                                               >  2:1254765/1‑150 (MQ=255)                        cGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCg                                                                                                                            >  1:623099/1‑150 (MQ=255)                        cGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGACGATTTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTcc                                                                                                                             >  2:73305/1‑150 (MQ=255)                             cTGAAGCGTTCGATACCGGAAGATAACGCCACCTCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTg                                                                                                                       >  2:714917/1‑150 (MQ=255)                             cTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTg                                                                                                                       >  2:691946/1‑150 (MQ=255)                              tGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCGGGCGGATCTCCCCTCCGCATTGa                                                                                                                      <  2:847602/150‑1 (MQ=255)                              tGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGa                                                                                                                      >  1:1580038/1‑150 (MQ=255)                                aaGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAgg                                                                                                                    <  2:1636153/150‑1 (MQ=255)                                aaGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAgg                                                                                                                    <  2:1416868/150‑1 (MQ=255)                                      tCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATg                                                                                                              <  1:705292/150‑1 (MQ=255)                                      tCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATg                                                                                                              >  1:1026413/1‑150 (MQ=255)                                      gcGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGGTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATg                                                                                                              <  2:1538617/149‑1 (MQ=255)                                        gATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCt                                                                                                            >  1:30036/1‑150 (MQ=255)                                          tACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTgg                                                                                                          <  2:623099/150‑1 (MQ=255)                                           aCCGGAAGATAACGTCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGa                                                                                                         <  2:1437950/150‑1 (MQ=255)                                           aCCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGa                                                                                                         >  1:44907/1‑150 (MQ=255)                                            ccGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGGGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGAt                                                                                                        <  1:1422666/150‑1 (MQ=255)                                            ccGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGAt                                                                                                        >  2:783851/1‑150 (MQ=255)                                            ccGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGAt                                                                                                        >  1:1358105/1‑150 (MQ=255)                                            ccGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGAt                                                                                                        >  2:1429778/1‑150 (MQ=255)                                            ccGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGAt                                                                                                        >  1:316526/1‑150 (MQ=255)                                             cGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATa                                                                                                       >  1:1335989/1‑150 (MQ=255)                                              ggAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATAc                                                                                                      >  2:980296/1‑150 (MQ=255)                                                 aGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTGATCGGCAGTTACGCGAGCTGTTGAATACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGtt                                                                                                   <  2:536948/150‑1 (MQ=255)                                                    tAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACa                                                                                                >  1:221584/1‑150 (MQ=255)                                                       cGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACaaaa                                                                                             >  2:505714/1‑150 (MQ=255)                                                        gCCACCGCACAAAGTATCATGGCCGA‑TTTTTCACAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAAt                                                                                            >  2:154123/1‑150 (MQ=255)                                                          cACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCa                                                                                          >  1:503236/1‑150 (MQ=255)                                                           aCCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCaa                                                                                         >  2:158527/1‑150 (MQ=255)                                                                     gTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCa                                                                               <  1:1563328/150‑1 (MQ=255)                                                                        tCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAc                                                                            <  1:691946/150‑1 (MQ=255)                                                                        tCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAc                                                                            <  1:1322254/150‑1 (MQ=255)                                                                         cATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACt                                                                           >  1:484346/1‑150 (MQ=255)                                                                            ggCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGAc                                                                        <  1:1585328/150‑1 (MQ=255)                                                                            ggCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGAc                                                                        <  1:1305943/150‑1 (MQ=255)                                                                             gCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTgc                                                                                                                                                 <  2:1452111/76‑1 (MQ=255)                                                                             gCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTgc                                                                                                                                                 >  1:1452111/1‑76 (MQ=255)                                                                             gCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCTCATTGAGGGAAATGCTTGATACGTTACAAAATCAAACCAGCGCCAAACTGAcc                                                                       >  2:103096/1‑150 (MQ=255)                                                                             gCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGAcc                                                                       >  1:1108360/1‑150 (MQ=255)                                                                                   tttttCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTATTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGa                                                                  >  2:1012308/1‑150 (MQ=255)                                                                                    ttttCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGAc                                                                 <  1:487969/150‑1 (MQ=255)                                                                                       tCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGc                                                              <  2:1270722/150‑1 (MQ=255)                                                                                       tCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGc                                                              >  2:291237/1‑150 (MQ=255)                                                                                        cccAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGcc                                                             >  2:1622625/1‑150 (MQ=255)                                                                                             gCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTg                                                        <  2:211310/150‑1 (MQ=255)                                                                                               gTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGcc                                                      >  1:880108/1‑150 (MQ=255)                                                                                               gTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGcc                                                      <  1:678524/150‑1 (MQ=255)                                                                                                    atgatgCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAAccc                                                 <  1:158527/150‑1 (MQ=255)                                                                                                    atgatgCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAAccc                                                 >  2:589515/1‑150 (MQ=255)                                                                                                        tgCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGc                                             <  2:947559/150‑1 (MQ=255)                                                                                                         gCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCa                                            <  2:1358105/150‑1 (MQ=255)                                                                                                          cTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCAc                                           <  1:1254765/150‑1 (MQ=255)                                                                                                                ggCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATg                                     >  1:1604430/1‑150 (MQ=255)                                                                                                                  cAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACTAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATgcg                                   >  1:329135/1‑150 (MQ=255)                                                                                                                  cAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATgcg                                   >  1:484831/1‑150 (MQ=255)                                                                                                                  cAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATgcg                                   <  2:614841/150‑1 (MQ=255)                                                                                                                   aGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATgcgc                                  <  2:503236/150‑1 (MQ=255)                                                                                                                    tttACGCGAGCTGTTGACTACTTTCCTCCTGACGCTGGAGCAGGCGGATCTCCCCTCCGCTTTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAAGCCTTGCACTGGATGCGCa                                 <  2:1026413/149‑1 (MQ=255)                                                                                                                        cgcgAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATg                             >  1:1266994/1‑150 (MQ=255)                                                                                                                        cgcgAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATg                             >  1:532632/1‑150 (MQ=255)                                                                                                                        cgcgAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATg                             >  1:770590/1‑150 (MQ=255)                                                                                                                           gAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAg                          <  1:1074532/150‑1 (MQ=255)                                                                                                                                  tGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCAtt                   >  2:117304/1‑150 (MQ=255)                                                                                                                                  tGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCAtt                   <  2:44907/150‑1 (MQ=255)                                                                                                                                  tGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCAtt                   >  1:181992/1‑150 (MQ=255)                                                                                                                                     ctactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgt                >  1:882694/1‑150 (MQ=255)                                                                                                                                     ctactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgt                <  1:422984/150‑1 (MQ=255)                                                                                                                                     ctactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgt                <  1:1429778/150‑1 (MQ=255)                                                                                                                                      tactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCTCAAATGCACCTGCATttgcg               <  2:221584/150‑3 (MQ=255)                                                                                                                                      ggagatCCGCCTGCTGCAGCGTCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgtt               >  2:1127593/6‑150 (MQ=255)                                                                                                                                       actTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgttg              >  1:1376173/1‑150 (MQ=255)                                                                                                                                        agatCCGCCTGCTGCAGCGTCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGc             >  1:422363/4‑150 (MQ=255)                                                                                                                                           tCCGCCTGCTGCAGCGTCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCATGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCaaa          >  2:1491288/1‑150 (MQ=255)                                                                                                                                            ccGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAat         <  2:1335989/150‑1 (MQ=255)                                                                                                                                                 tGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTc    >  2:217657/1‑150 (MQ=255)                                                                                                                                                  gACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTcg   >  2:785023/1‑150 (MQ=255)                                                                                                                                                   aCGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCTTCTGCAAAGCCTGGCACTGTATGCGCAAATGAAGGTGAATTTGTTGCAAATTTTTctc  <  2:1580038/150‑3 (MQ=255)                                                                                                                                                   |                                                                                                                                                      TTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGA‑TTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGC  >  NC_000913/2586861‑2587155

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