Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,133,235	T→G	I221S (ATT→AGT)	flgE →	flagellar hook protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,133,235	0	T	G	100.0%	271.5 / NA	73	I221S (ATT→AGT)	flgE	flagellar hook protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (44/29);  total (44/29)

GCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATAC  >  NC_000913/1133087‑1133382                                                                                                                                                     |                                                                                                                                                    gcCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTg                                                                                                                                                   <  1:486625/151‑1 (MQ=255) gcCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTg                                                                                                                                                   >  2:467477/1‑151 (MQ=255)  ccAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGc                                                                                                                                                  <  1:500213/151‑1 (MQ=255)     gCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGa                                                                                                                                                >  1:13446/1‑150 (MQ=255)       aaTGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGa                                                                                                                                                >  2:555373/1‑148 (MQ=255)       aaTGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGa                                                                                                                                                >  1:359268/1‑148 (MQ=255)        aTGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGAc                                                                                                                                            >  2:673557/1‑151 (MQ=255)         tGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGAc                                                                                                                                            >  1:86795/1‑150 (MQ=255)          gCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACa                                                                                                                                           <  2:245518/150‑1 (MQ=255)           cGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGc                                                                                                                                         <  1:417617/151‑1 (MQ=255)            ggATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCg                                                                                                                                        >  1:668017/1‑151 (MQ=255)                aGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGacaa                                                                                                                                    <  1:301325/151‑1 (MQ=255)                  cTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGacaac                                                                                                                                   >  1:364549/1‑150 (MQ=255)                   tataACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAacac                                                                                                                                 >  2:458753/1‑151 (MQ=255)                   tataACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAacac                                                                                                                                 >  2:398471/1‑151 (MQ=255)                    ataACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGacaaca                                                                                                                                  <  2:379034/149‑1 (MQ=255)                    ataACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACt                                                                                                                                <  1:625167/151‑1 (MQ=255)                    ataACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACt                                                                                                                                >  1:68930/1‑151 (MQ=255)                      aaCAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTgg                                                                                                                              >  1:514551/1‑151 (MQ=255)                          aaaaaGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAAtt                                                                                                                          >  2:465570/1‑151 (MQ=255)                          aaaaaGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAAtt                                                                                                                          >  2:900670/1‑151 (MQ=255)                            aaaGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAAttt                                                                                                                         >  1:30118/1‑150 (MQ=255)                            aaaGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAAttt                                                                                                                         >  1:679052/1‑150 (MQ=255)                            aaaGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTa                                                                                                                        <  1:636090/151‑1 (MQ=255)                               ggTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATg                                                                                                                     >  1:475351/1‑151 (MQ=255)                                gTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGc                                                                                                                    >  1:406579/1‑151 (MQ=255)                                   cGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGCCTACACCCAGGATAGAAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTa                                                                                                                  >  2:742828/1‑150 (MQ=255)                                    ggTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAAt                                                                                                                <  2:443146/151‑1 (MQ=255)                                    ggTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAAt                                                                                                                >  2:899606/1‑151 (MQ=255)                                     gTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATg                                                                                                               >  1:705733/1‑151 (MQ=255)                                      tGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATgg                                                                                                              >  1:675823/1‑151 (MQ=255)                                           gTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGcac                                                                                                           <  1:49019/149‑1 (MQ=255)                                           gTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACAt                                                                                                         >  1:801502/1‑151 (MQ=255)                                             tttCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACAtt                                                                                                        >  1:586286/1‑150 (MQ=255)                                             tttCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACAtt                                                                                                        <  1:580136/150‑1 (MQ=255)                                              ttCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACAtt                                                                                                        >  1:71773/1‑149 (MQ=255)                                                cGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTg                                                                                                    <  1:374101/151‑1 (MQ=255)                                                 gACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTgg                                                                                                   >  2:502546/1‑151 (MQ=255)                                                    aGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATg                                                                                                >  2:110511/1‑151 (MQ=255)                                                         gggTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCg                                                                                           <  1:721377/151‑1 (MQ=255)                                                              aTGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGc                                                                                      <  2:825808/151‑1 (MQ=255)                                                                 cTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCg                                                                                     >  2:713680/1‑149 (MQ=255)                                                                    atgacatgaGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGaataa                                                                                <  1:502559/151‑1 (MQ=255)                                                                       acatgaGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAatat                                                                             <  1:686367/151‑1 (MQ=255)                                                                         atgaGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCg                                                                           <  2:335737/151‑1 (MQ=255)                                                                             gCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCaa                                                                        <  1:664987/150‑1 (MQ=255)                                                                              cGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAAcc                                                                      <  2:484687/151‑1 (MQ=255)                                                                              cGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGCGGATGGTGCGATGGCGAATAATATCGCAAcc                                                                      >  2:371514/1‑151 (MQ=255)                                                                                   aCTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGgc                                                                   <  2:335111/149‑1 (MQ=255)                                                                                     tttGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCa                                                                >  1:629078/1‑150 (MQ=255)                                                                                       tgtgAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAAtt                                                             >  2:478112/1‑151 (MQ=255)                                                                                         tgAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTa                                                            >  1:290483/1‑150 (MQ=255)                                                                                              aCCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTg                                                      <  1:667128/151‑1 (MQ=255)                                                                                                 ggggATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAg                                                   <  2:383964/151‑1 (MQ=255)                                                                                                    gATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAAc                                                >  1:710741/1‑151 (MQ=255)                                                                                                     ataataACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAAcc                                               >  1:864177/1‑151 (MQ=255)                                                                                                       aataaCTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAAccc                                              <  1:371676/150‑1 (MQ=255)                                                                                                        ataaCTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGc                                            >  1:358340/1‑151 (MQ=255)                                                                                                          aaCTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGcc                                           <  1:790311/150‑1 (MQ=255)                                                                                                             tGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGt                                       >  1:801625/1‑151 (MQ=255)                                                                                                                    tCTACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCt                                >  1:517970/1‑151 (MQ=255)                                                                                                                      tACACCCAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTg                               >  1:331085/1‑150 (MQ=255)                                                                                                                          ccaaGGATAGGAGGGATCCAAACGGCGGTGGGAAGACGGCGACAACACGGGAATTTAATGCTAATGGCACATTAGTGGATGGGGCGATGGGGAATAATATGGCAACCGGCGCAATTAGCGGGGCAGAGCCCGCCACGTTTAGTCTGAGc                            <  2:364541/146‑1 (MQ=255)                                                                                                                            cAGGATAGCAGTGATCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTcc                        >  1:451099/1‑151 (MQ=255)                                                                                                                                     aGTGATCAAAACAGCAGTGCGAAGACGGCAACAACACTGAAATTTAATGCTAAGGGCACATTAGTGGTTGGTGCGATGGCGAATAATATCGCAACCGGGGCAATTAACGGTCCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTcc                <  1:555388/150‑1 (MQ=255)                                                                                                                                         aTCCAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATgc            >  2:105360/1‑150 (MQ=255)                                                                                                                                         aTCCAAACAGCAGTGCGAAGACAGCGACAACACTGAAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTCCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATg             <  1:252115/149‑1 (MQ=255)                                                                                                                                           ccAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATAAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATgcagc         <  2:406570/151‑1 (MQ=255)                                                                                                                                            cAAACAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATgcagc         >  2:735493/1‑150 (MQ=255)                                                                                                                                              aaCAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCaa       >  2:753488/1‑150 (MQ=255)                                                                                                                                               aCAGCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCaaa      >  2:121531/1‑150 (MQ=255)                                                                                                                                                cagcagTGCGAAGACAGCGAAAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAAAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAAt    <  1:158038/151‑1 (MQ=255)                                                                                                                                                  gCAGTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATAc  >  1:219016/1‑151 (MQ=255)                                                                                                                                                  gCAGTGAGAAGACAGCGAAAACACGGGAATTAAATGCTAAGGCCACATTAGTGGATGGTGGTATGGCGAATAGTATCGCAGCCGGCGCATTTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATa   <  2:640270/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATAC  >  NC_000913/1133087‑1133382

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