Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,308,625	A→T	11.5%	N100K (AAT→AAA) ‡	nlpI ←	lipoprotein NlpI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,308,625	0	A	T	11.5%	214.8 / 16.9	78	N100K (AAT→AAA) ‡	nlpI	lipoprotein NlpI
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base T (4/5);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

TGCGCTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAG  >  NC_000913/3308476‑3308774                                                                                                                                                       |                                                                                                                                                      tGCGCTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTaa                                                                                                                                                       >  1:341945/1‑150 (MQ=255)   cgcTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAAtt                                                                                                                                                     <  2:485767/150‑1 (MQ=255)    gcTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTg                                                                                                                                                    >  2:890761/1‑150 (MQ=255)    gcTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTg                                                                                                                                                    <  2:889566/150‑1 (MQ=255)    gcTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTg                                                                                                                                                    <  1:1265798/150‑1 (MQ=255)    gcTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTg                                                                                                                                                    <  2:233199/150‑1 (MQ=255)     cTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGa                                                                                                                                                   >  2:1029311/1‑150 (MQ=255)           tGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACtt                                                                                                                                             >  2:950428/1‑150 (MQ=255)              cACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAg                                                                                                                                          >  1:21482/1‑150 (MQ=255)              cACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAg                                                                                                                                          >  1:65951/1‑150 (MQ=255)              cACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAg                                                                                                                                          >  2:1148367/1‑150 (MQ=255)               aCGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAgg                                                                                                                                         <  2:858843/150‑1 (MQ=255)               aCGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAgg                                                                                                                                         <  2:720013/150‑1 (MQ=255)                 gACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCa                                                                                                                                       <  1:574883/150‑1 (MQ=255)                 gACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCa                                                                                                                                       <  1:1240671/150‑1 (MQ=255)                   ccgccgTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTGAAGTATTCAATTACTT‑AGGCata                                                                                                                                     <  1:550843/150‑1 (MQ=255)                      ccgTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTGAAGTATTCAATTACTT‑AGGCatatat                                                                                                                                  <  2:1328677/150‑3 (MQ=255)                          aaTATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTc                                                                                                                              <  1:229694/150‑1 (MQ=255)                          aaTATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTc                                                                                                                              <  2:831540/150‑1 (MQ=255)                           atatAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTGAAGTATTCAATTACTT‑AGGCatatatttaac                                                                                                                             <  1:1273598/150‑8 (MQ=38)                              taATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGAt                                                                                                                          <  2:653981/150‑1 (MQ=255)                               aaTGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGAtt                                                                                                                         >  2:537723/1‑150 (MQ=255)                                 tGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGc                                                                                                                       >  1:209536/1‑150 (MQ=255)                                  gCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGTATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGcc                                                                                                                      <  2:1149856/150‑1 (MQ=255)                                     ttCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTTGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAgc                                                                                                                   <  2:1236224/149‑1 (MQ=255)                                       cccGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAgcgc                                                                                                                 >  2:1207549/1‑150 (MQ=255)                                           cgATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGc                                                                                                             >  2:423414/1‑150 (MQ=255)                                              ttCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGaa                                                                                                          >  2:787509/1‑150 (MQ=255)                                                      cgcgTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGtt                                                                                                  <  1:1033988/150‑1 (MQ=255)                                                      cgcgTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGtt                                                                                                  <  1:890761/150‑1 (MQ=255)                                                       gagTAGTTCTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTa                                                                                                 >  2:772051/3‑150 (MQ=255)                                                        cgTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTAc                                                                                                <  1:1197631/150‑1 (MQ=255)                                                        cgTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTAc                                                                                                <  2:201655/150‑1 (MQ=255)                                                             ttGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCt                                                                                           <  2:181539/150‑1 (MQ=255)                                                              tGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTa                                                                                          >  2:750235/1‑150 (MQ=255)                                                                      ggATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTGAAGTATTCAATTACTT‑AGGCata                                                                                                                                     <  2:1176308/99‑1 (MQ=25)                                                                      ggATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTGAAGTATTCAATTACTT‑AGGCata                                                                                                                                     >  1:1176308/1‑99 (MQ=25)                                                                          cAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGAc                                                                              >  2:101300/1‑150 (MQ=255)                                                                            aGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCt                                                                                           <  2:318781/135‑1 (MQ=255)                                                                            aGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCt                                                                                           >  1:318781/1‑135 (MQ=255)                                                                            aGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTAAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCg                                                                            >  2:70413/1‑150 (MQ=255)                                                                             gCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAGAAATCGTTACGCGCTAATGCCCTCAGACCga                                                                           >  2:720041/1‑150 (MQ=255)                                                                              cTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCgag                                                                          >  1:830199/1‑150 (MQ=255)                                                                              cTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCgag                                                                          <  2:21482/150‑1 (MQ=255)                                                                                cAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGAc                                                                        <  1:1148367/150‑1 (MQ=255)                                                                                cAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGAc                                                                        <  1:1153076/150‑1 (MQ=255)                                                                                 aaGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCga                                                                           <  1:577227/146‑1 (MQ=255)                                                                                 aaGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCga                                                                           >  2:577227/1‑146 (MQ=255)                                                                                  aGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTa                                                                      >  1:1391230/1‑150 (MQ=255)                                                                                    tACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATc                                                                    >  2:930378/1‑150 (MQ=255)                                                                                     aCAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCa                                                                   <  1:1029311/150‑1 (MQ=255)                                                                                     aCAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCa                                                                   >  2:1172567/1‑150 (MQ=255)                                                                                        gAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATac                                                                <  1:500671/150‑1 (MQ=255)                                                                                        gAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATac                                                                >  1:106886/1‑150 (MQ=255)                                                                                              aaCGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCACACCGAGACTATCATACAACACt                                                          >  2:1247689/1‑150 (MQ=255)                                                                                                  cTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCcgc                                                      <  2:127114/150‑1 (MQ=255)                                                                                                  cTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCcgc                                                      >  2:64101/1‑150 (MQ=255)                                                                                                  cTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCcgc                                                      <  1:1207549/150‑1 (MQ=255)                                                                                                    tCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCgcgc                                                    >  1:93406/1‑150 (MQ=255)                                                                                                     cATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCt                                                   <  1:1172567/150‑1 (MQ=255)                                                                                                           cagcaTCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCAtata                                             <  2:65951/150‑1 (MQ=255)                                                                                                              caTCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATaaaa                                          >  1:44436/1‑150 (MQ=255)                                                                                                              caTCAAAATTGCCTGCCCGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCa                                                 <  2:958308/143‑1 (MQ=255)                                                                                                              caTCAAAATTGCCTGCCCGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCa                                                 >  1:958308/1‑143 (MQ=255)                                                                                                                tCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGc                                        >  1:1363698/1‑150 (MQ=255)                                                                                                                    aaTTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCtgtg                                    <  1:423414/150‑1 (MQ=255)                                                                                                                      ttGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTgcg                                  >  2:506137/1‑150 (MQ=255)                                                                                                                           tgcctgcGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTc                             >  2:987357/1‑150 (MQ=255)                                                                                                                            gcctgcGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCg                            <  1:950428/150‑1 (MQ=255)                                                                                                                                  gTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCg                      >  2:1275505/1‑150 (MQ=255)                                                                                                                                  gTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCg                      >  2:1019511/1‑150 (MQ=255)                                                                                                                                     ccggatatGCCTGAAGTATTCAATTACTT‑AGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAg                                         <  2:600035/124‑1 (MQ=37)                                                                                                                                     ccggatatGCCTGAAGTATTCAATTACTT‑AGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAg                                         >  1:600035/5‑128 (MQ=37)                                                                                                                                      cggatatGCCTGAAGTATTCAATTACTT‑AGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTa                  >  2:1090129/4‑150 (MQ=38)                                                                                                                                      cggatatGCCTGAAGTATTCAATTACTT‑AGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTa                  >  2:1048397/4‑150 (MQ=38)                                                                                                                                      aaTATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTa                  <  1:772051/150‑1 (MQ=255)                                                                                                                                            tGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGccc            >  2:1391701/1‑150 (MQ=255)                                                                                                                                               cT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGAc         >  2:1343525/1‑150 (MQ=255)                                                                                                                                                     tAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCa    >  2:276070/1‑150 (MQ=255)                                                                                                                                                       aTTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAg  >  1:1273462/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                      TGCGCTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAG  >  NC_000913/3308476‑3308774

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: ‑