Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,343,821	0	C	A	42.4%	87.6 / 32.5	59	A117S (GCT→TCT)	yfaL	adhesin
Reads supporting (aligned to +/- strand):  ref base C (28/6);  new base A (0/25);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 4.54e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCAAAAATACGGCCCCACCGCTGGCGTTATTGGTCATATCC  >  NC_000913/2343672‑2343966                                                                                                                                                       |                                                                                                                                                    gATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGc                                                                                                                                                     >  1:1001545/1‑150 (MQ=255)       aTTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAatat                                                                                                                                              <  1:1108576/150‑1 (MQ=255)             aaaCATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGG‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑cccc                                                                                                                                       <  2:795919/150‑1 (MQ=255)                cATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTCCCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGAAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCccgccg                                                                                                                                   <  2:44654/151‑1 (MQ=255)                      gTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACAAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGttatt                                                                                                                              >  2:398498/1‑150 (MQ=255)                               acGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCt                                                                                                                    >  1:228069/1‑151 (MQ=255)                                 gTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGt                                                                                                                  >  1:529861/1‑151 (MQ=255)                                   aaaCCGACGGCACCATTTTCGTTAGAACAAAAAAAATAGATTCCGCCCCAATACCCCCCTGCGTTGTAACCAAAAAAAACAACATCATA‑AAGATTAGGTTTTGATTTTTTTTAAGA‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGt                                                                                                                  <  2:344095/149‑1 (MQ=255)                                             cACCGGTTTCGTTAGTACAAAAAAAAAAGATCGCCCCACCATGCCCGCCTGCGACTTTACCAAAAAAAATAACAAAAGT‑ATTATTCTGTTTTGATTTTTCTTTAAA‑AAATAGGG‑CCCCGCCGTTATTATATTCTCCTGTAACGgtgt                                                                                                          <  2:228020/147‑1 (MQ=255)                                              acaattatttttAGAACAAAAAAAATAGATCCCCCCCCCATCCCCGCTTGCGACTTTACGAAAAAAAATAACATCAGA‑AATTTTCGGCTTTGATTTTTCTTTAAA‑AAATATGC‑CCCCGCCGTTATTATATTCTCCTGTAACGgtgt                                                                                                          <  1:313876/136‑1 (MQ=255)                                                cAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGc                                                                                                   >  1:515670/1‑151 (MQ=255)                                                      cGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAAc                                                                                               >  1:804752/1‑149 (MQ=255)                                                          aGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTc                                                                                         >  1:467522/1‑151 (MQ=255)                                                             aCCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCAtt                                                                                      >  1:93102/1‑151 (MQ=255)                                                                agaagaATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTcc                                                                                    >  2:694769/1‑150 (MQ=255)                                                                  aaaaaTAGATCGCCCCCCAATAGCCGCCTGCGACTTTACCGAAAAAAAAAACATCAGA‑AAGATCAGGTTTTGTTTTTTCTTAAAA‑AAATATGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTcc                                                                                    <  1:712669/145‑1 (MQ=255)                                                                     aaTAGTTGGCCCCCCCATCGCCGCCGGCGACTTTACCGAAAAAAAAAACATAAGA‑AAGATTCGGCGTTGTTTTTTTTTTAAA‑AAATAGGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTc                                                                              <  1:544556/151‑1 (MQ=255)                                                                     aaTAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTc                                                                              >  1:610657/1‑151 (MQ=255)                                                                     aaTAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTc                                                                              >  1:798282/1‑151 (MQ=255)                                                                     aaTAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTc                                                                              >  2:590361/1‑151 (MQ=255)                                                                          aTTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGTCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCAtttt                                                                         <  1:154577/151‑1 (MQ=255)                                                                           ttGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCATGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCAtttt                                                                         >  2:388002/1‑150 (MQ=255)                                                                            gccgccccaaacgccgccTGCGATGTTACGAAAAAAAATAACACTAAT‑ATATTATAGCTTTGTTTTTTTTATATA‑AAATAGCG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCAt                                                                            <  2:153248/134‑1 (MQ=255)                                                                              cgcCACAATGGCCGCCTGCGACTTTACAAAAAAAAATAACAATAGT‑AAGATTCGGTGTTGTTTTTTCTTTAAA‑AAATTTGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCt                                                                       <  2:265647/149‑1 (MQ=255)                                                                              ccccccaatcccccccTGCGACTTTACCAAAAAAAAAAACATAAAA‑AATATTCAGTTTTGTTTTTTTTTAAAA‑AAATGGGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTc                                                                        <  1:293411/134‑1 (MQ=255)                                                                               ccccccatccccgccgGCGACTTTCCCAAAAAAAAAAACATAAAA‑AAAATTCGGCGTTGTTTTTTTTTAAAA‑AAATGTGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTc                                                                        <  1:827490/137‑1 (MQ=255)                                                                               cccaccacgcccgccTGCGTCGTTACCAAAAAAAATAACACCAGA‑AATATTCAGGTTTTTTTTTTCTTTAAA‑AAATAGGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGt                                                                    <  2:474513/141‑1 (MQ=255)                                                                                     aTAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGt                                                              >  1:476591/1‑151 (MQ=255)                                                                                         ccgccgGCGACTTTACCAAAAAAAAAAACATCAGA‑AAGATTCGGGTTTGTTTTTTTTTAAGA‑AAATAGGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAg                                                               <  1:385455/146‑1 (MQ=255)                                                                                         ccgccTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGaa                                                           >  2:241489/1‑150 (MQ=255)                                                                                          cgccTGCGACTTTACGAAAAAAAATAACATCAGT‑AAGATTCAGCGTTGTTTTTTTTTTAAA‑AAATATGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGa                                                            <  1:1102256/148‑1 (MQ=255)                                                                                           gccTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAAt                                                        >  2:873673/1‑151 (MQ=255)                                                                                            ccTGTGACTTTACAGAAAAAAAAAACATCAAT‑AATTTTCTGTTTTGTTTTTTTTTTAAA‑AAATACGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTaa                                                                  <  2:14264/140‑1 (MQ=255)                                                                                            ccTGCGACTTTCCCGAAAAAAAAAACATAAGT‑AAGTTTCAGTTTTGTTTTTTTTTAAGA‑AAATATGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGaa                                                           <  1:753257/147‑1 (MQ=255)                                                                                            ccTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAAtt                                                       >  2:1051696/1‑151 (MQ=255)                                                                                              tccgcctttACC‑GAAAAAAAAACCACAGTAAACATTCGGGGTTG‑TTTTTTTTTAGA‑AAAATCGGCCCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTg                                                             <  1:735201/136‑1 (MQ=255)                                                                                               gcgtcttcaccaaaaaaaaaaACAACAAA‑AAGTTTCAGTTTTGTTTTTTCTTGAAA‑AAATATGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCg                                                    <  1:707308/138‑1 (MQ=255)                                                                                                cGACGTTACCGGAAAAAATAACATCAGA‑ACGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGc                                                   >  1:738875/1‑151 (MQ=255)                                                                                                          ggAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGc                                         >  1:337018/1‑151 (MQ=255)                                                                                                           gAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGc                                         >  1:324969/1‑150 (MQ=255)                                                                                                           aaaaaaaaaaaATAAGT‑AAGTTTCAGCGTTGTTTTTTTTTAAGA‑AAATAGGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCttg                                             <  1:1018257/145‑1 (MQ=255)                                                                                                           aaaaaaaaaACATAAGT‑AAAATTCGGGTTTTTTTTTTCTTAAAA‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTcc                                                     <  1:504175/137‑1 (MQ=255)                                                                                                           aaaaaaaTAACATAAGT‑AAGTTTCAGTTTTGATTTTTTTTAAAA‑AAATATGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCa                                        <  2:1042990/150‑1 (MQ=255)                                                                                                           aaaaaaaTAAAATAAGT‑AAGATTCGGTGTTGTTTTTTCTTAAGA‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCa                                        <  1:32146/150‑1 (MQ=255)                                                                                                            aaaaaaaaACATCATA‑ATTTTTAAGTGTTGTTTTTTTTTTAAA‑AAATGGGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCaa                                       <  2:140596/151‑1 (MQ=255)                                                                                                              aaaaTAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCaaa                                      >  1:509469/1‑150 (MQ=255)                                                                                                               aaaTAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCaaaa                                     >  1:1125396/1‑150 (MQ=255)                                                                                                                 aaaacaacAGT‑AATATTCAGTTTTGTGTTTTTTTAACA‑AAATATGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCAAAAATa                                  <  2:17365/149‑1 (MQ=255)                                                                                                                   aaCATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCAAAAATACg                                >  2:330900/1‑151 (MQ=255)                                                                                                                    aacaccatgaataTTCGGCTTTGATTTTTTTTTAGA‑AAAATTGC‑CTCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCAAAAATa                                  <  1:449748/137‑1 (MQ=255)                                                                                                                              aatgttcatggttgtttttttttACCAAAATACGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTTCAAAAATCCGGCACCACCg                       <  2:655994/136‑1 (MQ=255)                                                                                                                              aGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCAAAAATACGGCCCCACCGc                      >  1:1103359/1‑151 (MQ=255)                                                                                                                                  tctgtttttgtttttttttAAA‑ATATATGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCAAAAATACGGCCCCAcc                        <  2:926761/134‑1 (MQ=255)                                                                                                                                    aGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCAAAAATACGGCCCCACCGCTGGCGt                >  1:478802/1‑151 (MQ=255)                                                                                                                                          gAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCAAAAATACGGCCCCACCGCTGGCGTTATTgg          >  2:82060/1‑151 (MQ=255)                                                                                                                                           gtttttttttAACAAAATAGGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCAAAAATACGGCCCCACCGCTGGCg                 <  1:1051961/142‑1 (MQ=255)                                                                                                                                             ttttCTTAAAA‑AAATACGC‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTGTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCAAAAATACGGCCCCACCGc                      <  2:429578/136‑1 (MQ=255)                                                                                                                                               ttCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCAAAAATACGGCCCCACCGCTGGCGTTATTGGTCat      >  2:796752/1‑150 (MQ=255)                                                                                                                                                 cTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCAAAAATACGGCCCCACCGCTGGCGTTATTGGTCATATc   >  1:508240/1‑151 (MQ=255)                                                                                                                                                  tttAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCAAAAATACGGCCCCACCGCTGGCGTTATTGGTCATATcc  >  2:500680/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                    GATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT‑AAGATTCAGCGTTGAGTTTTCTTTAGC‑AAATATGG‑CCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCTTGTTGCAAAAATACGGCCCCACCGCTGGCGTTATTGGTCATATCC  >  NC_000913/2343672‑2343966

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