Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,653,190	G→T	10.3%	G11C (GGC→TGC)	PP_5659 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,653,190	0	G	T	10.3%	206.5 / 11.9	68	G11C (GGC→TGC)	PP_5659	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base T (4/3);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

AAGGCTTCGACCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTGGCAATCGCGT  >  NC_002947/4653046‑4653336                                                                                                                                                 |                                                                                                                                                   aaGGCTTCGACCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGc                                                                                                                                                  >  2:273486/1‑147 (MQ=255)    gCTTCGACCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAggg                                                                                                                                             >  2:87010/1‑149 (MQ=255)        cGACCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTTGCCAGGGccc                                                                                                                                          <  7:128748/148‑2 (MQ=255)          aCCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCg                                                                                                                                       <  2:264154/149‑1 (MQ=255)                 gCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGgcgcg                                                                                                                                >  4:302211/1‑149 (MQ=255)                   cTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCt                                                                                                                              >  4:158597/1‑149 (MQ=255)                      gAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTc                                                                                                                             <  1:87010/147‑1 (MQ=255)                        aaCTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCtt                                                                                                                          <  6:168018/148‑1 (MQ=255)                                 ccAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAgg                                                                                                                                              >  8:245834/1‑119 (MQ=255)                                 ccAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAgg                                                                                                                                              <  7:245834/119‑1 (MQ=255)                                        gCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGgcgc                                                                                                          <  1:143683/148‑1 (MQ=255)                                        gCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGcc                                                                                                         <  3:25016/149‑1 (MQ=255)                                        gCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGcc                                                                                                         >  3:298848/1‑149 (MQ=255)                                         cATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCt                                                                                                        >  1:190703/1‑149 (MQ=255)                                          aTTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTc                                                                                                       <  3:302211/149‑1 (MQ=255)                                                ccAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTc                                                                                                  >  4:101113/1‑148 (MQ=255)                                                 cAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCta                                                                                                <  7:151149/149‑1 (MQ=255)                                                 cAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCta                                                                                                >  7:767/1‑149 (MQ=255)                                                 cAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCta                                                                                                >  8:102376/1‑149 (MQ=255)                                                 cAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCta                                                                                                >  4:293830/1‑149 (MQ=255)                                                    aGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATAc                                                                                             <  3:158597/149‑1 (MQ=255)                                                          agtaAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGccc                                                                                        <  4:298848/148‑1 (MQ=255)                                                                   aaaaaaGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGACCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGCGGCGCCTCGTGTCTATACGGCCCCCAAGGCGGt                                                                              >  6:236226/1‑149 (MQ=255)                                                                    aaaaaGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGtt                                                                             <  1:273486/149‑1 (MQ=255)                                                                     aaaaGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACaa                                                                                    >  3:242457/1‑141 (MQ=255)                                                                     aaaaGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACaa                                                                                    <  4:242457/141‑1 (MQ=255)                                                                        aGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGAc                                                                         <  7:102376/149‑1 (MQ=255)                                                                         gCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACg                                                                        <  2:190703/149‑1 (MQ=255)                                                                             tACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCa                                                                                                                  <  4:90537/103‑1 (MQ=255)                                                                             tACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCa                                                                                                                  >  3:90537/1‑103 (MQ=255)                                                                               cAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGAtt                                                                   <  2:100318/148‑1 (MQ=255)                                                                               cAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTg                                                                  >  2:313317/1‑149 (MQ=255)                                                                                 aaaCATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCTGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGc                                                                 >  3:28079/1‑148 (MQ=255)                                                                                  aaCATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGa                                                               >  3:60308/1‑149 (MQ=255)                                                                                               gTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGACTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACgtagaa                                                  <  5:225042/149‑1 (MQ=255)                                                                                                    tttCCCAAGGGCATGATGGCCTCGCTCCTGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTg                                              <  4:28079/148‑1 (MQ=255)                                                                                                        ccAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTccc                                          >  3:245884/1‑148 (MQ=255)                                                                                                            gggCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAgg                                      >  3:131482/1‑148 (MQ=255)                                                                                                             ggCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGa                                    >  2:57820/1‑149 (MQ=255)                                                                                                              gCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGAc                                   <  1:313317/149‑1 (MQ=255)                                                                                                               cATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGtt                                               >  7:298472/1‑136 (MQ=255)                                                                                                               cATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGtt                                               <  8:298472/136‑1 (MQ=255)                                                                                                                atgatgGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACg                                  >  4:99072/1‑148 (MQ=255)                                                                                                                   atgGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGaa                               <  3:99072/148‑1 (MQ=255)                                                                                                                     ggCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAgcg                            >  3:179759/1‑149 (MQ=255)                                                                                                                     ggCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAgc                             <  4:131482/148‑1 (MQ=255)                                                                                                                     ggCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGACCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAgc                             >  5:270268/1‑148 (MQ=255)                                                                                                                     ggCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGACCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAgc                             <  6:270268/148‑1 (MQ=255)                                                                                                                      gCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAgcg                            <  1:57820/148‑1 (MQ=255)                                                                                                                        cTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCAt                         >  5:234258/1‑149 (MQ=255)                                                                                                                            cTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGAGTGTGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCAcc                     <  5:332779/149‑1 (MQ=255)                                                                                                                               ccGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCg                    <  4:60308/147‑1 (MQ=255)                                                                                                                                cGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAg                 <  3:101113/149‑1 (MQ=255)                                                                                                                                  cTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTa               >  6:103533/1‑149 (MQ=255)                                                                                                                                     cccgaaggCCCTTGCCAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTg            >  7:271983/7‑149 (MQ=255)                                                                                                                                      ccgaaggCCCTTGCCAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTg            >  6:139081/6‑148 (MQ=255)                                                                                                                                      ccgaaggCCCTTGCCAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTg            <  5:139081/143‑1 (MQ=255)                                                                                                                                       cgaaggCCCTTGCCAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTg            <  5:270815/143‑1 (MQ=255)                                                                                                                                       cgaaggCCCTTGCCAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTg            >  6:270815/5‑147 (MQ=255)                                                                                                                                        gaaggCCCTTGCCAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGAGCATCCCCGCAGTACTGGc          >  5:263266/4‑148 (MQ=255)                                                                                                                                        aGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTGCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTGGc          >  1:295538/1‑148 (MQ=255)                                                                                                                                        aGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTGGc          >  5:249411/1‑148 (MQ=255)                                                                                                                                               cTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAgg                                      <  6:331761/113‑1 (MQ=255)                                                                                                                                               cTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAgg                                      >  5:331761/1‑113 (MQ=255)                                                                                                                                               cTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAgcgc                           <  5:244470/124‑1 (MQ=255)                                                                                                                                               cTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAgcgc                           >  6:244470/1‑124 (MQ=255)                                                                                                                                               cTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTGGCAATCGCGt  <  4:245884/149‑1 (MQ=255)                                                                                                                                               cTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTGGCAATCGCGt  <  3:293830/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   AAGGCTTCGACCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTGGCAATCGCGT  >  NC_002947/4653046‑4653336

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