Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,308,627	T→C	11.5%	N100D (AAT→GAT) ‡	nlpI ←	lipoprotein NlpI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,308,627	0	T	C	11.5%	207.7 / 16.3	78	N100D (AAT→GAT) ‡	nlpI	lipoprotein NlpI
Reads supporting (aligned to +/- strand):  ref base T (37/32);  new base C (4/5);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

CGCTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAG  >  NC_000913/3308478‑3308774                                                                                                                                                       |                                                                                                                                                    cgcTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAAtt                                                                                                                                                     <  2:485767/150‑1 (MQ=255)  gcTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTg                                                                                                                                                    <  2:233199/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)   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                                                                                                                                          >  2:1148367/1‑150 (MQ=255)            cACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAg                                                                                                                                          >  1:65951/1‑150 (MQ=255)             aCGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAgg                                                                                                                                         <  2:858843/150‑1 (MQ=255)             aCGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAgg                                                                                                                                         <  2:720013/150‑1 (MQ=255)               gACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCa                                                                                                                                       <  1:1240671/150‑1 (MQ=255)               gACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCa                                                                                                                                       <  1:574883/150‑1 (MQ=255)                 ccgccgTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTGAAGTATTCAATTACTT‑AGGCata                                                                                                                                     <  1:550843/150‑1 (MQ=255)                    ccgTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTGAAGTATTCAATTACTT‑AGGCatatat                                                                                                                                  <  2:1328677/150‑3 (MQ=255)                        aaTATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTc                                                                                                                              <  2:831540/150‑1 (MQ=255)                        aaTATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTc                                                                                                                              <  1:229694/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)                                                                                                                                                       |                                                                                                                                                    CGCTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAG  >  NC_000913/3308478‑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: ‑