Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,054,949	A→G	5.1%	Q60Q (CAA→CAG)	PP_3574 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,054,949	0	A	G	5.1%	233.4 / 2.0	78	Q60Q (CAA→CAG)	PP_3574	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/39);  new base G (2/2);  total (37/41)
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.19e-01

CCTGGACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTTGGTGGTGGCC  >  NC_002947/4054805‑4055096                                                                                                                                                 |                                                                                                                                                    ccTGGACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGc                                                                                                                                                 >  6:384696/1‑149 (MQ=255) ccTGGACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGc                                                                                                                                                 >  3:82695/1‑149 (MQ=255)   tGGACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCtt                                                                                                                                               >  5:303256/1‑149 (MQ=255)        aaCGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAGTTGCTTgacgt                                                                                                                                          <  3:76/149‑4 (MQ=255)        aaCGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTgaag                                                                                                                                           >  4:125032/1‑148 (MQ=255)         aCGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCGACTGCTTgaag                                                                                                                                           >  8:250355/1‑147 (MQ=255)            aCTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGcc                                                                                                                                      <  3:146923/149‑1 (MQ=255)              ttGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCgg                                                                                                                                    >  6:391096/1‑149 (MQ=255)                gtgtCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCGACTGCTTGAAGAAGCCgg                                                                                                                                    <  7:250355/147‑1 (MQ=255)                 tgtCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCt                                                                                                                                  <  7:132381/148‑1 (MQ=255)                 tgtCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCt                                                                                                                                  <  2:90985/148‑1 (MQ=255)                 tgtCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTc                                                                                                                                 >  3:49797/1‑149 (MQ=255)                   tCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCTGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCgg                                                                                                                               <  3:365503/149‑1 (MQ=255)                   tCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCg                                                                                                                                >  3:356803/1‑148 (MQ=255)                   tCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCg                                                                                                                                <  4:49797/148‑1 (MQ=255)                       gCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCt                                                                                                                           >  4:306517/1‑149 (MQ=255)                          gTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTc                                                                                                                        <  2:92815/149‑1 (MQ=255)                                gcgggcAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACAt                                                                                                                  >  7:332120/1‑149 (MQ=255)                                    gcAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGa                                                                                                                              <  4:180600/133‑1 (MQ=255)                                    gcAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGa                                                                                                                              >  3:180600/1‑133 (MQ=255)                                      aaTACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCaa                                                                                                            <  1:166379/149‑1 (MQ=255)                                       aTACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAg                                                                                                           <  5:384696/149‑1 (MQ=255)                                        tACCATCTACGTTCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGa                                                                                                          <  7:197144/149‑1 (MQ=255)                                         aCCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGac                                                                                                         >  3:217783/1‑149 (MQ=255)                                           cATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGAc                                                                                                                             >  7:104998/1‑127 (MQ=255)                                           cATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGAc                                                                                                                             <  8:104998/127‑1 (MQ=255)                                            aTCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGaccac                                                                                                      >  3:127079/1‑149 (MQ=255)                                            aTCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGaccac                                                                                                      >  3:120784/1‑149 (MQ=255)                                              cTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGaccacca                                                                                                    >  5:331159/1‑149 (MQ=255)                                                 cGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGaccacc                                                                                                     <  3:306517/145‑1 (MQ=255)                                                 cGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGaccacc                                                                                                     <  5:50305/145‑1 (MQ=255)                                                   tGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTa                                                                                                >  2:154687/1‑148 (MQ=255)                                                   tGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTAc                                                                                               <  4:217783/149‑1 (MQ=255)                                                     cgcgGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACa                                                                                              >  4:69402/1‑148 (MQ=255)                                                         gCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCAcc                                                                                         <  4:356803/149‑1 (MQ=255)                                                            aGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGAc                                                                                      <  2:24902/149‑1 (MQ=255)                                                            aGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGAc                                                                                      <  8:332120/149‑1 (MQ=255)                                                                  gTACCGACTTCGAAGGCAAGCTGTTGCGCCTTGGTGACCAGCCCGCACAAACCCAACAGGCGATGAAGAACGTCAAGCCACTGCTAGAATAAGCCGGCTCGGACCTGTCGCACATCGACAAGACCACCACATACATCACCGTCCcgccc                                                                                >  5:411528/1‑147 (MQ=255)                                                                      cGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGCACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTcc                                                                            >  5:190301/1‑149 (MQ=255)                                                                      cGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCg                                                                                                                                     <  1:48274/92‑1 (MQ=255)                                                                      cGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCg                                                                                                                                     >  2:48274/1‑92 (MQ=255)                                                                       gACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTc                                                                                                                                 <  6:361439/95‑1 (MQ=255)                                                                       gACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTc                                                                                                                                 >  5:361439/1‑95 (MQ=255)                                                                          ttCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCcg                                                                           >  8:318405/1‑146 (MQ=255)                                                                          ttCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCcg                                                                           <  7:318405/146‑1 (MQ=255)                                                                            cGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGa                                                                        <  4:120784/147‑1 (MQ=255)                                                                            cGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGa                                                                        <  3:125032/147‑1 (MQ=255)                                                                             gAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGc                                                                      <  7:163601/148‑1 (MQ=255)                                                                              aaGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCg                                                                    <  1:23957/149‑1 (MQ=255)                                                                              aaGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCg                                                                    <  1:245925/149‑1 (MQ=255)                                                                                ggCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGgt                                                                  >  8:164122/1‑149 (MQ=255)                                                                                       gggtggGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGAGGAAGACCGTCACGCAACTGCTTGAAGACGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCCCCTACATCACCGCCCCGCGCGTACGCGAGCCGGTGTCCAag                                                           >  8:123922/2‑149 (MQ=255)                                                                                         gtggGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAagg                                                          <  4:82695/148‑1 (MQ=255)                                                                                              ccTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACAGCTTGAAGAAGCCGGCTCGGACCTTTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTgg                                                    <  2:210541/149‑1 (MQ=255)                                                                                               cTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCcgc                                                                                  >  5:158703/1‑118 (MQ=255)                                                                                               cTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCcgc                                                                                  <  6:158703/118‑1 (MQ=255)                                                                                               cTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTggg                                                   <  5:391096/149‑1 (MQ=255)                                                                                                  ggTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGc                                                  >  8:9890/1‑147 (MQ=255)                                                                                                  ggTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCGCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGc                                                  <  6:303256/147‑1 (MQ=255)                                                                                                   gTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTTGGCAAg                                               <  4:127079/149‑1 (MQ=255)                                                                                                        ccGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCt                                          >  7:4288/1‑149 (MQ=255)                                                                                                           cgcgcACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAagg                                                          >  1:218411/1‑130 (MQ=255)                                                                                                           cgcgcACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAagg                                                          <  2:218411/130‑1 (MQ=255)                                                                                                             cgcACAAACCGAACAGGCGATGAAGAACGTCAATAAACTGCTTGAAGAAGCCGTCTCGGACCTGTCTCACATCGTCAAGACCACCACCTACATCAACGACCCGCGCTACAGCGACCCGCTGTACAAGGAGGTGGGCAACAGGCACAAg                                      >  5:399293/1‑148 (MQ=255)                                                                                                                  aaaCCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTc                                         <  4:204937/140‑1 (MQ=255)                                                                                                                  aaaCCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTc                                         >  3:204937/1‑140 (MQ=255)                                                                                                                     ccGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTAAATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTc                             >  4:377032/1‑149 (MQ=255)                                                                                                                        aaCAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTccc                           >  8:208636/1‑148 (MQ=255)                                                                                                                            ggCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCGCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGAttt                      <  3:69402/149‑1 (MQ=255)                                                                                                                                    agaaCGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGa                         >  6:3311/1‑138 (MQ=255)                                                                                                                                    agaaCGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGa                         <  5:3311/138‑1 (MQ=255)                                                                                                                                     gaaCGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTggg                                                   >  6:304065/1‑111 (MQ=255)                                                                                                                                     gaaCGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTggg                                                   <  5:304065/111‑1 (MQ=255)                                                                                                                                       aCGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTtg           <  3:377032/149‑1 (MQ=255)                                                                                                                                        cttCAAGCAGTTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTcc                            <  2:124733/129‑1 (MQ=255)                                                                                                                                        cttCAAGCAGTTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTcc                            >  1:124733/3‑131 (MQ=255)                                                                                                                                        cttCAAGCAGTTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTtgg          >  2:266243/3‑149 (MQ=255)                                                                                                                                          tCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTtggtg        <  6:190301/149‑1 (MQ=255)                                                                                                                                                aaCTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTTGGTGGTGGcc  >  2:310372/1‑149 (MQ=255)                                                                                                                                                aaCTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTTGGTGGTGGcc  >  3:94316/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CCTGGACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTGAAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTTGGTGGTGGCC  >  NC_002947/4054805‑4055096

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