Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,780,943	C→A	5.7%	P81T (CCC→ACC)	PP_2434 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,780,943	0	C	A	5.7%	223.6 / 2.4	70	P81T (CCC→ACC)	PP_2434	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/31);  new base A (2/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

GTGTATGCCGCCACCAGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCAGA  >  NC_002947/2780796‑2781090                                                                                                                                                    |                                                                                                                                                     gtgtATGCCGCCACCAGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGcc                                                                                                                                                     <  8:167600/149‑1 (MQ=255)       gccgccACCAGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTg                                                                                                                                               >  2:264743/1‑149 (MQ=255)       gccgccACCAGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTg                                                                                                                                               >  7:44867/1‑149 (MQ=255)         cgccACCAGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGTAGCCCGGCTGGt                                                                                                                                             <  3:148491/149‑1 (MQ=255)           ccaccaGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCt                                                                                                                                           >  5:258928/1‑149 (MQ=255)              ccaGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCt                                                                                                                                        <  5:245110/149‑1 (MQ=255)                aGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCa                                                                                                                                      <  6:85685/149‑1 (MQ=255)                          gggCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGtt                                                                                                                             >  4:259258/1‑148 (MQ=255)                            gCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGAc                                                                                                                          <  2:93713/149‑1 (MQ=255)                                   tGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACg                                                                                                                   >  6:191606/1‑149 (MQ=255)                                    ggAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGc                                                                                                                  >  1:179446/1‑149 (MQ=255)                                     gtattgtTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTCGCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCt                                                                                                                 <  8:157270/145‑1 (MQ=255)                                        gtgtTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTggg                                                                                                              <  1:264743/149‑1 (MQ=255)                                           ttGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCaa                                                                                                           >  8:174368/1‑149 (MQ=255)                                           ttGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCaa                                                                                                           >  2:373773/1‑149 (MQ=255)                                              gTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑g                                                                                                       >  5:348602/1‑149 (MQ=255)                                                  aCCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCttt                                                                                                   >  8:276071/1‑149 (MQ=255)                                                   ccTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTc                                                                                                  >  5:290821/1‑149 (MQ=255)                                                    cTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCt                                                                                                 >  2:222875/1‑149 (MQ=255)                                                      ggAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑Gc                                                                                                      >  5:60084/1‑142 (MQ=255)                                                      ggAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑Gc                                                                                                      <  6:60084/142‑1 (MQ=255)                                                       gagaTCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTgg                                                                                               <  6:258928/148‑1 (MQ=255)                                                       gagaTCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGa                                                                                              <  8:49708/149‑1 (MQ=255)                                                        agaTCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGAt                                                                                             <  3:259258/149‑1 (MQ=255)                                                           tCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGAc                                                                                          >  7:69598/1‑149 (MQ=255)                                                               ccccGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGtt                                                                                      <  8:44867/149‑1 (MQ=255)                                                                cccGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTc                                                                                     >  5:100783/1‑149 (MQ=255)                                                                  cGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCgg                                                                                   >  6:380628/1‑149 (MQ=255)                                                                   gAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGCACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTGTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCGTTCTGGATGCCTGTTCgg                                                                                   >  6:315739/1‑148 (MQ=255)                                                                        tttCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATtgc                                                                             >  4:45158/1‑149 (MQ=255)                                                                                ccaccaTGCGCTTGCGTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGCGTTGCCCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCg                                                                     >  8:370510/1‑149 (MQ=255)                                                                                 caccaTGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGATTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCgg                                                                    >  2:317280/1‑149 (MQ=255)                                                                                     aTGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGcc                                                                <  6:100783/149‑1 (MQ=255)                                                                                         gcttgcttCGCATCGAACTCCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTATGCAGCCCGGCTGGTATGCTCAGCCTGACTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCTGCATTGCTGCTGCCGGTGCCAAgc                                                            >  3:273301/1‑149 (MQ=255)                                                                                              cttCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGAt                                                       <  3:45158/149‑1 (MQ=255)                                                                                                  gCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCAt                                                    <  7:174368/148‑1 (MQ=255)                                                                                                  gCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCAt                                                    <  7:276071/148‑1 (MQ=255)                                                                                                  gCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATg                                                   <  1:222875/149‑1 (MQ=255)                                                                                                   cATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGc                                                  <  1:373773/149‑1 (MQ=255)                                                                                                     tCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTTCTGCTGCCGGTGCCAAGCGCGATCATGccc                                                <  6:290821/149‑1 (MQ=255)                                                                                                     tCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGccc                                                <  8:69598/149‑1 (MQ=255)                                                                                                       gAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCa                                               <  2:179446/148‑1 (MQ=255)                                                                                                       gAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAg                                              <  8:331241/149‑1 (MQ=255)                                                                                                          cTGCACGACACGGTTTCCCAGGCGCAGTTGCTCGCTTTGCAGCTCGGCTTGTCTGCTCAGCCTGAGTTGACCCGAACGCTTGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCTGGTGCCAAGCGCGATCATGCCCAGcac                                           <  4:281345/149‑1 (MQ=255)                                                                                                             cccGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGcaccac                                        >  3:347040/1‑149 (MQ=255)                                                                                                              ccGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGcaccac                                        >  6:290548/1‑148 (MQ=255)                                                                                                                     ggTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTAcc                                <  5:191606/149‑1 (MQ=255)                                                                                                                     ggTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTAAACATTGCTTCAGCCGGTGCCAAGCGCGATAATTCCCAGCACCACCAATTAcc                                >  1:323860/1‑149 (MQ=255)                                                                                                                      gTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCt                               <  1:162701/149‑1 (MQ=255)                                                                                                                           ccAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTc                           >  2:196133/1‑148 (MQ=255)                                                                                                                           ccAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTc                           >  6:37502/1‑148 (MQ=255)                                                                                                                            cAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAAtt                                   <  8:284940/139‑1 (MQ=255)                                                                                                                            cAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAAtt                                   >  7:284940/1‑139 (MQ=255)                                                                                                                                   aGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACgg                  <  4:347040/149‑1 (MQ=255)                                                                                                                                     ttGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATCTTCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGc                 >  1:309389/1‑149 (MQ=255)                                                                                                                                       gCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCAcc              >  5:267289/1‑149 (MQ=255)                                                                                                                                        cccGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCACCCCACGGCAccc             >  6:293457/1‑149 (MQ=255)                                                                                                                                        cccGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCAcc              <  1:374106/148‑1 (MQ=255)                                                                                                                                        cccGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCAcc              <  5:37502/148‑1 (MQ=255)                                                                                                                                         ccGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCAACAATTACCTGTTCAACCAACGGCAcc              >  2:374106/1‑147 (MQ=255)                                                                                                                                          cGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCt                               <  5:146063/129‑1 (MQ=255)                                                                                                                                          cGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCt                               >  6:146063/1‑129 (MQ=255)                                                                                                                                           gCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCg            >  2:103345/1‑147 (MQ=255)                                                                                                                                           gCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCg            >  8:337623/1‑147 (MQ=255)                                                                                                                                             tgagCAGACCAGCCGGGCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCAGTGCCAAGCGCGATCATGCCCAGca                                            >  5:227004/4‑113 (MQ=255)                                                                                                                                             tgagCAGACCAGCCGGGCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCAGTGCCAAGCGCGATCATGCCCAGca                                            <  6:232367/110‑1 (MQ=255)                                                                                                                                             tgagCAGACCAGCCGGGCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCAGTGCCAAGCGCGATCATGCCCAGca                                            <  6:227004/110‑1 (MQ=255)                                                                                                                                             tgagCAGACCAGCCGGGCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCAGTGCCAAGCGCGATCATGCCCAGca                                            >  5:232367/4‑113 (MQ=255)                                                                                                                                                gCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGgcg      >  1:111752/1‑148 (MQ=255)                                                                                                                                                 cAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCa    <  3:89375/149‑1 (MQ=255)                                                                                                                                                   gCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCaga  <  5:290548/149‑1 (MQ=255)                                                                                                                                                    cccGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGcaccac                                        >  1:370033/1‑110 (MQ=255)                                                                                                                                                    cccGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGcaccac                                        <  2:370033/110‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GTGTATGCCGCCACCAGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCAGA  >  NC_002947/2780796‑2781090

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