Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,042,604	0	A	C	45.9%	104.8 / 15.8	61	V182G (GTG→GGG)	yqfA	hemolysin III family HylIII inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (26/7);  new base C (0/28);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 3.60e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGT  >  NC_000913/3042455‑3042745                                                                                                                                                       |                                                                                                                                               ggCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCa                                                                                                                                                <  1:61378/150‑1 (MQ=255) ggCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCAc                                                                                                                                               <  1:569210/151‑1 (MQ=255)    gcgTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCAcccc                                                                                                                                            >  2:18716/1‑151 (MQ=255)    gcgTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCAcccc                                                                                                                                            <  2:683153/151‑1 (MQ=255)          tATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGc                                                                                                                                       <  2:373670/149‑1 (MQ=255)             cAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATa                                                                                                                                  <  1:837280/151‑1 (MQ=255)              aGGCCTACAAAGGCATCCCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATa                                                                                                                                  <  1:511899/150‑1 (MQ=255)              aGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATa                                                                                                                                  >  2:15897/1‑150 (MQ=255)                gCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAAc                                                                                                                               >  2:271034/1‑151 (MQ=255)                    aCAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCAc                                                                                                                            <  2:253701/150‑1 (MQ=255)                      aaaGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCAcgcc                                                                                                                         >  1:500807/1‑151 (MQ=255)                       aaGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCAcgccg                                                                                                                        >  2:813474/1‑151 (MQ=255)                           cATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACgccgcc                                                                                                                      >  1:235322/1‑149 (MQ=255)                                ccATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGcc                                                                                                               >  2:69802/1‑151 (MQ=255)                                      cGCCAGCCTAATATCAAAATAGATCACCAGAAGTTGGCACCCACACCCGCCGATCCCGAGCCCGTGCGATATGGCCTGATTGTATGGAAT‑GCGTTTGCAGACGTAAAAAATCCCCCC‑GCGCGAATAACCCACGCCGCCTACCGCCAGt                                                                                                            <  1:578768/148‑1 (MQ=255)                                           gCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAAc                                                                                                    >  1:605663/1‑151 (MQ=255)                                            cccAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAAc                                                                                                    >  2:474891/1‑150 (MQ=255)                                              cAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCt                                                                                                 >  2:94372/1‑151 (MQ=255)                                                atatACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGc                                                                                               >  1:786000/1‑151 (MQ=255)                                                 tataCAAATAGCCCGCAAGAAAGTGGCACACACTACCGCCGAGCTCGAATCCGTGCCAGATGGCATGGTTGTATGGAAT‑GGGTTTGCAGACGTAAAAAATCACCCC‑GAGCGAATAACCCGCGCCGCCTACCGCCAGTAAGGTAAc                                                                                                    <  1:418532/145‑1 (MQ=255)                                                     cAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCCCCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCccgc                                                                                           <  2:234201/150‑1 (MQ=255)                                                       aaTAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCcgccgc                                                                                        >  1:731871/1‑151 (MQ=255)                                                        aTAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCcgcg                                                                                       >  2:563126/1‑151 (MQ=255)                                                              cGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCtt                                                                                  >  1:389223/1‑150 (MQ=255)                                                              cGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCtt                                                                                  >  2:534283/1‑150 (MQ=255)                                                                ccAGAAGGGGCCACCCACTCCCGCCGAGCACGAAGCCGGCCCAGATGGCTGGGTTGTAGGGAAG‑GCGTTGGCAGACGTAAAAATTCCCCCC‑GGGCGAAAAAACCCCCCCCCTTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTa                                                                                 <  1:168914/149‑1 (MQ=255)                                                                  aGAAAGGGCCACACTCTACCGCCAGGCACGAAGCCGTGCCAGAGGGCATGTTTGTATGGAAT‑GCGTTTGCAGACGAAAAAAACCCCCCC‑GAGCAAAAAACCCCCCCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACt                                                                              <  2:401956/150‑1 (MQ=255)                                                                    aaggggccacaccCTCCCCCGGACCACGAGCCGGTGCAAGAGGCCTGGGTTTTATGGAAG‑GCGTTTGAAGCCGGAAAAAACCCCCCC‑GAGCAAAAAACCCCCCCCCCTTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAg                                                                                     <  1:358223/134‑1 (MQ=255)                                                                     aagggccacacccccccgccgAGCACGAAGCCGTGCAAGATGGCATGGTTGTATGGAAT‑GCTTTTGCAGACGTAAAAAATCCCCCC‑GAGCAAATAACCCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCa                                                                          <  1:177155/136‑1 (MQ=255)                                                                      gggggCCCACCCTCCCGCCGAGCACAAACCCGTGCAAGAGGGCATGGTTGTAGGGAAT‑GCGTTTGAAGACGTAAAAAACCCCCCC‑GAGCAAAAAACCCCCCCCGCTTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCAt                                                                         <  1:354111/148‑1 (MQ=255)                                                                             cacacTCCCGCCGAGCACGAAGCCGTGCCAGATGGCATGTTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCCCCCC‑GAGCGAATAAACCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATaa                                                                  <  2:677961/151‑1 (MQ=255)                                                                                acTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAAtt                                                               >  2:916835/1‑151 (MQ=255)                                                                                    ccgccgAGCCCGAAGCCGTGCCAGAGGGCATGTTTGTAGGGAAT‑GCGTTTGCAGACGTAGAAAACCCCCCC‑GAGCAAAAAACCCACGCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAAt                                                                <  2:700139/146‑1 (MQ=255)                                                                                    ccgccgAGCACGAAGCCGTGCCAGATGGCATGTTTTTATGAAAT‑GCTTTTGCAGACGTAGAAAATCCCCCC‑GAGCAAATAACCCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAAtt                                                               <  2:808127/147‑1 (MQ=255)                                                                                       ccgAGCCCGAAGCCGTGCCAGATGGCTTGGTTTTATGGAAT‑GCTTTTGCAGACGTAAAAAACCCCCCC‑GAGCAAATAACCCCCGCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTaccacc                                                         <  2:445824/150‑1 (MQ=255)                                                                                             acaacgCCGAGCCAGATGGCAGGTTTATATGAAAT‑GCTTTGGCAGACGAAAAAAACCCCCCC‑GGGGAAAAAACCCCCCCCCCCTCCCCCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCg                                                     <  1:383153/143‑1 (MQ=255)                                                                                              cGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAc                                                   >  2:231125/1‑149 (MQ=255)                                                                                                aggcggTGCCAGAGGGCTTGGTTTTAGGAAAT‑GCTTTTGCAGACGTAAAAAACCCCCCC‑GAGCAAAAAACCCCCCCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGc                                                      <  1:758369/139‑1 (MQ=255)                                                                                                aggCCGGGCCAGATGGCTGGTTTGTA‑GGAATGGCGTTTGCAGGTGTAGAAAATACCCCC‑GGGCGAATAACCCCCGCCCCCTCCCGCCAGTAAGGAAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAc                                                   <  2:777865/145‑1 (MQ=255)                                                                                                  gCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagcc                                               >  1:613952/1‑149 (MQ=255)                                                                                                   cggtcccAGATGCCATGGTTGTATGAAAT‑GCGTTGGAAGCCGTAGAAAACCCCCCC‑GACCGAAAAAACCACGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagcca                                              <  2:203905/144‑1 (MQ=255)                                                                                                   ccGTGCAAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGCCGAAAAAAACCCCCCC‑GAGCAAAAAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAg                                                       <  1:456861/140‑1 (MQ=255)                                                                                                      tGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGccc                                          >  2:778509/1‑150 (MQ=255)                                                                                                             tggctgggTTTTATGGAAT‑GCTTTTGCGGCCTTAGAAAATCCCCCC‑GGGCAAATAACCCCCCCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCa                                   <  1:906783/144‑1 (MQ=255)                                                                                                             tGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGcc                                    >  1:400784/1‑149 (MQ=255)                                                                                                              ggCATGGTTGTATGAAAT‑GCTTTTGAAAACTAAAAAAACCCCCCC‑GAGCAAAAAACCCCCCCCCCCTCCCGCCAGAAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCa                                         <  1:793845/143‑1 (MQ=255)                                                                                                              ggCATGGTTGTAGGGAAT‑GCGTTTGCAGACGTAAAAAACCCCCCC‑GAGCGAATAACCCACGCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGc                                     <  2:503239/147‑1 (MQ=255)                                                                                                                cATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGAt                                >  2:221572/1‑150 (MQ=255)                                                                                                                  tGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAgg                             >  2:401758/1‑151 (MQ=255)                                                                                                                   gtttttacgaaaT‑GTGTTTGCAGATGTAGAAAATCCCCCC‑GCGCGAATAACCCCCCCCGCTTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGAt                                <  2:31484/137‑1 (MQ=255)                                                                                                                   ggTTGTATGAAAT‑GGTTTTGCAGACGTAAAATTCCCCCCC‑GGGCAAATAACCCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATa                               <  2:547443/148‑1 (MQ=255)                                                                                                                   ggTTGTATGAAAT‑GCGTTTGCAGACGTAAAAAACCCCCCC‑GAGCGAATAACCCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAg                              <  2:882819/149‑1 (MQ=255)                                                                                                                   ggTTGTAGGAAAT‑GCGTTTGAAGACGTAAAAAACCCCCCC‑GGGCAAAAAACCCCCGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAg                              <  1:467361/149‑1 (MQ=255)                                                                                                                       gTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAcc                        >  1:568657/1‑151 (MQ=255)                                                                                                                        tatgataT‑GCGTTTCCAGCGGTAAAAAACCCCCCC‑GAGCGAAAAACCCCCCCCCCCACCCGCCAGTAAGGTAACGCCGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGa                                 <  1:512897/135‑1 (MQ=255)                                                                                                                        tATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAcc                        >  2:96679/1‑150 (MQ=255)                                                                                                                           ggAAT‑GCGTTTGCAGACGTAAAAAACCCCCCC‑GAGCGAAAAACCCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCagag                    <  2:219817/151‑1 (MQ=255)                                                                                                                           ggAAT‑GCGTTTGAAGCCGTAAAAAACCCCCCC‑GAGCAAATAACCCCCCCCCCTTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAcc                        <  2:634478/147‑1 (MQ=255)                                                                                                                                 gCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATaa                >  2:845726/1‑150 (MQ=255)                                                                                                                                  cGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATaaa               >  2:911168/1‑150 (MQ=255)                                                                                                                                               aaaaaaCCCCCCC‑GGGCGAATAACCCCCGCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGt  <  2:67688/148‑1 (MQ=255)                                                                                                                                                  aaaTCCCCCCGGAGCAAAAAACCCCCGCCGCTTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGa        <  1:49784/142‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                               GGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAT‑GCGTTTGCAGACGTAGAAAATCACCCC‑GAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGT  >  NC_000913/3042455‑3042745

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