Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,224,502	T→G	6.1%	intergenic (‑614/+207)	PP_3703 ← / ← PP_3704	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,224,502	0	T	G	6.1%	205.0 / 4.7	66	intergenic (‑614/+207)	PP_3703/PP_3704	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/33);  new base G (2/2);  total (31/35)
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.31e-01

CCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGGA  >  NC_002947/4224358‑4224648                                                                                                                                                   |                                                                                                                                                   ccTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCtgt                                                                                                                                                  >  5:185862/1‑147 (MQ=255) ccTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCtgt                                                                                                                                                  <  6:185862/147‑1 (MQ=255) ccTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGc                                                                                                                                                <  5:83605/149‑1 (MQ=255)  cTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGTGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCa                                                                                                                                               >  4:174753/1‑149 (MQ=255)     tCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCAt                                                                                                                                              <  7:269979/147‑1 (MQ=255)        aTCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTg                                                                                                                                         >  2:349752/1‑149 (MQ=255)          ccAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGt                                                                                                                                       >  2:330052/1‑149 (MQ=255)            aaaCGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCAACCAGCAATGCACAGGTGCACTGCTggg                                                                                                                                       <  7:44630/149‑2 (MQ=255)              aCGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGcc                                                                                                                                   <  5:218880/149‑1 (MQ=255)                   aaTGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACgt                                                                                                                               <  2:308573/148‑1 (MQ=255)                    aTGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGTGAATACCCCGAGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACgtgt                                                                                                                             <  5:258926/149‑1 (MQ=255)                     tGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACgtgt                                                                                                                             >  8:276664/1‑148 (MQ=255)                         ccTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTaaa                                                                                                                          >  7:49978/1‑147 (MQ=255)                         ccTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTaaa                                                                                                                          <  8:49978/147‑1 (MQ=255)                         ccTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGA‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAgg                                                                                                                        >  5:177506/1‑149 (MQ=255)                          cTAAGCCTGGCGGTGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGc                                                                                                                       <  3:174753/149‑1 (MQ=255)                          cTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGc                                                                                                                       >  5:163279/1‑149 (MQ=255)                          cTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGc                                                                                                                       >  6:350996/1‑149 (MQ=255)                           tAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAATGACGGCCAGATACGTGATAAGCGGAGGGTTGCGGGAAATGCCCGCGGAAAATAATCGTAGGTGGC‑‑CGAGCAGTGCTCCTGTGCATTGCTg                                                                                                                                         <  5:348967/130‑1 (MQ=255)                           tAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCt                                                                                                                      >  6:38833/1‑149 (MQ=255)                              gCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTAt                                                                                                                   <  5:38833/149‑1 (MQ=255)                              gCCTGGCGGCGGGGAACAACCTCACTCCCGTAACAAAGACGGCCAGATACATGATAAGCGGAGCGATGCGGGAAAACCCCGAGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTAt                                                                                                                   <  2:314283/149‑1 (MQ=255)                                     ggcggGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTaa                                                                                                            <  2:67305/149‑1 (MQ=255)                                      gcggGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTaaa                                                                                                           <  5:350996/149‑1 (MQ=255)                                            aacaacCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAACCGGAGCGATGCTGGAAAACCCCGAGGAAAATAAACCTAGGTGGC‑‑ACAGCAGAGCACCTGAGCATTCCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAg                                                                                                      >  3:100878/1‑148 (MQ=255)                                                acCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTg                                                                                                 <  6:159115/149‑1 (MQ=255)                                                 ccTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTgg                                                                                                <  2:338511/149‑1 (MQ=255)                                                 ccTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTgg                                                                                                >  8:190469/1‑149 (MQ=255)                                                 ccTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAATGCACAGGTGCACTGCTgg                                                                                                                                        >  3:293588/1‑109 (MQ=255)                                                 ccTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAATGCACAGGTGCACTGCTgg                                                                                                                                        <  4:293588/109‑1 (MQ=255)                                                  cTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTggg                                                                                               >  3:271682/1‑149 (MQ=255)                                                  cTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTggg                                                                                               <  8:140106/149‑1 (MQ=255)                                                    cGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGc                                                                                                                                    <  1:182062/110‑1 (MQ=255)                                                    cGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGc                                                                                                                                    >  2:182062/1‑110 (MQ=255)                                                       tCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAg                                                                                                                >  4:276361/1‑127 (MQ=255)                                                       tCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAg                                                                                                                <  3:276361/127‑1 (MQ=255)                                                         ccGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATcgcg                                                                                        <  8:177701/149‑1 (MQ=255)                                                           gTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGcc                                                                                      <  1:41290/149‑1 (MQ=255)                                                                aaaGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCt                                                                                 >  1:281247/1‑149 (MQ=255)                                                                   gACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTg                                                                                >  3:169565/1‑147 (MQ=255)                                                                            aTACGTGATAAGCTGAGCTATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCGCATCTGTAGCGGAgat                                                                      <  4:100878/148‑1 (MQ=255)                                                                                gTGATAAGCTGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATAc                                                                 <  1:349752/149‑1 (MQ=255)                                                                                gTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAgatgat                                                                   <  7:226177/147‑1 (MQ=255)                                                                                gTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAgatgat                                                                   >  8:226177/1‑147 (MQ=255)                                                                                      aGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAAcc                                                           >  3:177130/1‑149 (MQ=255)                                                                                         ggAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCt                                                          >  7:255913/1‑147 (MQ=255)                                                                                           aGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAgg                                                      <  2:281247/149‑1 (MQ=255)                                                                                                  cGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTc                                                <  7:190469/148‑1 (MQ=255)                                                                                                   gggAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTGAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCa                                              >  5:195612/1‑149 (MQ=255)                                                                                                    ggAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCAt                                             >  7:111286/1‑149 (MQ=255)                                                                                                       aaaCCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATa                                                                  <  5:112727/125‑1 (MQ=255)                                                                                                       aaaCCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATa                                                                  >  6:112727/1‑125 (MQ=255)                                                                                                        aaCCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTaa                                         <  7:276664/149‑1 (MQ=255)                                                                                                              gcgGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTgg                                   <  8:111286/149‑1 (MQ=255)                                                                                                                    aaTAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAgc                             <  8:255913/149‑1 (MQ=255)                                                                                                                     aTAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTa                                                        <  1:281018/121‑1 (MQ=255)                                                                                                                     aTAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTa                                                        >  2:281018/1‑121 (MQ=255)                                                                                                                      tAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGt                                                 >  7:282789/1‑127 (MQ=255)                                                                                                                      tAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGt                                                 <  8:282789/127‑1 (MQ=255)                                                                                                                        aaCCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAgcgtc                          <  8:58606/148‑3 (MQ=255)                                                                                                                            gttctGGC‑‑CCATCATTGCACCTGTGCATTTCTAGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGG‑AAAGCGCCaaggaa                    <  2:333055/145‑1 (MQ=255)                                                                                                                                  gC‑‑CCAGCAGTGCACCTGTGCATTGCTAGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGG‑AAAGCGCCaaggaa                    >  1:333055/1‑143 (MQ=255)                                                                                                                                   c‑‑ccAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTCGGTATCGCGCCCATCTGTAGCGGACATGATACAGAACCTTATGATTGTCCAGTTAAACATGGAAAAGCCCCAAGGAAGGTTCt              >  8:244278/1‑149 (MQ=255)                                                                                                                                      ccAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCtg             >  5:252310/1‑149 (MQ=255)                                                                                                                                      ccAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAACACAAGACTTGGGATCGCGCCCATCTGTATCGGAGATGATACAGAAACTTATGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAATGTtgtg             >  1:101735/1‑149 (MQ=255)                                                                                                                                      ccAGCAATGCACAGGTGCACTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCtg             >  5:136261/1‑149 (MQ=255)                                                                                                                                        agcagTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCtgtg           <  1:330052/149‑1 (MQ=255)                                                                                                                                                  cTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGGa  >  2:315947/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGGA  >  NC_002947/4224358‑4224648

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