Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,190,000	T→C	D40D (GAT→GAC)	yqiH →	putative fimbrial chaperone YqiH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,190,000	0	T	C	100.0%	165.5 / NA	46	D40D (GAT→GAC)	yqiH	putative fimbrial chaperone YqiH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/20);  total (26/20)

ATAACAAACTTGCTTTACGCTGGGGTGACAAATCTTGTTTTATTCAACCACCCAATAGCAGCAACTTAACCACCGGAACGGCTATTTTACCGTGTATTAGCCAAAATTAATTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCA  >  NC_000913/3189756‑3190238                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               ggactcgtgaaaaactggcgctgcttgaacaaagaagcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctctccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATAGCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCAGCCCGGATCGTACCCGTACCATAATGAATTAGTCAGACAAAGGc                                                                                                                                                                                                                                          <  2:49015‑M2/141‑1 (MQ=255)   actcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGa                                                                                                                                                                                                                                        <  1:427112‑M2/143‑1 (MQ=255)                   cgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACt                                                                                                                                                                                                                        >  1:192153‑M2/93‑251 (MQ=255)                       gcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATc                                                                                                                                                                                                                    <  2:15544‑M2/163‑1 (MQ=255)                         ttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCa                                                                                                                                                                                                                   >  1:335093‑M2/87‑250 (MQ=255)                                               ggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTACAAACCTGCCCTATCTgg                                                                                                                                                                                            >  2:129187‑M2/65‑251 (MQ=255)                                                   aatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCa                                                                                                                                                                                        <  1:389951‑M2/191‑1 (MQ=255)                                                      cgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAAc                                                                                                                                                                                     <  2:335093‑M2/194‑1 (MQ=255)                                                        cagtactgagatcaaaaagcagccagccccagagactctccccccaacgctgtaaTTCAGGACGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAACGCGAAGCATTAAAACTTACTAATAAAAGTAAAAACCTGCCCTATCTGGACCAATCAt                                                                                                                                                                                   >  2:51536‑M2/56‑251 (MQ=255)                                                           tactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATgg                                                                                                                                                                                 >  2:179862‑M2/53‑250 (MQ=255)                                                                     aaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACCAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATa                                                                                                                                                                       >  2:64105‑M2/43‑250 (MQ=255)                                                                       aaagcagccagatccagagactctcccaccaacgctgtatTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAGCTTACTAATCAAAGTAAAAAACTGCCCAATCTGGCCCAAACATGGATAGAAGATACg                                                                                                                                                                     >  2:327061‑M2/41‑250 (MQ=255)                                                                                     cagagactctcccccgaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAAACTGCCCTATCTGGCCCAAACTTGGATAGAAGATACGAAAGGCAATAAATcg                                                                                                                                                      >  2:449098‑M2/27‑251 (MQ=255)                                                                                      agagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTa                                                                                                                                                                                                 >  2:216374‑M2/26‑207 (MQ=255)                                                                                       gagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTAGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCAGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGAAATAAATcgcg                                                                                                                                                    >  2:416699‑M2/25‑251 (MQ=255)                                                                                        ggacgatccccccaacgctgtaattCAGGATGAAACTAATGCGCTACTGAAATTCAAAGAAAATCCTTGCAGCAGGAGTATTACTCTCCCGTATGAGCATGAATGCCTGGGGTGCCATCAGCCCGGATCATACCCGTAACATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCATAGTTAAAACCTGCCCTTTCGCGCCCAAACGTGGATAGAAGATACGAAATGCAATAAATCGCGt                                                                                                                                                   <  2:210160/231‑1 (MQ=255)                                                                                                cccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTAt                                                                                                                                           >  1:452552‑M2/16‑251 (MQ=255)                                                                                                               ttCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTtcctcc                                                                                                                            >  2:61096/1‑251 (MQ=255)                                                                                                                         aaCTAATGCGCTACTTAAATACAAAGAATATCCTTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGACGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACCTGGATAGAAGATACGAAAGGCAATAAATCGCGTGCCTTTATTGTCACTGTTCCTCCACTGGTGCGt                                                                                                                  >  1:288072/1‑251 (MQ=255)                                                                                                                                           aTACAACGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGAATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGAGAGCa                                                                                                 >  1:23330/1‑250 (MQ=255)                                                                                                                                                aaGAATATCATGGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTCGGGTGCCATCATCCCAGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTTCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGCGCGTTTAAATCCCAGTGAGCAAATTc                                                                                            >  1:56607/1‑250 (MQ=255)                                                                                                                                                         aTTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAAGCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAAt                                                                                  <  1:219974/251‑1 (MQ=255)                                                                                                                                                               gcagGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCAATAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGTCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATAGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAAc                                                                            <  2:468675/251‑1 (MQ=255)                                                                                                                                                                       aTTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGaa                                                                     <  1:28894/250‑1 (MQ=255)                                                                                                                                                                        ttACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGaaa                                                                    >  1:320165/1‑250 (MQ=255)                                                                                                                                                                         tACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCAGGATCGTACCCGTATCATAATGAATGAGTCAGACATAGGCGGAGCCTTAAAACTTACTAATCAAAGGAAAAACCGGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGAGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGaaaa                                                                   >  1:340455/1‑250 (MQ=255)                                                                                                                                                                               ccTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAAGCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCt                                                            <  2:192153/251‑1 (MQ=255)                                                                                                                                                                                  gTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCa                                                          >  1:339989/1‑250 (MQ=255)                                                                                                                                                                                    aTGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGc                                                                                                                                                                                                     <  2:81429/109‑1 (MQ=255)                                                                                                                                                                                    aTGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGc                                                                                                                                                                                                     >  1:81429/1‑109 (MQ=255)                                                                                                                                                                                    aTGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAAc                                                        <  1:179862/250‑1 (MQ=255)                                                                                                                                                                                       agcagcATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTcc                                                    >  1:197801/1‑251 (MQ=255)                                                                                                                                                                                       agcagcATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTcc                                                    >  1:343681/1‑251 (MQ=255)                                                                                                                                                                                                   tGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATc                                                                                                                                                                                               >  1:154426/1‑100 (MQ=255)                                                                                                                                                                                                   tGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATc                                                                                                                                                                                               <  2:154426/100‑1 (MQ=255)                                                                                                                                                                                                    ggggTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAAt                                                                                        <  2:26475/202‑1 (MQ=255)                                                                                                                                                                                                    ggggTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAAt                                                                                        >  1:26475/1‑202 (MQ=255)                                                                                                                                                                                                         gCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTCCTAAAGACAGAGAAACGct                                  >  1:421352/1‑251 (MQ=255)                                                                                                                                                                                                         gCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTCCTAAAGACAGAGAAACGcc                                  >  2:451080/1‑250 (MQ=255)                                                                                                                                                                                                            atcatcCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTCCTAAAGACAGAGAAACGCtctt                               <  1:129187/251‑1 (MQ=255)                                                                                                                                                                                                                  ccGGATCGTAACCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTCCTAACGACAGAGAAACGCTCTTCTAttt                         <  1:449098/251‑1 (MQ=255)                                                                                                                                                                                                                   cGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTCCTAAAGACAGAGAAACGCTCTTCTAtttt                        <  1:331669/251‑1 (MQ=255)                                                                                                                                                                                                                           aCCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGc                 <  1:419998/250‑1 (MQ=255)                                                                                                                                                                                                                                  acATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATa         <  2:485503/250‑1 (MQ=255)                                                                                                                                                                                                                                    aTAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATAcc       <  1:451080/251‑1 (MQ=255)                                                                                                                                                                                                                                          aatgaGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGACATACCGccct  >  2:55465/1‑249 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               ATAACAAACTTGCTTTACGCTGGGGTGACAAATCTTGTTTTATTCAACCACCCAATAGCAGCAACTTAACCACCGGAACGGCTATTTTACCGTGTATTAGCCAAAATTAATTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAATGGTGCGTTTAAATCCCAGTGAGCAAATTCAAATAAGAATGATAACTCAGGAAAAAATTGCTCAACTTCCTAAAGACAGAGAAACGCTCTTCTATTTTAACGTGCGAGAAATACCGCCCA  >  NC_000913/3189756‑3190238

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