Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,185,560	C→G	7.8%	L446L (CTC→CTG)	PP_3681 →	putative Helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,185,560	0	C	G	7.8%	236.2 / 9.4	77	L446L (CTC→CTG)	PP_3681	putative Helicase
Reads supporting (aligned to +/- strand):  ref base C (38/33);  new base G (3/3);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

AGGACTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTGGGGCAA  >  NC_002947/4185421‑4185703                                                                                                                                            |                                                                                                                                                aGGACTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTGAGCAGCgc                                                                                                                                         >  6:87053/1‑147 (MQ=255) aGGACTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTGAGCAGCgc                                                                                                                                         <  6:163288/148‑2 (MQ=255) aGGACTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTGAGCAGCgc                                                                                                                                         <  5:87053/148‑2 (MQ=255)   gACTGGTGGAAGCTCACTAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGaaaa                                                                                                                                      >  6:67846/1‑149 (MQ=255)           gAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTc                                                                                                                              >  7:36646/1‑149 (MQ=255)            aaGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTcc                                                                                                                             >  7:291487/1‑149 (MQ=255)             aGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTAAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTccc                                                                                                                            >  1:93443/1‑149 (MQ=255)             aGCTCACAAGCAGCTTCATCCGTGATGTAGGATATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTccc                                                                                                                            >  1:220067/1‑149 (MQ=255)                    aaGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGccc                                                                                                                     >  2:254864/1‑149 (MQ=255)                      gcagcTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCtt                                                                                                                   >  7:278674/1‑149 (MQ=255)                          cTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGccc                                                                                                               <  8:291487/149‑1 (MQ=255)                             cATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCgg                                                                                                             <  2:151870/148‑1 (MQ=255)                               tCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCtt                                                                                                          >  6:153288/1‑149 (MQ=255)                               tCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCtt                                                                                                          >  3:285073/1‑149 (MQ=255)                               tCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCtt                                                                                                          >  3:280162/1‑149 (MQ=255)                                 cGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTc                                                                                                        <  2:220067/149‑1 (MQ=255)                                   tGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAg                                                                                                      <  8:36646/149‑1 (MQ=255)                                     aTGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGa                                                                                                     <  2:93443/148‑1 (MQ=255)                                        tAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACAATGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGAt                                                                                                 >  2:259010/1‑149 (MQ=255)                                               tCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTAAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCg                                                                                          >  2:135167/1‑149 (MQ=255)                                                ccAATCCTGACTGCGCACTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGACCGGCTTTCAGACGATCAGGCCg                                                                                          >  3:50587/1‑148 (MQ=255)                                                 cAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGc                                                                                                                 >  1:60181/1‑124 (MQ=255)                                                 cAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGc                                                                                                                 <  2:60181/124‑1 (MQ=255)                                                  aaTCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGa                                                                                                  >  5:70909/1‑138 (MQ=255)                                                  aaTCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGa                                                                                                  <  6:70909/138‑1 (MQ=255)                                                  aaTCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGAt                                                                                                 <  3:237642/139‑1 (MQ=255)                                                  aaTCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGAt                                                                                                 >  4:237642/1‑139 (MQ=255)                                                  aaTCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTAc                                                                                       >  6:50818/1‑149 (MQ=255)                                                      cTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCt                                                                                    <  4:285073/148‑1 (MQ=255)                                                      cTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCt                                                                                    <  4:280162/148‑1 (MQ=255)                                                      cTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCACGGCTTTAAGACGATCAGGCCGTACGCtt                                                                                   >  4:8818/1‑149 (MQ=255)                                                         aCTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTtc                                                                                 <  5:350835/148‑1 (MQ=255)                                                             cgcTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTcac                                                                             >  1:62721/1‑148 (MQ=255)                                                              gcTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACaa                                                                           <  6:39234/149‑1 (MQ=255)                                                                tctcAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCttt                                                                                  >  3:334548/1‑140 (MQ=255)                                                                tctcAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCttt                                                                                  <  4:334548/140‑1 (MQ=255)                                                                      gTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACaa                                                                    >  3:110121/1‑148 (MQ=255)                                                                             cGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCaa                                                              <  8:278674/147‑1 (MQ=255)                                                                              gAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACCATCCACAAAAGCAAggg                                                           <  4:358568/149‑1 (MQ=255)                                                                              gAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAggg                                                           >  1:360336/1‑149 (MQ=255)                                                                                 tCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGGACGCTTTCTCACAATCCCCAACAGCAAGGGGCt                                                        >  6:284886/1‑149 (MQ=255)                                                                                   gAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTcac                                                                             <  1:184894/126‑1 (MQ=255)                                                                                   gAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTcac                                                                             >  2:184894/1‑126 (MQ=255)                                                                                   gAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCtt                                                       >  4:233148/1‑148 (MQ=255)                                                                                   gAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCtt                                                       >  6:147104/1‑148 (MQ=255)                                                                                   gAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTCGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCACTTGCCCGGCTTTCAGACGATCAGGCCGTACCCTTTCTCACAATCCACAAAAGCAAGGGGCtt                                                       >  8:248654/1‑148 (MQ=255)                                                                                    aaaCTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGa                                                     <  1:259010/149‑1 (MQ=255)                                                                                          cGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTca                                               <  5:147104/149‑1 (MQ=255)                                                                                             ttGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGaaaa                                                                                                                                      >  7:248195/1‑59 (MQ=255)                                                                                               gCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTTCGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACg                                           <  1:254864/148‑1 (MQ=255)                                                                                               gCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACgg                                          >  8:303138/1‑149 (MQ=255)                                                                                                cACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACgg                                          <  5:67846/148‑1 (MQ=255)                                                                                                   gAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTca                                               <  5:152774/140‑1 (MQ=255)                                                                                                   gAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTca                                               >  6:152774/1‑140 (MQ=255)                                                                                                   gAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGtca                                       >  3:222967/1‑148 (MQ=255)                                                                                                    aaGTCGTTCGTGAGGTCAGAGATCGATTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGtcatc                                     <  4:50587/149‑1 (MQ=255)                                                                                                      gTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTcatcat                                   <  1:135167/149‑1 (MQ=255)                                                                                                             gtGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCgg                             <  2:360336/148‑1 (MQ=255)                                                                                                                aGGTCAGAGATCGACTTGGAGCCCTGCTCAGAAGCGAAAAAACCATTCCCGGTGCCCTTGCCCGGCATTCAGAACATCAGGCCCTACGCAAACACACAAACCACAAAAGCAAGGGGCTTGAGTTTCACAACGTCATTATTATCGGGGt                          >  2:323008/1‑148 (MQ=255)                                                                                                                  gTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCggg                            <  4:183675/144‑1 (MQ=255)                                                                                                                  gTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCggg                            >  3:183675/1‑144 (MQ=255)                                                                                                                     agagATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGaaa                     <  5:50818/148‑1 (MQ=255)                                                                                                                     agagATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGaaa                     <  4:222967/148‑1 (MQ=255)                                                                                                                     agagATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAAc                    >  4:345166/1‑149 (MQ=255)                                                                                                                     agagATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAAc                    >  8:277724/1‑149 (MQ=255)                                                                                                                         aTCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCaa                <  3:233148/149‑1 (MQ=255)                                                                                                                                cttttcgcTGCTGAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGaaa                     <  4:71970/132‑1 (MQ=255)                                                                                                                                cttttcgcTGCTGAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGaaa                     >  3:71970/6‑137 (MQ=255)                                                                                                                                 ttttcgcTGCTGAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAACGCAAGGGGCTTGAGTTTCACACGGGCATCATTCTCGGGGTCGAAGCCCAAGCTTTcgg        >  4:257047/5‑147 (MQ=255)                                                                                                                                  gAGCGTTGATCAGCATCGAAAAGACCATTCACGGTGCCCTTTCCCGGCTTTCAGACGATCAGGACGTACGCTTTCTCTCAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTg        <  7:248654/148‑1 (MQ=255)                                                                                                                                  gAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAAt                                                                          >  4:276160/1‑82 (MQ=255)                                                                                                                                  gAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAAt                                                                          <  3:276160/82‑1 (MQ=255)                                                                                                                                  gAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTgg       <  7:277724/149‑1 (MQ=255)                                                                                                                                   aGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTggg      <  7:303138/149‑1 (MQ=255)                                                                                                                                       ctgctCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTGGGGCaa  >  1:191975/1‑149 (MQ=255)                                                                                                                                          ctCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAAcc                   <  2:163237/129‑1 (MQ=255)                                                                                                                                          ctCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAAcc                   >  1:163237/1‑129 (MQ=255)                                                                                                                                            |                                                                                                                                                AGGACTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTGGGGCAA  >  NC_002947/4185421‑4185703

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: ‑