Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,303,672	0	C	A	40.7%	92.8 / 29.2	55	intergenic (‑175/‑233)	napF/eco	ferredoxin‑type protein, role in electron transfer to periplasmic nitrate reductase NapA/ecotin, a serine protease inhibitor
Reads supporting (aligned to +/- strand):  ref base C (4/28);  new base A (22/0);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.92e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAATTGAGAGCGACCAAATAAACCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTCTGCTGGTTTGCTCTCTCATTAGAATTTAACACTAAAAGAGCAGGTAAAATTGTCTGAATGTTCTTTAAGTTATTCATAAA  >  NC_000913/2303525‑2303821                                                                                                                                                             |                                                                                                                                                        gAAAATTGAGAGCGACCAAATAAACCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAAc                                                                                                                                                         <  2:430465/148‑1 (MQ=255) gAAAATTGAGAGCGACCAAATAAACCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAAc                                                                                                                                                         <  1:694505/148‑1 (MQ=255)     aTTGAGAGCGACCAAATAAACCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCttt                                                                                                                                                  <  1:624775/151‑1 (MQ=255)        gagaGCGACCAAATAAACCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTg                                                                                                                                                <  1:415056/150‑1 (MQ=255)           agCGACCAAATAAACCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑G‑AAAAAAAGGCGTTTTCt                                                                                                                                              >  2:546795/1‑149 (MQ=255)                   aaTAAACCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑aa                                                                                                                                  <  2:58703/151‑1 (MQ=255)                   aaTAAACCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑aa                                                                                                                                  <  1:1077295/151‑1 (MQ=255)                        aCCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑GA‑GAG‑G‑G‑GG‑G‑AAAAAAGGGGtttttttcccttttaaa                                                                                                                                    >  1:1028129/1‑134 (MQ=255)                           gCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑G‑G‑C‑GGAC‑AAAAAAGGGGtttgct                                                                                                                                               >  1:969388/1‑130 (MQ=255)                            cATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑G‑GG‑G‑AAAAAAAGCCGTGTGCTTACA‑TTTC‑taa                                                                                                                                 >  1:53395/1‑143 (MQ=255)                              taattaatAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑G‑G‑G‑GG‑G‑AAAAAAGGGGGTTTTTTTACA‑TTTT‑AAAttggt                                                                                                                            >  1:1016242/1‑142 (MQ=255)                               aattaatAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑GA‑G‑C‑G‑G‑GG‑G‑AAAAAAGGGCGTTTTCTCCCA‑TTTT‑CAATGGTTTTaaa                                                                                                                       >  2:649493/1‑150 (MQ=255)                               aattaatAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑G‑G‑G‑GG‑C‑GAAAACGGGGTTTTTTTTACC‑TTTT‑AAATGGGTATtaa                                                                                                                       >  2:495092/1‑150 (MQ=255)                               aattaatAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATaaa                                                                                                                       <  1:509276/150‑1 (MQ=255)                                   aataaGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAatatat                                                                                                                  <  2:1078046/151‑1 (MQ=255)                                      aaGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGGAAA‑G‑G‑G‑G‑GG‑A‑AAAAAAGGGGTTTTG‑TTACA‑TTTT‑AAAATGGTATaaaaaaaat                                                                                                                 >  1:734016/1‑143 (MQ=255)                                      aaGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTta                                                                                                               <  1:933566/151‑1 (MQ=255)                                      aaGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑A‑GAG‑GG‑G‑AAAAAAGGGGTTTT‑TTTTCA‑TTTT‑AAATAGTGATAAAAATAatt                                                                                                                >  2:142262/1‑150 (MQ=255)                                      aaGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGACGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑A‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAACATATTta                                                                                                               <  1:535370/151‑1 (MQ=255)                                        gCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTtata                                                                                                             <  1:545801/151‑1 (MQ=255)                                              tttATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑G‑G‑TGGG‑G‑AAAAAAGGGC‑TTTGCTTCCA‑TTTT‑CAAAGGTGATaaaaa                                                                                                                     >  1:38599/1‑135 (MQ=255)                                               ttATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCga                                                                                                       <  1:1071479/150‑1 (MQ=255)                                                tataGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCgat                                                                                                      <  2:745158/150‑1 (MQ=255)                                                 ataGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCgatt                                                                                                     <  2:1013917/150‑1 (MQ=255)                                                    gCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCgattgatt                                                                                                 <  2:86918/151‑1 (MQ=255)                                                       gCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑GAG‑G‑G‑GG‑G‑AAAAAAGGGCGTTT‑TTTAAA‑TTTT‑AAAAAGTTATAAAAATATTTATATAACgaat                                                                                                     >  1:496431/1‑142 (MQ=255)                                                           aGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCaga                                                                                          <  2:600071/151‑1 (MQ=255)                                                             atatTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑G‑GG‑C‑CAAAAAAGGCGTTTTTTCACC‑TTTT‑AAATAGTTATAAATATATTTTTTTAATgattgttt                                                                                                 >  2:599140/1‑142 (MQ=255)                                                               atTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGCGGGG‑AA‑C‑G‑G‑G‑GG‑G‑AAAAAAGGGGTTTT‑CTTCCA‑TTTT‑AAAAGGGTATTAATATAATTTTTTAAAgattgattttacccaga                                                                                        >  2:606443/1‑140 (MQ=255)                                                               atTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑A‑G‑G‑G‑GG‑G‑AAAAAAAGGGTTTTTCTCCCA‑TTTT‑CAAAGGTTTTTAAAAAATTTTTATAACgattgatttaac                                                                                             >  2:371385/1‑140 (MQ=255)                                                                  aaaGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑G‑GG‑G‑AAAAAAAGGGGTTTTTTTTCA‑TTTT‑AAATTGTTTTAAAAATATTTTTATAACGATTTATTTAACACAGAAAttt                                                                                   >  2:593828/1‑151 (MQ=255)                                                                   aaGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑A‑G‑G‑GG‑C‑AAAAAAAGGCCTTTTTTCTCA‑TTTT‑TAATAGGTATAAATATATTTATTTATAGATTTATTTACCAAAGAAAtttt                                                                                  >  2:920701/1‑150 (MQ=255)                                                                       aTGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTctgct                                                                              <  1:1068549/151‑1 (MQ=255)                                                                        tgtgTCAAAGATGCATGCCCCGAGCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTctgct                                                                              <  1:292248/150‑1 (MQ=255)                                                                         gtgtCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTctgctg                                                                             <  1:581532/150‑1 (MQ=255)                                                                          tgtCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑G‑AAAAAAAGGGGTTTGTTCACA‑TTTT‑AAAAGGGTATAAATATATTTTTATATCGATTTATTTACCCGAGATAttttttcttg                                                                            >  1:545218/1‑143 (MQ=255)                                                                            tCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑G‑G‑G‑GG‑CGAAAAAAGGGATTTT‑TTCACA‑ATTC‑AAATATTTTTAAAAATATTTTTATATCGATTTGTTCACCATAGAAATTTTTTTTGGtt                                                                          >  1:1101898/1‑150 (MQ=255)                                                                             cAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTCTGCTGGttt                                                                         <  2:916477/150‑1 (MQ=255)                                                                                      aTACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAGACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTCTGCTGGTTTGCTCTCTCAt                                                               <  2:872574/151‑1 (MQ=255)                                                                                      aTACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTCTGCTGGTTTGCTCTCTCAt                                                               <  2:961722/151‑1 (MQ=255)                                                                                         ccccGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTCTGCTGGTTTGCTCTCTCATTAg                                                            <  2:694861/151‑1 (MQ=255)                                                                                         ccccGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑G‑GG‑G‑AAAAAAAGGGTTTTTTTCTCA‑TTTT‑AAATAGGTATAAAAATATTTTTATATAGATTGATTTTCCACAAATATTTTTTTTTGTTTTTTCTCTCAtt                                                              >  2:1002358/1‑149 (MQ=255)                                                                                           ccGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGGAAA‑G‑G‑G‑G‑GG‑C‑AAAAA‑GGGGTTTTTTTCCCC‑TTTC‑AAATAGGTTTAAATAAATTTTTATAAAGCTTTTTTTACCAAACTTATTTTTGCTGGgtttgttttttttttta                                                           >  1:263121/1‑136 (MQ=255)                                                                                                cGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑GA‑G‑CGG‑G‑GG‑G‑AAAAAAGGGCTTTT‑TTCACA‑TTTT‑AAAATGTTATAAAAAAAATTTTATAACGATTGATTTTCCAGAGAAATTTTTTTTTGTTTTCTTTTTCATTATAAatt                                                       >  1:10927/1‑149 (MQ=255)                                                                                                  gggTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AAGG‑G‑G‑G‑GG‑G‑AAAAAAGGGGTTTT‑TTCACA‑ATTT‑AAAAGGGTATAAAAATATTTTTTTAACGATTTGTTTTCCCGAGAAAATTTTTTTGGTTTTTTCTCTCTTTTTAAttt                                                       >  1:271120/1‑147 (MQ=255)                                                                                                      aaaaCCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTCTGCTGGTTTGCTCTCTCATTAGAATTTAACACTaa                                               <  1:934539/151‑1 (MQ=255)                                                                                                       aaaCCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTCTGCTGGTTTGCTCTCTCATTAGAATTTAACACTaaa                                              <  2:587931/151‑1 (MQ=255)                                                                                                        aaCCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑A‑C‑G‑G‑GG‑G‑CAAAACGGGCTTTTTTTTACATTTTC‑AAAT‑GTTATAAAAATATTTTTATAATGATTTATTTACCAAAGATATTTTTTCTTGATTTCTCTCTCCATACAATTTAACAATaaa                                              >  2:267824/1‑150 (MQ=255)                                                                                                           cTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑G‑GG‑G‑CAAAAAAGGCTTTTTTTCACC‑TTTCAAAATGGTT‑TAAAAATATTTTTATAACGCTTGATTTTCCCAAGAAATTTTTTCTGGGTTTCTCTCTCCTTTTAATTTAAAACTaaaaaa                                           >  2:930886/1‑148 (MQ=255)                                                                                                                gCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑G‑GG‑G‑AAAAAAAGGCCTTTTGTCACA‑TTTT‑TAATGGGTATTAAAAAATTTTTTTAAAAATTTATTTAACCGAAATATTTTTTTTGGTTTTTTTTCTCATTTTAATTTTAAACAaaaaaa                                           >  2:1011625/1‑143 (MQ=255)                                                                                                                                        aa‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTCTGCTGGTTTGCTCTCTCATTAGAATTTAACACTAAAAGAGCAGGTAAAATTGTCTGAATGTTCTTTa               <  1:657146/151‑1 (MQ=255)                                                                                                                                         a‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTCTGCTGGTTTGCTCTCTCATTAGAATTTAACACTAAAAGAGCAGGTAAAATTGTCTGAATGTTCTTTaa              <  1:476314/151‑1 (MQ=255)                                                                                                                                           g‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTCTGCTGGTTTGCTCTCTCATTAGAATTTAACACTAAAAGAGCAGGTAAAATTGTCTGAATGTTCTTTAAg             <  2:2212/151‑1 (MQ=255)                                                                                                                                             c‑G‑T‑GG‑C‑AAGAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTCTGCTGGTTTGCTCTCTCATTAGAATTTAACACTAAAAGAGCAGGTAAAATTGTCTGAATGTTCTTTAAGt            <  1:195430/151‑1 (MQ=255)                                                                                                                                             c‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTCTGCTGGTTTGCTCTCTCATTAGAATTTAACACTAAAAGAGCAGGTAAAATTGTCTGAATGTTCTTTAAGt            >  2:178599/1‑151 (MQ=255)                                                                                                                                                        aaaaaCGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTCTGCTGGTTTGCTCTCTCATTAGAATTTAACACTAAAAGAGCAGGTAAAATTGTCTGAATGTTCTTTAAGTTATTc       >  2:870050/1‑150 (MQ=255)                                                                                                                                                            aCGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTCTGCTGGTTTGCTCTCTCATTAGAATTTAACACTAAAAGAGCAGGTAAAATTGTCTGAATGTTCTTTAAGTTATTCATaaa  <  1:593735/151‑1 (MQ=255)                                                                                                                                                             |                                                                                                                                                        GAAAATTGAGAGCGACCAAATAAACCGCATAATTAATAAGCCATTTTTATAGCCGCTAAGATATTAAAGGATGTGTCAAAGATGCATACCCCGATCGGGGTAAAACCTGTAGCAGGATCAAAAGAGTGGC‑GGG‑AA‑G‑C‑G‑T‑GG‑C‑AAAAACGGGCTTTTGCTCACA‑TTTC‑AAATGGTTATAAATATATTTATATAGCGATTGATTCACCAGAGATATTTCTGCTGGTTTGCTCTCTCATTAGAATTTAACACTAAAAGAGCAGGTAAAATTGTCTGAATGTTCTTTAAGTTATTCATAAA  >  NC_000913/2303525‑2303821

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