Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	715,016	Δ1 bp	12.5%	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,016	0	A	.	12.5%	187.6 / 22.3	64	coding (478/483 nt)	PP_0607	type IV pili biogenesis protein FimT
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base . (4/4);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGCACGCC  >  NC_002947/714869‑715160                                                                                                                                                      |                                                                                                                                                  gACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTagcag                                                                                                                                                  >  1:46048/1‑149 (MQ=255)   cGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGCGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑g                                                                                                                                               >  1:255626/1‑149 (MQ=255)      gcgcGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGt                                                                                                                                             >  3:348908/1‑148 (MQ=255)       cgcgAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCa                                                                                                                                           <  5:181553/149‑1 (MQ=255)         cgAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGCTGCTAGC‑GC‑GGTCAtt                                                                                                                                         <  5:269161/149‑1 (MQ=255)            agagAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCt                                                                                                                                      <  8:230880/149‑1 (MQ=255)              agaACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCa                                                                                                                                    <  3:325700/149‑1 (MQ=255)               gaACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCAGCAAT‑GACCGCTGCTAGC‑GC‑GGTCATTGcagga                                                                                                                                   <  5:230921/149‑5 (MQ=255)               gaACATCACCTTATCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAg                                                                                                                                   <  7:63015/149‑1 (MQ=255)                    tCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGtcct                                                                                                                              >  7:306998/1‑149 (MQ=255)                    tCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGtcct                                                                                                                              <  2:46048/149‑1 (MQ=255)                     cACCTTTTCAATGCCCAAAGAGCTGCTGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCAGCAAT‑GACCGCTGCTAGC‑GC‑GGTCAtt                                                                                                                                         <  5:194068/137‑1 (MQ=255)                     cACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCAGCAAT‑GACCGCTGCTAGC‑GC‑GGTCAtt                                                                                                                                         >  6:194068/1‑137 (MQ=255)                           ttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTaa                                                                                                                        <  6:162495/148‑1 (MQ=255)                           ttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTaa                                                                                                                        <  7:329760/148‑1 (MQ=255)                                  ccgatagagGTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑TGTCATTGCTCAGGTCCTCCGTAAGc                                                                                                                      <  8:610/138‑1 (MQ=255)                                  cccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGTGGCTGCTGTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGATCAGGTCCTCCGTAAGc                                                                                                                      >  7:610/1‑143 (MQ=255)                                   ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGTCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGa                                                                                                               <  8:183061/149‑1 (MQ=255)                                    cAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAc                                                                                                              >  6:360816/1‑149 (MQ=255)                                     aaaGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAc                                                                                                              >  5:180201/1‑148 (MQ=255)                                      aaGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGActc                                                                                                            <  7:303279/149‑1 (MQ=255)                                             gCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCt                                                                                                     <  4:3380/149‑1 (MQ=255)                                              cgtgTTCTCTGGCTGATCTGACGGGCAGTACCGCGGACAATGCAGTGGCTGCTTTGGAGCCGTTGCGAGACACAAGACCTCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAATCATCCTGACTCTGCCGCtt                                                                                                    <  4:348908/146‑1 (MQ=255)                                              cGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGc                                                                                                      <  2:255626/147‑1 (MQ=255)                                              cGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCtt                                                                                                    <  3:199167/149‑1 (MQ=255)                                                     gctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACAAAAGACCGCTCCTGCAAT‑GACCCC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGc                                                                                               >  8:151279/1‑147 (MQ=255)                                                      ctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGCGGCTGCTTTGCAGCCCTTGCGCGACACAACACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAACCAACCTCACTCTGCCGCTTGATCCCac                                                                                            >  8:291972/1‑148 (MQ=255)                                                         gctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGca                                                                                          <  7:291972/148‑1 (MQ=255)                                                                   cGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGatat                                                                               >  6:140626/1‑149 (MQ=255)                                                                          tACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGc                                                                        <  6:180201/149‑1 (MQ=255)                                                                          tACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCCGc                                                                        >  4:176040/1‑149 (MQ=255)                                                                            ccgccgCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGcaccac                                                                                      <  8:220725/133‑1 (MQ=255)                                                                            ccgccgCCAATGCAGGGGCTGCTTGGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGcaccac                                                                                      >  7:220725/1‑133 (MQ=255)                                                                             cgccgcCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGa                                                                                                               <  1:302283/107‑1 (MQ=255)                                                                             cgccgcCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGa                                                                                                               >  2:302283/1‑107 (MQ=255)                                                                                  ccAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAgcg                                                                <  5:360816/149‑1 (MQ=255)                                                                                    aaTGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTGCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCCGCACCTGATATCGGCTGCGAGCCgcggc                                                              >  2:45645/1‑149 (MQ=255)                                                                                    aaTGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCt                                                                          <  7:69807/137‑1 (MQ=255)                                                                                    aaTGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCt                                                                          >  8:69807/1‑137 (MQ=255)                                                                                         aagggCTGCAAAGCAGCCCCTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGCTGGACATTGCTCAGGTCCTCCGTAAGCAACATGACTCTGCCGCTTGATGCCAGCACAACCTGATATCGGATGCGAGCAGCGGCGCGt                                                          >  1:286289/3‑149 (MQ=255)                                                                                           gggCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCg                                                       >  7:71173/1‑149 (MQ=255)                                                                                                    tGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATAGACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCACCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTc                                               >  1:239654/1‑149 (MQ=255)                                                                                                     gCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCgg                                                                            <  5:225660/118‑1 (MQ=255)                                                                                                     gCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCgg                                                                            >  6:225660/1‑118 (MQ=255)                                                                                                     gCAGCCCTAGCGCGACACAAGCCCGCTAATGCAAA‑GACCGC‑GCGAGCAGC‑GGTCATTGCTCAGGGCCTCCGTAACCAACCTGCCTCTCCCCCTAGATCCcaccac                                                                                         >  8:153745/1‑105 (MQ=255)                                                                                                           cTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGc                                                      >  1:282972/1‑134 (MQ=255)                                                                                                           cTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGc                                                      <  2:282972/134‑1 (MQ=255)                                                                                                           cTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTg                                       >  5:175718/1‑149 (MQ=255)                                                                                                            ttGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCGTGATGCCAGCCCCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTgc                                      >  7:168532/1‑149 (MQ=255)                                                                                                              gcgcGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTgcg                                     <  8:306998/148‑1 (MQ=255)                                                                                                              gcgcGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTgc                                      <  2:239654/147‑1 (MQ=255)                                                                                                                     acaAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAAc                               >  2:274110/1‑147 (MQ=255)                                                                                                                            cGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTg                      >  2:164546/1‑149 (MQ=255)                                                                                                                               tCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGccc                   <  8:71173/149‑1 (MQ=255)                                                                                                                               ctgagCAAT‑GACCGCTGCTAGC‑GC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGccc                   >  5:345322/5‑149 (MQ=255)                                                                                                                                     aaT‑GACCGCTGCTAGC‑GC‑GGTCATTGCTCAGGTCCTCCGAAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGt              >  4:142695/1‑148 (MQ=255)                                                                                                                                     aaT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAAc                                                                                                                   <  3:136198/47‑1 (MQ=255)                                                                                                                                     aaT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAAc                                                                                                                   >  4:136198/1‑47 (MQ=255)                                                                                                                                           ccGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATggg        >  3:171810/1‑149 (MQ=255)                                                                                                                                           ccGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATggg        >  8:67305/1‑149 (MQ=255)                                                                                                                                           ccGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATggg        >  8:72187/1‑149 (MQ=255)                                                                                                                                            cgc‑gCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGc       <  2:286289/149‑1 (MQ=255)                                                                                                                                            cgc‑gCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGc       <  4:171810/149‑1 (MQ=255)                                                                                                                                               tgcTAGC‑GC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGc                     >  3:20908/2‑132 (MQ=255)                                                                                                                                               tgcTAGC‑GC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGc                     <  4:20908/131‑1 (MQ=255)                                                                                                                                                   agcagc‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCAGCTTGATGCCAGCACCACCTGATATATGCTGCGAGCAGCGGCGCGTTAGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCAAGAAGTCATGGGCACGcc  >  1:200686/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                  GACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAAT‑GACCGC‑GCTAGCAGC‑GGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGCACGCC  >  NC_002947/714869‑715160

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