Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	104,876	G→C	7.0%	P345P (CCG→CCC)	PP_0099 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	104,876	0	G	C	7.0%	217.4 / 4.2	69	P345P (CCG→CCC)	PP_0099	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/33);  new base C (3/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

TATGACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAGTTCCCGCCTTT  >  NC_002947/104733‑105022                                                                                                                                                |                                                                                                                                                     tATGACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTgg                                                                                                                                                 <  1:160273/149‑1 (MQ=255)     aCATCCGTTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCa                                                                                                                                              <  7:164366/148‑1 (MQ=255)     aCAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTggg                                                                                                                                                >  2:239143/1‑146 (MQ=255)        gCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATg                                                                                                                                          >  5:22825/1‑149 (MQ=255)           ggTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGgcg                                                                                                                                       <  7:12414/149‑1 (MQ=255)                  cGTGCGCATGCGTTGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCt                                                                                                                                <  1:199546/149‑1 (MQ=255)                  cGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGTAGCGGGATTGCTCGGCGGTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCt                                                                                                                                <  1:81823/149‑1 (MQ=255)                         aTGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATg                                                                                                                          >  8:298891/1‑148 (MQ=255)                                   aaGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAacc‑cca‑ccc                                                                                                              <  7:298891/148‑1 (MQ=255)                                   aaGGAGTATTTGCCCAGTTGCTGGCGCAGCGCATTGCCAAGGCAATGAAACGGCTGGAGTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGTCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAacc‑cca‑cccc                                                                                                             <  6:196086/149‑1 (MQ=255)                                    aGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑ccccc                                                                                                            >  6:224740/1‑149 (MQ=255)                                     ggAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑cccccc                                                                                                           <  6:22825/149‑1 (MQ=255)                                      gAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCca                                                                                                          >  1:203455/1‑149 (MQ=255)                                         tATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAg                                                                                                       >  4:14868/1‑149 (MQ=255)                                          aTTTGCCGAGTTGCTGGAGAAGCGCTTTGCGAAGGCAATGAAACGGCTGGATTTCGAAGGGCGGGAGGCGCAGGAGCTGGATTTCTCGGCGTTTTGCCCGCCGGGTCGGCAGATGGCGCTGTTCTGACATGACCC‑CCC‑ACCCCCACAgg                                                                                                      >  8:149932/1‑149 (MQ=255)                                           tttGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAgg                                                                                                      <  2:203455/148‑1 (MQ=255)                                                 gAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTgc                                                                                               <  5:59223/149‑1 (MQ=255)                                                 gAGGTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTgc                                                                                               <  4:223222/149‑1 (MQ=255)                                                      gCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCAtt                                                                                          >  7:110729/1‑149 (MQ=255)                                                      gCTGGAGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCAtt                                                                                          >  2:43730/1‑149 (MQ=255)                                                       cTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCACTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCAtt                                                                                          <  6:203052/148‑1 (MQ=255)                                                        tGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGc                                                                                        >  7:119930/1‑149 (MQ=255)                                                          gcgcAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGc                                                                                                                                               >  5:113084/1‑94 (MQ=255)                                                          gcgcAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGc                                                                                                                                               <  6:113084/94‑1 (MQ=255)                                                               gcgcTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCTGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAAc                                                                                 >  7:214033/1‑149 (MQ=255)                                                                   cttcccAAGGCAATCAAACGGCTGGCTTTCGAAGGGCGCCAGGCGCACGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGATCACAACGCaa                                                                             <  8:50442/145‑1 (MQ=255)                                                                     tGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAAc                                                                            <  1:43730/148‑1 (MQ=255)                                                                      gCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACt                                                                           >  3:248596/1‑148 (MQ=255)                                                                         aaGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACtaat                                                                        >  5:276925/1‑148 (MQ=255)                                                                         aaGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACtaat                                                                        >  1:312311/1‑148 (MQ=255)                                                                         aaGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACtaat                                                                        >  1:291084/1‑148 (MQ=255)                                                                          aGGCAATGAAACGGGTGGATTTCGAAGGGCGCGAGGGGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACtaataa                                                                      <  5:312482/149‑1 (MQ=255)                                                                          aGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACtaataa                                                                      <  2:254568/149‑1 (MQ=255)                                                                           ggCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCTCCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACtaataa                                                                      <  1:239143/149‑1 (MQ=255)                                                                            gCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATaac                                                                     >  4:231145/1‑148 (MQ=255)                                                                              aaTGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATaacaac                                                                  <  5:224740/149‑1 (MQ=255)                                                                              aaTGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATaacaac                                                                  >  6:28987/1‑149 (MQ=255)                                                                                tGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAAcc                                                                 >  1:248486/1‑148 (MQ=255)                                                                                    aCGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCaataa                                                            <  8:119930/149‑1 (MQ=255)                                                                                       gCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTg                                                         <  3:231145/149‑1 (MQ=255)                                                                                       gCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTg                                                         >  2:66022/1‑149 (MQ=255)                                                                                        cTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTgc                                                        <  8:151188/149‑1 (MQ=255)                                                                                        cTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTgc                                                        <  5:28987/149‑1 (MQ=255)                                                                                         tGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTgcg                                                       >  1:96187/1‑149 (MQ=255)                                                                                           gATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTgcgc                                                      <  2:312311/148‑1 (MQ=255)                                                                                              ttCGATGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCCGGTGGGc                                                                                                                                               <  8:204679/58‑1 (MQ=255)                                                                                              ttCGAAGGGCGCGAGGCGCAGGTAATGGATTGCTCGGCGTTTTGCCCACACGGTGGGc                                                                                                                                               >  7:204679/1‑58 (MQ=255)                                                                                                cGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGTGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATa                                                 >  5:125638/1‑148 (MQ=255)                                                                                                cGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCACCCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTgcgc                                                      >  8:108885/1‑144 (MQ=255)                                                                                                cGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCACCCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTgcgc                                                      <  7:108885/144‑1 (MQ=255)                                                                                                 gAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACa                                               >  2:71731/1‑149 (MQ=255)                                                                                                       cgcgAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTgg                                         <  7:54/149‑1 (MQ=255)                                                                                                           aGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGTGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAAc                                     <  6:125638/149‑1 (MQ=255)                                                                                                               gcAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGAt                                 <  6:276925/149‑1 (MQ=255)                                                                                                               gcAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGAt                                 <  2:291084/149‑1 (MQ=255)                                                                                                               gcAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCAGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGAt                                 >  5:148365/1‑149 (MQ=255)                                                                                                                cAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATg                                >  6:261915/1‑149 (MQ=255)                                                                                                                         gATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGc                       <  3:14868/149‑1 (MQ=255)                                                                                                                               tCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCACCCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTCGCAACTGATGTTTTTTTGCtatta                  >  8:137300/1‑149 (MQ=255)                                                                                                                               tCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTgcgc                                                      >  7:71179/1‑112 (MQ=255)                                                                                                                               tCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTgcgc                                                      <  8:71179/112‑1 (MQ=255)                                                                                                                               tCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCtattat                 >  4:72196/1‑149 (MQ=255)                                                                                                                                  gccatctGCCCACCCGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGGAACTGATGTTTTTTTGCTATTATCCa              >  7:270984/7‑149 (MQ=255)                                                                                                                                     ttttGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAGtt           >  1:249888/1‑149 (MQ=255)                                                                                                                                         gCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAGTTccc        <  2:249888/148‑1 (MQ=255)                                                                                                                                          cccACCCGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATaacaac                                                                  <  1:139930/89‑1 (MQ=255)                                                                                                                                          cccACCCGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATaacaac                                                                  >  2:139930/1‑89 (MQ=255)                                                                                                                                            cACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAGTTCCCGCCt    >  5:4852/1‑149 (MQ=255)                                                                                                                                              ccGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAGTTCCCGCCttt  <  6:4852/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     TATGACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACC‑CCA‑CCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAGTTCCCGCCTTT  >  NC_002947/104733‑105022

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