Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,308,664	T→C	Q87Q (CAA→CAG)	nlpI ←	lipoprotein NlpI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,308,664	0	T	C	96.9%	214.0 / ‑2.6	65	Q87Q (CAA→CAG)	nlpI	lipoprotein NlpI
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (30/33);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGTTCCATACGTGCCAGAA  >  NC_000913/3308531‑3308797                                                                                                                                      |                                                                                                                                      cgTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGaa                                                                                                                                 >  2:98635/1‑140 (MQ=255)    aGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGaaaaa                                                                                                                              >  1:747893/1‑140 (MQ=255)                gATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCt                                                                                                                  >  2:609640/1‑140 (MQ=255)                gATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCt                                                                                                                  >  2:356937/1‑140 (MQ=255)                gATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCt                                                                                                                  >  1:744711/1‑140 (MQ=255)                gATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCt                                                                                                                  >  1:29682/1‑140 (MQ=255)                gATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCt                                                                                                                  >  1:842516/1‑140 (MQ=255)                gATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTa                                                                                                                 >  2:23106/1‑140 (MQ=255)                  tCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTaa                                                                                                                <  2:747893/140‑1 (MQ=255)                   cAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAAt                                                                                                               <  1:800958/140‑1 (MQ=255)                     aGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGc                                                                                                             >  2:279751/1‑140 (MQ=255)                        tCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCt                                                                                                          <  2:174290/140‑1 (MQ=255)                        tCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCAAAAAATCGTTACGCGCTAATGCCCt                                                                                                          <  1:391671/140‑1 (MQ=255)                          aaGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCa                                                                                                        >  1:446371/1‑140 (MQ=255)                           aGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAg                                                                                                       >  1:71795/1‑140 (MQ=255)                             tACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGAc                                                                                                     <  2:211273/140‑1 (MQ=255)                             tACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGAc                                                                                                     <  1:41479/140‑1 (MQ=255)                                aGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCga                                                                                                  <  1:793658/140‑1 (MQ=255)                                aGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCga                                                                                                  <  2:29682/140‑1 (MQ=255)                                aGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCga                                                                                                  <  1:162840/140‑1 (MQ=255)                                 gAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCgag                                                                                                 <  2:446371/140‑1 (MQ=255)                                 gAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCga                                                                                                  >  1:110826/1‑139 (MQ=255)                                 gAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCga                                                                                                  <  2:110826/139‑1 (MQ=255)                                  aaTCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCgaga                                                                                                <  2:744711/140‑1 (MQ=255)                                  aaTCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCgaga                                                                                                <  2:725286/140‑1 (MQ=255)                                  aaTCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCgaga                                                                                                <  2:71795/140‑1 (MQ=255)                                  aaTCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCgaga                                                                                                >  1:512766/1‑140 (MQ=255)                                   aTCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGAc                                                                                               <  2:812334/140‑1 (MQ=255)                                   aTCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGAc                                                                                               <  2:512766/140‑1 (MQ=255)                                    tCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTa                                                                                             <  1:23106/140‑1 (MQ=255)                                     cAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTa                                                                                             <  1:606936/140‑1 (MQ=255)                                     cAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTa                                                                                             <  1:278870/140‑1 (MQ=255)                                     cAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTa                                                                                             <  1:98635/140‑1 (MQ=255)                                            ttCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATaca                                                                                      >  1:404040/1‑140 (MQ=255)                                                      gcaTCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCcgcg                                                                            <  2:842516/140‑1 (MQ=255)                                                       caTCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTc                                                                                >  1:581291/1‑135 (MQ=255)                                                       caTCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTc                                                                                <  2:581291/135‑1 (MQ=255)                                                        aTCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCt                                                                          <  1:837152/140‑1 (MQ=255)                                                                 gcctgcctGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATaaaa                                                                 >  2:857424/1‑140 (MQ=255)                                                                  cctgcctgCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCt                                                                          <  2:206155/130‑1 (MQ=255)                                                                  cctgcctgCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCt                                                                          >  1:206155/1‑130 (MQ=255)                                                                  cctgcctgCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAg                                                                >  2:509866/1‑140 (MQ=255)                                                                                         tAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTa                                         >  1:452443/1‑140 (MQ=255)                                                                                              aaTTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGcc                                    <  2:381143/140‑1 (MQ=255)                                                                                                       cTTCAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCa                           >  1:427809/1‑140 (MQ=255)                                                                                                          cAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCa                           <  1:287139/137‑1 (MQ=255)                                                                                                          cAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCa                           >  2:287139/1‑137 (MQ=255)                                                                                                          cAGGCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAgg                        >  2:448078/1‑140 (MQ=255)                                                                                                            ggCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGAt                      <  2:452443/140‑1 (MQ=255)                                                                                                             gCATATCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGAtt                     <  1:609640/140‑1 (MQ=255)                                                                                                                 atCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGtt                 >  1:817528/1‑140 (MQ=255)                                                                                                                 atCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGtt                 <  2:817528/140‑1 (MQ=255)                                                                                                                  tCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGa                                 >  1:787867/1‑123 (MQ=255)                                                                                                                  tCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGa                                 <  2:787867/123‑1 (MQ=255)                                                                                                                  tCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGAtt                     <  2:584348/135‑1 (MQ=255)                                                                                                                  tCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGAtt                     >  1:584348/1‑135 (MQ=255)                                                                                                                  tCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGTTc                >  1:858041/1‑140 (MQ=255)                                                                                                                  tCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGTTc                >  2:304104/1‑140 (MQ=255)                                                                                                                  tCCGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGTTc                >  2:845807/1‑140 (MQ=255)                                                                                                                   ccGGTCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGTTcc               <  2:591361/140‑1 (MQ=255)                                                                                                                       tCGGATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGTTCCATAc           >  2:295475/1‑140 (MQ=255)                                                                                                                         ggATTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGTTCCATACGt         <  1:509866/140‑1 (MQ=255)                                                                                                                           aTTGCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGTTCCATACGTGc       >  1:255686/1‑140 (MQ=255)                                                                                                                              gCCAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGTTCCATACGTGCCAg    >  1:522858/1‑140 (MQ=255)                                                                                                                                cAGCGCCTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGTTCCATACGTGCCAGaa  <  2:404040/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                      CGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTAAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGTTCCATACGTGCCAGAA  >  NC_000913/3308531‑3308797

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