Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,866,767	C→T	W262* (TGG→TAG)	rpoS ←	RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,866,767	0	C	T	100.0%	241.2 / NA	64	W262* (TGG→TAG)	rpoS	RNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/37);  total (27/37)

CAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGAC  >  NC_000913/2866622‑2866911                                                                                                                                                  |                                                                                                                                                 cAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTg                                                                                                                                             <  2:246476/151‑1 (MQ=255)       ttCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATg                                                                                                                                       >  1:676746/1‑151 (MQ=255)                  tCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGcatcat                                                                                                                            >  1:856585/1‑151 (MQ=255)                   cTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCat                                                                                                                            >  2:1112161/1‑150 (MQ=255)                     ggCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCat                                                                                                                            <  2:154878/148‑1 (MQ=255)                      gCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCg                                                                                                                        <  1:28035/151‑1 (MQ=255)                         aaCACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCa                                                                                                                     <  1:311420/151‑1 (MQ=255)                           cacGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATc                                                                                                                   >  2:1034763/1‑151 (MQ=255)                            acGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATc                                                                                                                   <  1:626608/150‑1 (MQ=255)                             cGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCtt                                                                                                                 >  1:1051130/1‑151 (MQ=255)                              gTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTg                                                                                                                >  2:225750/1‑151 (MQ=255)                              gTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTg                                                                                                                <  2:838267/151‑1 (MQ=255)                                tCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGc                                                                                                               <  1:329016/150‑1 (MQ=255)                                 cACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGt                                                                                                             <  1:55706/151‑1 (MQ=255)                                  aCGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTg                                                                                                            <  2:567968/151‑1 (MQ=255)                                  aCGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTg                                                                                                            <  1:1049001/151‑1 (MQ=255)                                    gtgTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCTATCAAACCGAATCGCCGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTCCATATCGTCATCTTGCGTg                                                                                                            <  1:957394/147‑1 (MQ=255)                                          ggCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTc                                                                                                    >  2:533405/1‑151 (MQ=255)                                           gCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTcc                                                                                                   >  2:96285/1‑151 (MQ=255)                                            ccAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCg                                                                                                  <  2:908293/151‑1 (MQ=255)                                             cAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCgg                                                                                                 >  2:988731/1‑151 (MQ=255)                                              aaTTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCgg                                                                                                 >  1:845128/1‑150 (MQ=255)                                                       aCCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTatct                                                                                       >  1:592522/1‑151 (MQ=255)                                                          tACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCtttt                                                                                    >  2:96276/1‑151 (MQ=255)                                                           aCATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCttttt                                                                                   <  2:191171/151‑1 (MQ=255)                                                           aCATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCtttt                                                                                    <  1:521140/150‑1 (MQ=255)                                                            cATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTc                                                                                  >  1:803941/1‑151 (MQ=255)                                                              tCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCAt                                                                                >  2:490416/1‑151 (MQ=255)                                                                 tCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCa                                                                                 <  2:68854/147‑1 (MQ=255)                                                                 tCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCg                                                                              <  1:255897/150‑1 (MQ=255)                                                                 tCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCg                                                                              <  2:1089877/150‑1 (MQ=255)                                                                   cAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGcc                                                                           <  1:907290/151‑1 (MQ=255)                                                                       gtTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAgg                                                                        <  2:226222/150‑1 (MQ=255)                                                                       gtTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGa                                                                       <  2:626640/151‑1 (MQ=255)                                                                            cGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTcc                                                                  >  1:188452/1‑151 (MQ=255)                                                                                 cGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCa                                                              <  2:866004/150‑1 (MQ=255)                                                                                    aCCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGc                                                          >  1:663517/1‑151 (MQ=255)                                                                                             aaCCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCgg                                                  >  2:1002237/1‑150 (MQ=255)                                                                                             aaCCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCgg                                                  >  1:1065096/1‑150 (MQ=255)                                                                                             aaCCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGa                                                 <  1:1022946/151‑1 (MQ=255)                                                                                                cGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATc                                              >  2:435223/1‑151 (MQ=255)                                                                                                cGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATc                                              <  2:621025/151‑1 (MQ=255)                                                                                                  aaTCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATcac                                            <  1:904296/151‑1 (MQ=255)                                                                                                   aTCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATcacc                                           >  2:679995/1‑151 (MQ=255)                                                                                                      gACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCaccacc                                        <  2:811913/151‑1 (MQ=255)                                                                                                      gACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCaccacc                                        >  2:972047/1‑151 (MQ=255)                                                                                                        cGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCa                                      <  1:1090217/151‑1 (MQ=255)                                                                                                          tGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGc                                    >  1:260260/1‑151 (MQ=255)                                                                                                           gCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGc                                    <  2:938759/150‑1 (MQ=255)                                                                                                                 aCTTCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGgtgt                             <  2:1008588/151‑1 (MQ=255)                                                                                                                    tCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCt                           <  1:1037853/150‑1 (MQ=255)                                                                                                                    tCACGCTGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCt                           <  1:553769/150‑1 (MQ=255)                                                                                                                     cACGCTGTTTGGCGTTCAGTTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTa                          <  1:236531/150‑1 (MQ=255)                                                                                                                          tGTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAgg                    >  1:236014/1‑151 (MQ=255)                                                                                                                           gTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAgg                    <  2:1052867/150‑1 (MQ=255)                                                                                                                           gTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAgg                    <  2:1047911/150‑1 (MQ=255)                                                                                                                           gTTTGGCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAgg                    <  1:573091/150‑1 (MQ=255)                                                                                                                                gCGTTCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATgc              >  2:686237/1‑151 (MQ=255)                                                                                                                                   ttCAGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCt           >  1:464141/1‑151 (MQ=255)                                                                                                                                      aGCTCGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGt        >  1:504685/1‑151 (MQ=255)                                                                                                                                          cGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTaa     <  2:750251/150‑1 (MQ=255)                                                                                                                                          cGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTa      >  2:88809/1‑149 (MQ=255)                                                                                                                                          cGAACAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAg    <  2:336059/151‑1 (MQ=255)                                                                                                                                            aaCAGCTATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGAc  <  2:1091611/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGAC  >  NC_000913/2866622‑2866911

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: ‑