Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,436,184	A→C	6.3%	N40T (AAC→ACC)	gpI →	tail formation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,436,184	0	A	C	6.3%	199.7 / 4.7	64	N40T (AAC→ACC)	gpI	tail formation protein
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base C (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

GTTACGCTGACGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGC  >  NC_002947/3436040‑3436331                                                                                                                                                  |                                                                                                                                                     gTTACGCTGACGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTc                                                                                                                                                  >  1:331140/1‑149 (MQ=255) gTTACGCTGACGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTc                                                                                                                                                  <  8:149899/149‑1 (MQ=255)   tACGCTGACGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑t                                                                                                                                               >  5:163954/1‑149 (MQ=255)     cGCTGACGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGc                                                                                                                                             >  6:196016/1‑149 (MQ=255)          aCGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGCGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCct                                                                                                                                          <  1:14810/147‑1 (MQ=255)          aCGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCctgg                                                                                                                                        <  6:163954/149‑1 (MQ=255)               gAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCctggctgg                                                                                                                                    <  4:11314/148‑1 (MQ=255)               gAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCtggctggc                                                                                                                                   <  3:298356/149‑1 (MQ=255)                         aCATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCAGCCC‑TCGACCTTCCATGCCCAGGCCAGCCAGGGCagga                                                                                                                         <  5:319179/149‑4 (MQ=37)                           aTGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATgg                                                                                                                         >  7:132172/1‑147 (MQ=255)                                   gTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCctggctgg                                                                                                                                    <  8:319054/128‑1 (MQ=255)                                   gTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCctggctgg                                                                                                                                    >  7:319054/1‑128 (MQ=255)                                      aCCGCACAACACCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCt                                                                                                            >  6:355225/1‑149 (MQ=255)                                        cGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCCGACTCGTCACGCTGCCCGCCGAACTTc                                                                                                                                                  <  1:208758/111‑1 (MQ=255)                                        cGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTc                                                                                                                                                  >  2:208758/1‑110 (MQ=255)                                        cGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTgg                                                                                                          <  2:256691/149‑1 (MQ=255)                                        cGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTCGGCATGGAAGGTCGAGGGCTgg                                                                                                          >  7:152688/1‑149 (MQ=255)                                             aaCGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggc                                                                                                     >  8:92122/1‑149 (MQ=255)                                             aaCGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggc                                                                                                     <  7:201739/149‑1 (MQ=255)                                                gCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggcggc                                                                                                  <  2:331140/149‑1 (MQ=255)                                                  cACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCt                                                                                                 >  6:92390/1‑148 (MQ=255)                                                    cccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTa                                                                                               >  6:197777/1‑148 (MQ=255)                                                     ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                             <  5:92390/149‑1 (MQ=255)                                                     ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                             <  8:132172/149‑1 (MQ=255)                                                     ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                             >  1:338598/1‑149 (MQ=255)                                                            gAGCGCGCCCTGGAGGTGGCCGACGACCTTGGGCTTGAGCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGGGCTGGCTGGGGTGGGCATGGAAGGTCGTGGGCTGGGAGGCGGCt                                                                                                 <  8:184519/138‑1 (MQ=255)                                                            gAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATACGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAAGGAAGGTCGAGGGCTGGGAGGCGGCt                                                                                                 >  7:184519/1‑138 (MQ=255)                                                              gcgcgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGAcc                                                                                    >  2:133000/1‑149 (MQ=255)                                                                    ccTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGcgcc                                                                              <  5:355225/149‑1 (MQ=255)                                                                           gcggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCt                                                                       >  6:70162/1‑149 (MQ=255)                                                                              gccgccGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTg                                                                                                                                >  3:289747/1‑89 (MQ=255)                                                                              gccgccGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTg                                                                                                                                <  4:289747/89‑1 (MQ=255)                                                                                cgccgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATcc                                                                  >  6:121314/1‑149 (MQ=255)                                                                                   cgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCcgc                                                                >  6:232090/1‑148 (MQ=255)                                                                                    gACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGa                                                              <  1:136321/149‑1 (MQ=255)                                                                                     aCCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAg                                                             <  7:92122/149‑1 (MQ=255)                                                                                            ccTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTc                                                      >  5:246400/1‑149 (MQ=255)                                                                                                gACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCg                                                    <  5:232090/147‑1 (MQ=255)                                                                                                 aCCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGtt                                                 <  1:133000/149‑1 (MQ=255)                                                                                                 aCCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGt                                                  <  5:70162/148‑1 (MQ=255)                                                                                                    cGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                              <  6:246400/149‑1 (MQ=255)                                                                                                    cGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                              <  5:197777/149‑1 (MQ=255)                                                                                                    cGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                              <  5:196016/149‑1 (MQ=255)                                                                                                    cGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                               >  6:337837/1‑148 (MQ=255)                                                                                                         aTTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGAc                                         >  8:49438/1‑149 (MQ=255)                                                                                                         aTTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGAc                                         >  7:105605/1‑149 (MQ=255)                                                                                                             tCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCgcgcgc                                                                        <  1:82390/114‑1 (MQ=255)                                                                                                             tCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCgcgcgc                                                                        >  2:82390/1‑114 (MQ=255)                                                                                                                    cGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACgg                              <  5:121314/149‑1 (MQ=255)                                                                                                                        g‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCg                           <  7:49438/148‑1 (MQ=255)                                                                                                                            gACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcgc                       <  8:105605/149‑1 (MQ=255)                                                                                                                                  tCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGTGCGATTcg                 >  7:163851/1‑149 (MQ=255)                                                                                                                                  tCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcg                 >  1:35385/1‑149 (MQ=255)                                                                                                                                    aCGCTGCCCGCCGAACTTCC‑TGCCATGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCAtt                                                       >  5:97226/1‑109 (MQ=255)                                                                                                                                    aCGCTGCCCGCCGAACTTCC‑TGCCATGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCAtt                                                       <  6:97226/109‑1 (MQ=255)                                                                                                                                      gCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAg                                                             <  5:356191/101‑1 (MQ=255)                                                                                                                                      gCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAg                                                             >  6:356191/1‑101 (MQ=255)                                                                                                                                      gCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg             >  2:315238/1‑149 (MQ=255)                                                                                                                                        tGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg             >  3:131811/1‑147 (MQ=255)                                                                                                                                        tGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg             <  4:131811/147‑1 (MQ=255)                                                                                                                                            cctcgACCTTCCATGCCCAGGCCAGCCAGGGCA‑GGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                     >  5:345476/4‑135 (MQ=25)                                                                                                                                            cctcgACCTTCCATGCCCAGGCCAGCCAGGGCA‑GGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                     <  6:345476/132‑1 (MQ=25)                                                                                                                                             ctcgACCTTCCATGCCCAGGCCAGCCAGGGCA‑GGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACAAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTc       >  7:252213/3‑148 (MQ=37)                                                                                                                                                gAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAgg   <  2:35385/149‑1 (MQ=255)                                                                                                                                                 aaCTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGc  <  1:315238/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GTTACGCTGACGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGG‑CCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGC  >  NC_002947/3436040‑3436331

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