Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,653,193	A→C	10.4%	K12Q (AAG→CAG)	PP_5659 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,653,193	0	A	C	10.4%	196.6 / 13.9	67	K12Q (AAG→CAG)	PP_5659	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base C (4/3);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

GCTTCGACCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTGGCAATCGCGT  >  NC_002947/4653049‑4653336                                                                                                                                                 |                                                                                                                                                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)                                                                                                                                                 |                                                                                                                                                GCTTCGACCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTGGCAATCGCGT  >  NC_002947/4653049‑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: ‑