Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,430,009	A→G	6.5%	M94V (ATG→GTG)	PP_3924 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,430,009	0	A	G	6.5%	177.1 / 2.5	61	M94V (ATG→GTG)	PP_3924	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/25);  new base G (2/2);  total (34/27)
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.44e-01

CGGCACCTGCGACCTGCAGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACA  >  NC_002947/4429861‑4430156                                                                                                                                                     |                                                                                                                                                    cGGCACCTGCGACCTGCAGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACatga                                                                                                                                                     <  6:161322/149‑1 (MQ=255) cGGCACCTGCGACCTGCAGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACatga                                                                                                                                                     <  6:181897/149‑1 (MQ=255)   gCACCTGCGACCTGCAGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACatgatg                                                                                                                                                   <  8:143911/149‑1 (MQ=255)        tGCGACCTGCAGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTc                                                                                                                                              >  5:124392/1‑149 (MQ=255)          cGACCTGCAGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGGAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGa                                                                                                                                            <  3:266554/149‑1 (MQ=255)          cGACCTGCAGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGa                                                                                                                                            <  6:124392/149‑1 (MQ=255)          cGACCTGCAGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGGGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTc                                                                                                                                              <  8:192128/147‑1 (MQ=255)            aCCTGCAGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATg                                                                                                                                          <  7:187771/149‑1 (MQ=255)              cTGCAGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGAt                                                                                                                                           <  7:192209/146‑1 (MQ=255)              cTGCAGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGAt                                                                                                                                           >  8:192209/1‑146 (MQ=255)              cTGCAGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGa                                                                                                                                         >  6:186116/1‑148 (MQ=255)               tGCAGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGaaa                                                                                                                                       >  3:126456/1‑149 (MQ=255)                 cAGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGaaa                                                                                                                                       <  4:126456/147‑1 (MQ=255)                  aGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGAGGCGGGTTGAAGGCGACATGATCCTGTCGATGAAACt                                                                                                                                     >  6:173346/1‑148 (MQ=255)                     aaGGCGAAAACCCGGATGCCAACTGCACGCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGc                                                                                                                                 <  2:217084/149‑1 (MQ=255)                              aCCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTc                                                                                                                        <  5:186116/149‑1 (MQ=255)                                         aaCTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAgg                                                                                                             >  7:92926/1‑149 (MQ=255)                                                 ccTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGaa                                                                                                     >  3:44244/1‑149 (MQ=255)                                                     gTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATAcc                                                                                                 >  8:194773/1‑149 (MQ=255)                                                      tGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGa                                                                                                                                            <  6:95244/105‑1 (MQ=255)                                                      tGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGa                                                                                                                                            >  5:95244/1‑105 (MQ=255)                                                        aTGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCACGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGga                                                                                              <  1:268375/149‑1 (MQ=255)                                                            aCGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTgg                                                                                                        <  5:35834/135‑1 (MQ=255)                                                            aCGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTgg                                                                                                        >  6:35834/1‑135 (MQ=255)                                                                cGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTg                                                                                                                                                >  8:156936/1‑91 (MQ=255)                                                                cGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTg                                                                                                                                                <  7:156936/91‑1 (MQ=255)                                                                   taCCCTTCAGGGTTTCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGaaaaa                                                                                   >  4:124279/2‑149 (MQ=255)                                                                      ccTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTGCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAAcc                                                                                 >  6:62401/1‑148 (MQ=255)                                                                       cTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTcc                                                                                                                  <  8:294990/114‑1 (MQ=255)                                                                       cTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTcc                                                                                                                  >  7:294990/1‑114 (MQ=255)                                                                        tGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGGTGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATAc                                                                                                  >  6:188912/1‑129 (MQ=255)                                                                        tGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGGTGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATAc                                                                                                  <  1:107890/129‑1 (MQ=255)                                                                        tGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGGTGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATAc                                                                                                  <  5:188912/129‑1 (MQ=255)                                                                        tGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGGTGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATAc                                                                                                  >  2:107890/1‑129 (MQ=255)                                                                              gTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGa                                                                        <  8:92926/149‑1 (MQ=255)                                                                                    tAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCtt                                                                  >  6:142089/1‑149 (MQ=255)                                                                                       gcggcgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGTTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCGTCAGGCTTCgg                                                               <  5:173346/149‑1 (MQ=255)                                                                                       gcggcgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCgg                                                               >  2:22649/1‑149 (MQ=255)                                                                                            gAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTc                                                                                                                        >  2:50615/1‑87 (MQ=255)                                                                                            gAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTc                                                                                                                        <  1:50615/87‑1 (MQ=255)                                                                                                  ccacccATGAAGGCCTGCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATAcc                                                                                                 >  5:197628/6‑104 (MQ=255)                                                                                                  ccacccATGAAGGCCTGCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATAcc                                                                                                 <  6:197628/99‑1 (MQ=255)                                                                                                      gCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGc                                                <  7:194773/149‑1 (MQ=255)                                                                                                         tGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCAtt                                             >  5:166975/1‑149 (MQ=255)                                                                                                          gCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGaa                                                                              >  8:155895/1‑115 (MQ=255)                                                                                                          gCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGaa                                                                              <  7:155895/115‑1 (MQ=255)                                                                                                           cAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGc                                                                                                                                 <  2:253141/63‑1 (MQ=255)                                                                                                           cAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGc                                                                                                                                 >  1:253141/1‑63 (MQ=255)                                                                                                            aGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGc                                                                    <  2:44393/123‑1 (MQ=255)                                                                                                            aGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGc                                                                    >  1:44393/1‑123 (MQ=255)                                                                                                                 ttCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAg                                     >  3:78480/1‑149 (MQ=255)                                                                                                                        gTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATcca                              >  8:228478/1‑149 (MQ=255)                                                                                                                        gTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGACGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGtt                                                             <  5:68851/118‑1 (MQ=255)                                                                                                                        gTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGACGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGtt                                                             >  6:68851/1‑118 (MQ=255)                                                                                                                               gCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCaa                       >  2:15140/1‑149 (MQ=255)                                                                                                                               gCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCaa                       <  3:124279/149‑1 (MQ=255)                                                                                                                                   cGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCTTTCGGTGCAGTGATCCACCAGCAacac                    <  1:15140/148‑1 (MQ=255)                                                                                                                                      gTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCt                 >  7:243913/1‑148 (MQ=255)                                                                                                                                           aGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTa           >  1:33194/1‑149 (MQ=255)                                                                                                                                             gCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTg                                        >  6:46169/1‑118 (MQ=255)                                                                                                                                             gCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTg                                        <  5:46169/118‑1 (MQ=255)                                                                                                                                               gACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTg       >  2:89259/1‑149 (MQ=255)                                                                                                                                               gACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTg       >  1:242111/1‑149 (MQ=255)                                                                                                                                                aCATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTgg      >  5:67244/1‑149 (MQ=255)                                                                                                                                                    gatgCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGt                                         >  2:178758/1‑110 (MQ=255)                                                                                                                                                    gatgCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGt                                         <  1:178758/110‑1 (MQ=255)                                                                                                                                                    gatgCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGcaca  >  4:48243/1‑149 (MQ=255)                                                                                                                                                     atgCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGcaca  <  2:33194/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CGGCACCTGCGACCTGCAGGAAGGCGAAAACCCGGATGCCAACTGCACCCTGGTGATGGACGCCGAAACCCTGAAGGGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACA  >  NC_002947/4429861‑4430156

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