Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,103,122	A→G	A25A (GCT→GCC)	csgD ←	csgBAC operon transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,103,122	0	A	G	100.0%	250.0 / NA	66	A25A (GCT→GCC)	csgD	csgBAC operon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/32);  total (34/32)

TATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGA  >  NC_000913/1102977‑1103269                                                                                                                                                  |                                                                                                                                                    tATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGt                                                                                                                                                >  2:631488/1‑151 (MQ=255)        ccAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGc                                                                                                                                         >  2:512774/1‑151 (MQ=255)        ccAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGc                                                                                                                                         >  1:321291/1‑151 (MQ=255)           ataatGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCaa                                                                                                                                       >  1:783091/1‑150 (MQ=255)           ataatGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCa                                                                                                                                        >  2:340826/1‑149 (MQ=255)             aatGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAaaga                                                                                                                                    <  1:525133/151‑1 (MQ=255)               tGGCTAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAaaga                                                                                                                                    >  2:680631/1‑149 (MQ=255)                 gATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGAt                                                                                                                                <  1:650604/151‑1 (MQ=255)                  aTAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGAtt                                                                                                                               >  2:737150/1‑151 (MQ=255)                  aTAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGAtt                                                                                                                               >  2:686279/1‑151 (MQ=255)                   tAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGAttt                                                                                                                              >  1:750662/1‑151 (MQ=255)                    aaGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGAttt                                                                                                                              <  2:716738/150‑1 (MQ=255)                         ttttATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTg                                                                                                                          <  2:473050/149‑1 (MQ=255)                         ttttATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGa                                                                                                                         <  2:878810/150‑1 (MQ=255)                           ttATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATc                                                                                                                       <  2:229544/150‑1 (MQ=255)                               ccGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATcaaca                                                                                                                   >  2:213798/1‑150 (MQ=255)                               ccGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATcaaca                                                                                                                   >  2:498106/1‑150 (MQ=255)                                 gCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACaata                                                                                                                <  2:609220/151‑1 (MQ=255)                                    tCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACaataa                                                                                                               <  2:770491/149‑1 (MQ=255)                                      catcatATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGt                                                                                                            >  1:875689/1‑150 (MQ=255)                                       atcatATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTAt                                                                                                          <  1:399128/151‑1 (MQ=255)                                        tcatATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATg                                                                                                         >  1:890930/1‑151 (MQ=255)                                         catATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATg                                                                                                         >  2:407745/1‑150 (MQ=255)                                         catATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATg                                                                                                         <  1:675952/150‑1 (MQ=255)                                            atCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGAcc                                                                                                      >  2:537620/1‑150 (MQ=255)                                                    gaaTAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACt                                                                                             <  2:812160/151‑1 (MQ=255)                                                      ataataGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTAt                                                                                           <  2:499467/151‑1 (MQ=255)                                                             agCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTc                                                                                    >  2:618559/1‑151 (MQ=255)                                                               ccTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCAt                                                                                  >  1:605205/1‑151 (MQ=255)                                                                  gAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTaa                                                                               <  2:700887/151‑1 (MQ=255)                                                                       gaTATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACatg                                                                          >  1:299834/1‑151 (MQ=255)                                                                        atatCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACatg                                                                          >  1:508862/1‑150 (MQ=255)                                                                          atCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACatgatg                                                                       <  2:383053/151‑1 (MQ=255)                                                                             gTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGaa                                                                     >  1:320027/1‑150 (MQ=255)                                                                              tCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAAc                                                                   >  1:581052/1‑151 (MQ=255)                                                                                 agagAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAAcccc                                                                >  2:469214/1‑151 (MQ=255)                                                                                 agagAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAAcccc                                                                >  1:62521/1‑151 (MQ=255)                                                                                 agagAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAAcccc                                                                <  2:308585/151‑1 (MQ=255)                                                                                     aaCGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCtt                                                            <  1:752640/151‑1 (MQ=255)                                                                                       cGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTGGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCtttt                                                          <  2:91190/151‑1 (MQ=255)                                                                                       cGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCtttt                                                          <  2:220379/151‑1 (MQ=255)                                                                                        gTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCttttt                                                         >  2:832314/1‑151 (MQ=255)                                                                                          tGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCtttttt                                                        <  2:88662/150‑1 (MQ=255)                                                                                            aaTATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTAt                                                     >  1:666782/1‑151 (MQ=255)                                                                                             atatTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTAtt                                                    <  1:504998/151‑1 (MQ=255)                                                                                                 tatGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGAt                                                 <  1:485369/150‑1 (MQ=255)                                                                                                  atGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATc                                                <  2:149773/150‑1 (MQ=255)                                                                                                  atGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATc                                                <  1:213824/150‑1 (MQ=255)                                                                                                      aaTTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGcac                                            <  1:109765/150‑1 (MQ=255)                                                                                                      aaTTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGcac                                            >  1:234801/1‑150 (MQ=255)                                                                                                       aTTCTCCTGTAACTCCCCGCGATTGTTTACGGTGCTGCAAGAGGGCTGTCGCCTGCCAAGAAGATTTAGTGATCAACATTAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCacac                                          >  2:480888/1‑151 (MQ=255)                                                                                                         tttCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCt                                        >  2:304173/1‑151 (MQ=255)                                                                                                          ttCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCt                                        <  1:545876/150‑1 (MQ=255)                                                                                                             cTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGAc                                     >  2:785113/1‑150 (MQ=255)                                                                                                               gTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAg                                   <  2:214376/150‑1 (MQ=255)                                                                                                                 aaTTGCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTg                                <  1:701530/151‑1 (MQ=255)                                                                                                                     gCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCctc                            >  1:78261/1‑151 (MQ=255)                                                                                                                     gCCAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCctc                            <  2:890859/151‑1 (MQ=255)                                                                                                                       cAGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCctct                           >  1:866946/1‑150 (MQ=255)                                                                                                                        aGCGATTGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTa                          >  1:672772/1‑150 (MQ=255)                                                                                                                              tGTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAg                    <  1:148793/150‑1 (MQ=255)                                                                                                                               gTTTAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGa                   <  1:404972/150‑1 (MQ=255)                                                                                                                                 ttAAGGTGCTGCAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGc                >  1:199030/1‑151 (MQ=255)                                                                                                                                        gctgcAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGa          <  1:831734/150‑1 (MQ=255)                                                                                                                                         ctgcAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGa          >  2:381462/1‑149 (MQ=255)                                                                                                                                           gcAAGAGGGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGt       <  2:872263/150‑1 (MQ=255)                                                                                                                                               gagggCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGa  <  1:631537/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    TATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGA  >  NC_000913/1102977‑1103269

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