Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	715,019	Δ1 bp	14.1%	coding (478/483 nt)	PP_0607 ←	type IV pili biogenesis protein FimT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	715,019	0	A	.	14.1%	180.9 / 26.9	64	coding (478/483 nt)	PP_0607	type IV pili biogenesis protein FimT
Reads supporting (aligned to +/- strand):  ref base A (28/27);  new base . (6/3);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

GGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGCACGCCCAGC  >  NC_002947/714871‑715167                                                                                                                                                      |                                                                                                                                                     ggACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTagcag                                                                                                                                                     >  1:316507/1‑150 (MQ=255)           gAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGc                                                                                                                                           >  1:68840/1‑150 (MQ=255)                      cACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGtcctcc                                                                                                                                >  2:1199295/1‑150 (MQ=255)                       aCCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCg                                                                                                                               >  2:1298094/1‑150 (MQ=255)                         cTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTa                                                                                                                             >  1:842018/1‑150 (MQ=255)                          ttttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGGCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTaa                                                                                                                            <  1:1670799/150‑1 (MQ=255)                           tttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAg                                                                                                                           <  2:68840/150‑1 (MQ=255)                                    ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCTACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAc                                                                                                                  <  1:301282/150‑1 (MQ=255)                                    ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGC‑GCGGTCATTGcag                                                                                                                                         <  2:951214/127‑3 (MQ=39)                                    ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGC‑GCGGTCATTGcag                                                                                                                                         >  1:951196/1‑125 (MQ=39)                                    ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGC‑GCGGTCATTGcag                                                                                                                                         >  1:951214/1‑125 (MQ=39)                                    ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGC‑GCGGTCATTGcag                                                                                                                                         <  2:951196/127‑3 (MQ=39)                                    ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAc                                                                                                                  <  2:471918/150‑1 (MQ=255)                                     cAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAct                                                                                                                 >  2:1031303/1‑150 (MQ=255)                                          agCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGcc                                                                                                            <  1:1298094/150‑1 (MQ=255)                                                 ggTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGAt                                                                                                     <  2:316507/150‑1 (MQ=255)                                                     cgctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCa                                                                                                 <  2:842018/150‑1 (MQ=255)                                                      gctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAg                                                                                                                                       >  2:1298189/1‑111 (MQ=255)                                                      gctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAg                                                                                                                                       <  1:1298189/111‑1 (MQ=255)                                                              aTCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCTTTGCGCGAAACAAGACCGCTCCTGCAATGACAGC‑GCTAGCAGCGGTCATTGCTCATGTCCTCCGTAAGCAACCTGACTCTCCCGCTTGATGCCAGCACCACCt                                                                                        <  2:220822/150‑1 (MQ=255)                                                              aTCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTc                                                                                                                                         <  2:647055/101‑1 (MQ=255)                                                              aTCTGACGGGCAGTACCGCCGCCAATGCAGGGGATGCTGTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTc                                                                                                                                         >  1:647055/1‑101 (MQ=255)                                                               tCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTg                                                                                       >  2:1516790/1‑150 (MQ=255)                                                                     gggCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCg                                                                                 >  1:1084229/1‑150 (MQ=255)                                                                      ggCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCgg                                                                                <  1:488900/150‑1 (MQ=255)                                                                         aGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTg                                                                             <  1:985926/150‑1 (MQ=255)                                                                         aGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTg                                                                             >  1:705840/1‑150 (MQ=255)                                                                          gTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGc                                                                            >  2:978382/1‑150 (MQ=255)                                                                          gTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGc                                                                            >  1:240938/1‑150 (MQ=255)                                                                          gTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATACCAGCACCACCTGATATCGGCTGc                                                                            >  1:1734360/1‑150 (MQ=255)                                                                           tACCGCCGCAATTGAAGGGGCTGCTTGGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCg                                                                           >  2:1392161/1‑150 (MQ=255)                                                                            cccgacgCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCTACACAAGACCGCTCCTGCAATGACCGC‑GCTAGGAGCAGTCATTTCTCAGGTCCTCTGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCGCCTGATATCGGCTGCGa                                                                          <  2:1740389/149‑1 (MQ=255)                                                                             ccgccgCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCg                                                                                                           <  1:1100887/116‑1 (MQ=255)                                                                             ccgccgCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCg                                                                                                           >  2:1100887/1‑116 (MQ=255)                                                                                 cgccAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGagcagc                                                                     <  2:625355/150‑1 (MQ=255)                                                                                      aTGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGgcgc                                                                <  2:705840/150‑1 (MQ=255)                                                                                        gCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGt                                                              <  1:1516790/150‑1 (MQ=255)                                                                                         cAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTCCACCGTTCCTCCGTAAGCGCCCGGAGTCGGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGtt                                                             <  1:236436/150‑1 (MQ=255)                                                                                          aGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGc                                                                            >  1:1434140/1‑134 (MQ=255)                                                                                          aGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGc                                                                            <  2:1434140/134‑1 (MQ=255)                                                                                              gctgctTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCGAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGcac                                                        >  2:1726961/1‑150 (MQ=255)                                                                                                tgctTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTGAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATGGGCTGCGAGCAGCGGCGCGTTAGCACAt                                                      >  1:1788923/1‑150 (MQ=255)                                                                                                tgctTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACAt                                                      >  1:1402980/1‑150 (MQ=255)                                                                                                     tGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCa                                                 <  1:1199295/150‑1 (MQ=255)                                                                                                           ccTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTg                                           >  2:646244/1‑150 (MQ=255)                                                                                                              tGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTgcgc                                        <  2:1402980/150‑1 (MQ=255)                                                                                                                      acaAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAg                                >  2:617095/1‑150 (MQ=255)                                                                                                                           aCCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCGGGTCCTCCGTAAGCAACCTGACTCTGCCGCtt                                                                                                        >  1:1851493/1‑73 (MQ=255)                                                                                                                           aCCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCtt                                                                                                        <  2:1851493/73‑1 (MQ=255)                                                                                                                           aCCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGtt                           <  1:1536393/150‑1 (MQ=255)                                                                                                                           aCCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGtt                           >  1:189375/1‑150 (MQ=255)                                                                                                                               cTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGccc                       >  1:376821/1‑150 (MQ=255)                                                                                                                                tCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCg                      <  1:1031303/150‑1 (MQ=255)                                                                                                                                 tgagCAATGACCGCTGCTAGC‑GCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGa                     >  1:1565167/4‑150 (MQ=255)                                                                                                                                   tGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGa                                                                          >  2:1897516/1‑95 (MQ=255)                                                                                                                                   tGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGa                                                                          <  1:1897516/95‑1 (MQ=255)                                                                                                                                   tGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACg                   >  1:498037/1‑150 (MQ=255)                                                                                                                                     cAATGACCGCTGCTAGC‑GCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACGACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTc                 >  2:771809/1‑150 (MQ=255)                                                                                                                                      aaTGACCGCTGCTAGC‑GCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCg                                                                                 <  1:1448624/85‑1 (MQ=38)                                                                                                                                      aaTGACCGCTGCTAGC‑GCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCg                                                                                 >  2:1448624/1‑85 (MQ=38)                                                                                                                                      aaTGACCGC‑GCTAGCAGCGGTCATTGCTCAGGtcctcc                                                                                                                                >  1:1501705/1‑38 (MQ=255)                                                                                                                                      aaTGACCGC‑GCTAGCAGCGGTCATTGCTCAGGtcctcc                                                                                                                                <  2:1501705/38‑1 (MQ=255)                                                                                                                                       aTGACCGCTGCTAGC‑GCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACTTTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCAt               >  2:118813/1‑150 (MQ=255)                                                                                                                                         gACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCAAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATgg             <  1:646244/150‑1 (MQ=255)                                                                                                                                                  tAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCAGTTGCCCGACGTCATGGgcccgccc     >  2:1850541/1‑150 (MQ=255)                                                                                                                                                   agcagcGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGCACGCCCa    <  1:1215897/150‑1 (MQ=255)                                                                                                                                                     cagcGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGCACGCCCAgc  <  1:617095/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGCACGCCCAGC  >  NC_002947/714871‑715167

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑