Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,308,630	A→T	12.3%	F99I (TTC→ATC) ‡	nlpI ←	lipoprotein NlpI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,308,630	0	A	T	12.3%	196.6 / 15.4	73	F99I (TTC→ATC) ‡	nlpI	lipoprotein NlpI
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base T (4/5);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

TGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATT  >  NC_000913/3308486‑3308778                                                                                                                                                  |                                                                                                                                                     tGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACtt                                                                                                                                                 >  2:950428/1‑150 (MQ=255)    cACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAg                                                                                                                                              >  1:65951/1‑150 (MQ=255)    cACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAg                                                                                                                                              >  2:1148367/1‑150 (MQ=255)    cACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAg                                                                                                                                              >  1:21482/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:890761/150‑1 (MQ=255)                                            cgcgTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGtt                                                                                                      <  1:1033988/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                                                                              <  2:21482/150‑1 (MQ=255)                                                                    cTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCgag                                                                              >  1:830199/1‑150 (MQ=255)                                                                      cAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGAc                                                                            <  1:1148367/150‑1 (MQ=255)                                                                      cAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGAc                                                                            <  1:1153076/150‑1 (MQ=255)                                                                       aaGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCga                                                                               >  2:577227/1‑146 (MQ=255)                                                                       aaGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCga                                                                               <  1:577227/146‑1 (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                                                     >  1:958308/1‑143 (MQ=255)                                                                                                    caTCAAAATTGCCTGCCCGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCa                                                     <  2:958308/143‑1 (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)                                                                                                                                                 aaTACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGAtt  >  1:140009/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATT  >  NC_000913/3308486‑3308778

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