Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,520,084	G→A	32.1%	intergenic (+20/‑178)	curA → / → mcbR	curcumin/dihydrocurcumin reductase, NADPH‑dependent/colanic acid and biofilm gene transcriptional regulator, MqsR‑controlled

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,520,084	0	G	A	32.1%	33.3 / 5.6	43	intergenic (+20/‑178)	curA/mcbR	curcumin/dihydrocurcumin reductase, NADPH‑dependent/colanic acid and biofilm gene transcriptional regulator, MqsR‑controlled
Reads supporting (aligned to +/- strand):  ref base G (9/17);  new base A (8/4);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 8.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAA‑GTTTCAA‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGAC  >  NC_000913/1519863‑1520299                                                                                                                                                                                                                                      |                                                                                                                                                                                                                         gCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCAAAGTGGTGATCCGCGTGGCGGGTGGTGATTAAA‑GTTTTAA‑AGGGGGC‑GTAAGACGCCa                                                                                                                                                                                                                    >  2:347670/1‑228 (MQ=255)         ttattaTCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCAAAGTGGTGATCCGCGTGGCGGGGGATGGTTAAA‑GTTTCAA‑CGCGGGC‑GGAAACCGCCATaaaaaac                                                                                                                                                                                                            >  2:656761/1‑226 (MQ=255)                 gCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑ATTGGGCGAAAGAGGATAAAATCCACTACCGCGAAGAAATTAATGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAA‑GTTTCAA‑AGGCGGG‑GTAAACCGGCATaaaaagaccgggaa                                                                                                                                                                                                     >  2:304876/1‑218 (MQ=255)                 gCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCACAGTGGAGATCAGCGGGGCGGGGGATGCGTAAA‑AATTTAA‑AGGCGGG‑GGAAGCGGCCCCaaaaaaaaagagaaa                                                                                                                                                                                                    >  1:368648/1‑218 (MQ=255)                            tGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCAAAGTGGTGATCCGCGTGGCGGGTGATGTTTAAA‑GTTTCAA‑CGGCGGC‑GTAACGCGCCATAAAAATCCAGGCGAAAAGCTtttttt                                                                                                                                                                                         >  2:12501/1‑228 (MQ=255)                            tGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCAAAGTGGTGATCCGCGTGGCGGGGGATGATTAAA‑GTTTAAA‑AGGGGGG‑GTAAACCGCCCTAAAAATCCCAGCAAAAAACTTAtttt                                                                                                                                                                                         >  2:472045/1‑228 (MQ=255)                                 aCCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCAAAGTGGTGATCCGCGTGGCGGGGGATGTTTAAA‑GTTTCAA‑CGGCAGC‑GGAAGCCGCCATAAAAAGCCTAGTCACAAGGTTAtttttacct                                                                                                                                                                                    >  2:229592/1‑223 (MQ=255)                                 aCCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCAAAGTGGTGATCCGCGTGGCGGGGGATGGTTAAA‑GTTTCCA‑CGGCGGCGGTAAACCGCCCAAAAAATCAGGGCAAAAACTTTTTTTTCAgt                                                                                                                                                                                     >  1:411308/1‑227 (MQ=255)                                 aCCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCAAAGTGGGGATCGGCGTGGCGGGTGAGGGTTAAA‑GTTTAAA‑AGGCGGG‑GTAAGGCGCCCTAAAAATTCAGGCAAAAAACTTATTTTTAGCt                                                                                                                                                                                    >  2:252361/1‑228 (MQ=255)                                      aTCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCAAAGTGGTGATCCGCGTGGCGGGGGGTGATTAAA‑GTTTCAA‑CGGCGGC‑GTAAACCGCCCTAAAAATCCACGCGAAAAGCTTATTTTCAGCTttaat                                                                                                                                                                               >  2:351581/1‑228 (MQ=255)                                            gggattCAGGGGGAAAAGGGGCAAAGGGGTAAAAGAGGAAAAAACCCCCTCCCGGGAAAAAATTCCGTACGGTTTA‑GAGATGGCCC‑CACAGCCGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTC‑CGGCAAAGTGGTGACCCGCGTGGCGGGTGATGATTAAA‑GTTTCAC‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaact                                                                                                                                                                          <  2:344224/224‑1 (MQ=255)                                               tttCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGGTTTCAA‑‑GGCGGC‑GGAAACCCGCACAAAAATCCACGACAAAAGTCTACATTCAACTTTCAATAAATAAcc                                                                                                                                                                      >  2:67376/1‑228 (MQ=255)                                                ttgagggggAGAGGGGGC‑AGGGGGTGAAAGAGGAAAAAACCCCACGGCGAGAAGAAATGGCGGGCGGTTGA‑GATAGGGACC‑AACAGGCGTTTACGGGCCTGGGGAAGTAAAAAAA‑TTC‑CGGCAAATGGGTGCCCCGCGTGGGGGGTAAGGATTAAA‑GTTTCAA‑CGGCGGC‑TTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACca                                                                                                                                                                     <  1:521049/222‑1 (MQ=255)                                                 tCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGAAGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCAAAGTGGTGATCAGCGAGGCGGGGGAGGGTTAAA‑GATTCAA‑CAGCGGC‑GTAAGGCCCCCTAAGAATCCAGGGAAAAAGATTAtttttatggtttattaacaaaacaa                                                                                                                                                                    >  1:266526/1‑207 (MQ=255)                                                         agaTGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCAAAGTGGTGATCCGCGTGGCGGGGGAGGGTTAAA‑GTTTAAA‑CGGCGGG‑GGAAAGCGCCCAAAAAAACCCGGGGAAAAACTTATTTTTACGCttaattaaacaaacaaaaaaattt                                                                                                                                                            >  1:57660/1‑212 (MQ=255)                                                                       ggTGAAAGGGGATAAATGCCACGACCGCAAAGAAATTCGTGACGGTTAA‑GGAAAGGCGC‑AACAGGCTTTTAAGGGCTGGTAGATGGATAAAAA‑CTG‑AGGAAAAGTGTTAATCGGCGTGGGGGTGGATGAGTAAA‑GTTGCAA‑CGGCGGC‑TTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCATCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                                                                                                                                                <  2:320717/227‑1 (MQ=255)                                                                                     aaaCCCACTCCCGCGAAAAATTTTCGGAGGGTTTA‑AAGAAGGCGC‑CACAGGCGTTTATCGGCCTGGTGAAGGGTAAAAC‑CTT‑CGGCAAAGGGGTGACCCGCGTGGCGGTTTATGAAAAAA‑GTTTAAA‑CGGCGGC‑TTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                                                                                                 <  1:550497/228‑1 (MQ=255)                                                                                      acccccccaCCGGGAAAAAATTTCTGACGGTTAA‑GAGAATCCGC‑CACGGACTTTTACGGGCCTGCGAAAGGATAAAAA‑CTT‑GGGCAAGGTGGTCACCCGCGTGGCGGGTGATGATTAAA‑GTTTCAA‑CGGCGGC‑ATAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                                                                                                 <  2:57660/219‑1 (MQ=255)                                                                                      aaTCCGCTACCGAAAAAAAATTATTG‑TGGTTTAGGAGATTCCAC‑GAGAGACGTTTACCGGCGTGCTGAAGTGTAAAAT‑CTA‑CGGCAATGTGGTGACCCGCGTGGCGGGTGATGATAAAA‑TTTTCAA‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                                                                                                 <  1:387509/227‑1 (MQ=255)                                                                                          cACTACCGAGAAAAAATTACTGACGGTTTA‑GGAAAGGCGC‑CACAGACGTTTATGGGCCTGGTAAAGGGTAAAAG‑CTT‑CGGAAAGGTGGTGATCCGGGTGGCGGGTGAGGATAAAA‑GTTTCAA‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGtt                                                                                                                            <  2:36057/228‑1 (MQ=255)                                                                                                  cgAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCAAAGTGGTGATCCGCGTGGCGGGGGATGATTAAA‑GTTTCAA‑CGGCGGC‑GTAAGCCGCCCTAAAAATCCAGGGAAAAAAGTTATTTTTAGCTTTTATTAACTAACCACATCATTGTTTTGCATCTAATAAAACCAAAAGtatggggg                                                                                                                             >  1:36057/1‑212 (MQ=255)                                                                                                    acgaatatACTGACGTTATA‑GAGAAGTGGA‑CCCAGTAGTTGAGCGACCGGCTGAAGGTTAAAAT‑CTG‑CGGGAAAGTGGTACTCCGCTGGGGGGGCGATGATTAAA‑TTTTAAA‑GAGCGGC‑ATAGACCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGattatt                                                                                                                  <  2:672892/221‑1 (MQ=255)                                                                                                    aagaaTTTCCTGATGTTTTA‑AGGAAT‑CGCACACGGGCGTTTACCGGCCTGCTAAAGGGAAAAAA‑CTC‑AGAAAAAGTGGTATACCGCGTGGCGGGTGAAGATTAAA‑GTTTAAA‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAATTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGattatt                                                                                                                  <  2:368648/228‑1 (MQ=255)                                                                                                            aCTGACGTTTAA‑GAGAATGCGC‑CACGGACGTTTAGCGGCAGGTTGAAGTGTAAAAA‑CTT‑CGAAAAGGTGGTAATCCGCGTGGCGGGTAATGATTAAA‑GTTTCAC‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTaa                                                                                                          <  2:158553/228‑1 (MQ=255)                                                                                                                aCGGTGTA‑GAGAAGGCGC‑CACAGACGTTTATGGGCGTGTTAAAGGGTAAAAT‑CTT‑CAGAAAAGTGGTGATCCGGGTGGCGGGTGATGATTAAA‑TTTTCAA‑CGGCGGC‑TAAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATaca                                                                                                      <  2:266526/228‑1 (MQ=255)                                                                                                                         gaaaaTGCGC‑CACAGGCTTTAAGCGGCGTGGTAAAGGATAAAAA‑CTG‑CGGCAAGGTGGTGACCCGGGGGGCGGGGGAAGA‑TAAA‑ATTTCAAGCGGCGGC‑GAAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGa                                                                                              <  2:564405/225‑1 (MQ=255)                                                                                                                                         cggttttacgggcgtggtaaaggTAAAAA‑CTTGCGACAAAGTGGTGACCCGGGTGGGGGGAAAGGTTAAAA‑TTTACAA‑CGGGGCA‑GAAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCa                                                                                 <  2:294927/206‑1 (MQ=255)                                                                                                                                          agagttaggggcaggcgAAAGGTAAAAAGCGT‑CGGCAAGTGGGTTCGCCGGGGGGGGGGGGATGTTTAAA‑TGTTCAC‑GGGCGAC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAAGTCCTTCTGCCATg                                                                               <  2:590425/212‑1 (MQ=255)                                                                                                                                             ttaccgtctggttgaatgtcaaaac‑CTG‑CGGCAAAGTGGTCAACGGCGTGGGTGGTAAGTTGTAAA‑TTTTCAA‑CGGGGTG‑CAAATCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCCCATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATgcagg                                                                           <  2:342320/208‑1 (MQ=255)                                                                                                                                                                       actgcgcaaaaGTGGTGATCCGAGGGGCGGGTATTGTTTAAA‑TTCTCAA‑CGGCGGG‑TAAATCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGa                                                  <  2:301026/220‑1 (MQ=255)                                                                                                                                                                        ttgcGGAAAAGGGGTGACCCGCGGGGCGGTGAATGATTAAG‑TTTTCAA‑GGGCGGC‑GAAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGa                                                  <  2:352664/224‑1 (MQ=255)                                                                                                                                                                        taagaaacagggtgtgACCCGCTTGGCGGTTGATAATTAAA‑GTTCCAA‑GGGCGGC‑TAACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTGTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGa                                                  <  2:581874/214‑1 (MQ=255)                                                                                                                                                                           cGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAA‑GTTTCAA‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTTTGCTGCTAATAACAACAAAGGGGAGTGGTTTATGATTATTTACTTTAATAACCCAGGTGAATTCCTTTTGCCATGCAGGCAGGGGTGGAC‑AGAAAACACCCGAAg                                                >  2:353137/1‑226 (MQ=255)                                                                                                                                                                            ggCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAA‑GTTTCAA‑CGGCGGC‑GTAAGCCGCCCTAAAAATCCAGGGAAAAAAGTTATTTTTAGCTTTTATTTAACAACCACATCCTTTGTTTTCTGCCAATAAAAACCAAAGTGAGTGGGTTATGATTTTTTACTTTTATATACCAGGGGCATTCTTTTTGCCATGTATGTAGGGTTGGGC‑aaaaaaacaccgcaa                                                >  1:411924/1‑216 (MQ=255)                                                                                                                                                                                   tggtgATCCGCGTGGCGGGTGATGATTAAA‑GTTTCAA‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACAAACCACATCATTGTTCTGCTGCTAATAAAAACAAAAGTGAGTGGGTTATTATTATTCACTTTAATACACCAGGGGAATTCCTTTTGCCATGCAGGCAGGGGTGGACAATAAAACAC‑CGAAGAAATTtaaaa                                      >  2:258757/1‑227 (MQ=255)                                                                                                                                                                                                ggggggTAATAATTAAA‑GTTCAAG‑CGGGGGT‑GAAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACggg                          <  2:411308/224‑1 (MQ=255)                                                                                                                                                                                                ggCGGGTTATGATTAAA‑GTTACAA‑CGGGGGG‑GAAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACggg                          <  2:411924/227‑1 (MQ=255)                                                                                                                                                                                                        atttatAAA‑TTTTCAA‑CGGCGGC‑ATAAACCGCCATAAAAATCCAGGCAAAAAGCTTATTTACAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGaa                 <  1:353137/223‑1 (MQ=255)                                                                                                                                                                                                             tAAA‑TTTCCAA‑AGGCGGC‑GTAGTCCGGCATAAATATCCAGGCAAAAAGCTTATTTTCATCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACAt             <  2:87492/227‑1 (MQ=255)                                                                                                                                                                                                                 attttttcagcggctatttAAACCGCCATAAAAATCCAGGCAAAAAACTTATTTACAGCTTTAAATAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCTATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGt           <  1:302889/210‑1 (MQ=255)                                                                                                                                                                                                                      cAA‑CGGCGGC‑GTAAACCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCGTTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTTTTCACTTTATAACACCAGGGGAATTCCCTTCGCCATtcag                                                                            >  1:38884/1‑157 (MQ=255)                                                                                                                                                                                                                         ctgcgggc‑tTAAACCGCCATAAAAATCCAGGAAAAAAGTTTATTTTCAGCTTTAAGTAACTAACCACATCATTGTTCTTCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGAc  <  1:252361/223‑1 (MQ=255)                                                                                                                                                                                                                          acctggtaatAACCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAAATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGAc  <  2:282201/219‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                         GCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTT‑CGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAA‑GTTTCAA‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGAC  >  NC_000913/1519863‑1520299

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