Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,796,432	T→C	L152L (TTG→CTG)	waaC →	ADP‑heptose:LPS heptosyltransferase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,796,432	0	T	C	97.8%	154.7 / NA	46	L152L (TTG→CTG)	waaC	ADP‑heptose:LPS heptosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (18/27);  minor base G (0/1);  total (18/28)
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

CCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGG‑ACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGAC  >  NC_000913/3796187‑3796663                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        ccATAAACGCTGAACGCAACGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGc                                                                                                                                                                                                                                    >  1:318111/1‑251 (MQ=255)        agctgAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAgc                                                                                                                                                                                                                              >  1:174758/5‑250 (MQ=255)        agctgAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAgc                                                                                                                                                                                                                              >  1:75312/5‑250 (MQ=255)              aCGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAgcaaaccgc                                                                                                                                                                                                                       >  2:347232/1‑251 (MQ=255)                 cAAAGCGAGTCGTAAAGCGCTACAAGCAGAGAACTATGACGCAGTTAGCGATGCTCAGGGGCTGGTAGAAAGCGGGGCGCTGGTGACGCGGCTGGCGCATGGCGTAACGCATGGATTGT‑ACTGGCAAACCGCACGGGAGCCGTGAGCCAGCCTGTGTGGTAATCGGAAGCATCATGTTGCAAATCAGGAGCACGCAGGAGAACGCACGCGCGAACGGTTGGCCAAAAGTGTGGGCTATTgcaatcggcaaa                                                                                                                                                                                                                    <  2:23009/251‑1 (MQ=255)                       gTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCACACCCAgg                                                                                                                                                                                                              <  1:401801/251‑1 (MQ=255)                             ggAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCGCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGAtt                                                                                                                                                                                                        <  2:157048/250‑1 (MQ=255)                                   gcTACAAGCGGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAACAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCCAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCt                                                                                                                                                                                                   >  2:181613/1‑250 (MQ=255)                                    cTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTAt                                                                                                                                                                                                 >  1:335658/1‑251 (MQ=255)                                      aCAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTAttg                                                                                                                                                                                               <  1:180347/251‑3 (MQ=255)                                         aGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGcac                                                                                                                                                                                            <  2:415183/251‑1 (MQ=255)                                          gCAGAGAACTATGACGCAGTTAGCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGcaca                                                                                                                                                                                           <  2:376643/251‑1 (MQ=255)                                              agaACTATGACGCAGTTATCGATGCACAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGGAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAGCCCAGGGCGATTATGCTATTGCACAGCAt                                                                                                                                                                                       <  2:416528/251‑1 (MQ=255)                                                    aTGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCAtttcct                                                                                                                                                                                  <  2:287234/250‑4 (MQ=255)                                                      gACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGGGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACTCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGa                                                                                                                                                                                <  1:187813/250‑1 (MQ=255)                                                               aTCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGAAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGcc                                                                                                                                                                      >  1:51684/1‑251 (MQ=255)                                                                            ggCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCAGGGCGTAAAGCATGGATTGG‑ACTGGCAACGCGTTTGCGACCCGTTAGCCTGCCTGTTTTACAATCGTAAGGATCATATTGCAAAACAGCGGCACGCCGTAGGACGCTCCAGCGAACTGTTTTCCAAAAGTCTGGGCTATATCAATCCGTAGACCCAGTGCGATTGGGCTATTTCTCAGCAGTTCCTGTGGAACCTGCGTACGGATGCTGGc                                                                                                                                                         <  2:193809/251‑1 (MQ=255)                                                                             gCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGGTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATGATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCg                                                                                                                                                        <  1:414378/251‑1 (MQ=255)                                                                                      aaGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCg                                                                                                                                               <  1:330714/251‑1 (MQ=255)                                                                                          gcggcgCTGATGACGCGTCTGGCGCATGGCGTACAGCATGGATTGG‑ACTGGCAAACCGCTCGCGCACCGTTAGCCAGCCTGTTTTACACTCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGACCTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCAGTAtt                                                                                                                                           >  2:359712/1‑251 (MQ=255)                                                                                           cggcgCTGGTGACGCGTCTGGCGCATGGCTGAAAGAATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTAttt                                                                                                                                          <  2:196702/251‑1 (MQ=255)                                                                                               gcTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATttcttc                                                                                                                                      <  1:283845/251‑1 (MQ=255)                                                                                                           tCTGGCGCATGGCGAAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGAAAGCCTGGTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACACGCGAACTGATTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGATATTGCACAGCATTTCCTGCCGAACCTGCCTACAGAAGCTGGCGAATCTGCCGAATTTCTTCATGCGAAGAccc                                                                                                                          >  1:334565/1‑251 (MQ=255)                                                                                                              ggCGCATGGCGGAAAGCATGG‑ATGGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCATGTTTTACTATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACTGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATTCCGTATTTGGTTATGCGATGACCCg                                                                                                                         <  2:366772/249‑1 (MQ=255)                                                                                                                 gcagggCGTAAAGCATGGAATGG‑ACGGGCAAACCTCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCATCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGtgatg                                                                                                                    <  2:68302/247‑1 (MQ=255)                                                                                                                        gTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAac                                                                                                              >  1:291820/1‑250 (MQ=255)                                                                                                                         tAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCAtt                                                                                                                                                                                      <  2:203011/177‑1 (MQ=255)                                                                                                                         tAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCAtt                                                                                                                                                                                      >  1:203011/1‑177 (MQ=255)                                                                                                                          aaaGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACt                                                                                                           <  2:93760/251‑1 (MQ=255)                                                                                                                               atggattgg‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACGCTGGCCg                                                                                                      <  2:270425/245‑1 (MQ=255)                                                                                                                                           tGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGCCCCGTGATTATAAACACTGGCCGGAAGACCACTg                                                                                           >  2:417850/1‑251 (MQ=255)                                                                                                                                               aaaccccaagcgAACCGTTAGCCAGCCTGTTTTACAACCGTAAGCATCATATTGCAACACCGCAGCACGCCGTAGAACGCACCCGCGAACTGCTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCCGACGAACCTGACTTCCGATGCTGGCGAATATGACGTATTTCTTCATGCGACGACCCGTGATGATAAACAATTGCCGGAAAAaccattgggg                                                                                       >  2:434817/10‑243 (MQ=255)                                                                                                                                                    gCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCTAGAatt                                                                                  <  2:69861/251‑1 (MQ=255)                                                                                                                                                                     cagccTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACAGTGGCCGGAAGAGCACTGGCGAGAATTGATTGGTTTACTGGCTg                                                                 <  1:359712/251‑1 (MQ=255)                                                                                                                                                                           gTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGACAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACTAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTTATGATAAACACTGGCCGGAATAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTc                                                             <  2:395261/249‑1 (MQ=255)                                                                                                                                                                           gTTTTAAAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCAAACGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCAGTATTTCTTCATGCGACGACCCGTGATGATAAACCCTGGCCGGACGAACAATGGAGAGAATTGATTGGTTTAATGGCTGCTTc                                                             >  2:416025/1‑249 (MQ=255)                                                                                                                                                                                        aaGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAAGACGGATTa                                               <  1:98133/250‑1 (MQ=255)                                                                                                                                                                                         aGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTaaa                                             <  1:181613/251‑1 (MQ=255)                                                                                                                                                                                               atatTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCg                                       <  2:318111/251‑1 (MQ=255)                                                                                                                                                                                                             cagcaCGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGcgccgc                          >  1:23469/1‑250 (MQ=255)                                                                                                                                                                                                              agcaCGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGAGTGGTTTACTGGCTGATGCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCAt                        <  2:291820/251‑1 (MQ=255)                                                                                                                                                                                                                  cGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAgg                    >  2:329839/1‑251 (MQ=255)                                                                                                                                                                                                                        aGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGCGCATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTggg                                   >  1:174101/1‑230 (MQ=255)                                                                                                                                                                                                                        aGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACCGATTAAACTTCCGTggg                                   <  2:174101/230‑1 (MQ=255)                                                                                                                                                                                                                               aCCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGa        >  2:283797/1‑250 (MQ=255)                                                                                                                                                                                                                                    cgAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGAc  <  1:329839/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        CCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGG‑ACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGAC  >  NC_000913/3796187‑3796663

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