Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,181,146	A→T	5.4%	L266H (CTC→CAC)	PP_4560 ←	ribonuclease BN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,181,146	0	A	T	5.4%	230.1 / 2.3	74	L266H (CTC→CAC)	PP_4560	ribonuclease BN
Reads supporting (aligned to +/- strand):  ref base A (35/35);  new base T (2/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GACGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGC  >  NC_002947/5181001‑5181288                                                                                                                                                   |                                                                                                                                               gACGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGag                                                                                                                                               <  5:65645/147‑1 (MQ=255) gACGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGag                                                                                                                                               >  5:103434/1‑147 (MQ=255)        tGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCAc                                                                                                                                      <  7:99913/149‑1 (MQ=255)           ggCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGc                                                                                                                                   >  2:214840/1‑149 (MQ=255)              ttGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCTGCCGTGCTGCTGCTCGGcgc                                                                                                                                 >  7:102804/1‑147 (MQ=255)                gCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCAATTCGGC‑GCCGAGCAGCAGCACGGCAGcc                                                                                                                               <  1:214840/148‑1 (MQ=255)                  tCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGa                                                                                                                             >  6:308191/1‑148 (MQ=255)                   cGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGaaa                                                                                                                           >  7:163389/1‑149 (MQ=255)                    gggTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGaaa                                                                                                                           <  6:103434/148‑1 (MQ=255)                       tttGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCTGCCGTGCTGCTGCTCGGC‑GCCGAAatgaa                                                                                                                       <  7:57282/149‑4 (MQ=38)                         tGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAg                                                                                                                      <  5:308191/148‑1 (MQ=255)                         tGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTACCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAg                                                                                                                      <  2:321725/148‑1 (MQ=255)                          gggTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCGGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCTGCCGTGCTGCTGCTCGGC‑GCCGAAatgaatg                                                                                                                     <  8:102804/148‑6 (MQ=37)                                  gTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagca                                                                                                            <  4:20302/149‑1 (MQ=255)                                  gTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagc                                                                                                             <  8:143587/148‑1 (MQ=255)                                     aCCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagcaaca                                                                                                         <  3:347159/149‑1 (MQ=255)                                      cccAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagcaacag                                                                                                        >  8:205369/1‑149 (MQ=255)                                       ccAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagcaacagc                                                                                                       >  7:345296/1‑149 (MQ=255)                                       ccAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagcaacag                                                                                                        >  3:105491/1‑148 (MQ=255)                                         aGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCAc                                                                                                     <  4:105491/149‑1 (MQ=255)                                                tCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGAtt                                                                                              >  1:223394/1‑149 (MQ=255)                                                  aaGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGaaa                                                                                                                   >  3:336905/1‑126 (MQ=255)                                                  aaGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTg                                                                                             <  1:330987/148‑1 (MQ=255)                                                    ggTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCAc                                                                                          <  2:205655/149‑1 (MQ=255)                                                         ccTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAAtg                                                                                     >  2:220842/1‑149 (MQ=255)                                                          cTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAAtgc                                                                                    >  6:11744/1‑149 (MQ=255)                                                           ttGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAg                                                                                                                      <  2:7409/114‑1 (MQ=255)                                                           ttGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAg                                                                                                                      >  1:7409/1‑114 (MQ=255)                                                                 ttttCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGTAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTa                                                                             >  5:479/1‑149 (MQ=255)                                                                   ttCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTa                                                                             <  3:346795/147‑1 (MQ=255)                                                                     ccccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACAt                                                                          <  2:223394/148‑1 (MQ=255)                                                                        cGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCg                                                                      >  2:268230/1‑149 (MQ=255)                                                                          tCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATaca                                                                                                               <  6:84718/106‑1 (MQ=255)                                                                          tCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATaca                                                                                                               >  5:84718/1‑106 (MQ=255)                                                                          tCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGtt                                                                    <  4:58512/149‑1 (MQ=255)                                                                          tCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGtt                                                                    <  7:205369/149‑1 (MQ=255)                                                                           cGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTg                                                                   <  8:345296/149‑1 (MQ=255)                                                                            gAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGATATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTg                                                                   >  6:29230/1‑148 (MQ=255)                                                                            gAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAAtgctg                                                                                  <  8:285419/133‑1 (MQ=255)                                                                            gAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAAtgctg                                                                                  >  7:285419/1‑133 (MQ=255)                                                                            gAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTg                                                                   >  2:332228/1‑148 (MQ=255)                                                                             aaGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagcaacag                                                                                                        >  5:199658/1‑110 (MQ=255)                                                                             aaGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagcaacag                                                                                                        <  6:199658/110‑1 (MQ=255)                                                                              aGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAg                                                                <  8:163389/149‑1 (MQ=255)                                                                               gTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCTCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGt                                                               <  5:220256/149‑1 (MQ=255)                                                                                  tttCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTg                                                                   <  7:69442/142‑1 (MQ=255)                                                                                  tttCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCgg                                                            >  3:93732/1‑149 (MQ=255)                                                                                  tttCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACAGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTg                                                                   >  8:69442/1‑142 (MQ=255)                                                                                    tctcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGTAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCgg                                                            <  6:479/147‑1 (MQ=255)                                                                                    tctcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCgg                                                            <  1:332228/147‑1 (MQ=255)                                                                                       cGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAg                                                       >  6:261121/1‑149 (MQ=255)                                                                                        gCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGTTGCCGTACATGGCGTTGTAGTCGGCAAAg                                                       <  1:220842/148‑1 (MQ=255)                                                                                        gCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAgg                                                      >  7:279197/1‑149 (MQ=255)                                                                                        gCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAg                                                       <  6:53210/148‑1 (MQ=255)                                                                                               ttCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGtt                                                                    >  1:200261/1‑128 (MQ=255)                                                                                               ttCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGtt                                                                    <  2:200261/128‑1 (MQ=255)                                                                                                  ttGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAAt                                                                                      <  4:241523/107‑1 (MQ=255)                                                                                                  ttGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAAt                                                                                      >  3:241523/1‑107 (MQ=255)                                                                                                           gCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCa                                     >  4:257920/1‑147 (MQ=255)                                                                                                                 aTGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAg                             <  5:261121/149‑1 (MQ=255)                                                                                                                  tgtgCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAgg                            <  5:11744/149‑1 (MQ=255)                                                                                                                      cTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAg                         <  8:279197/148‑1 (MQ=255)                                                                                                                        cGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGACCAGGGAGGc                       >  2:75394/1‑148 (MQ=255)                                                                                                                         gATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCa                     >  6:81318/1‑149 (MQ=255)                                                                                                                         gATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCGTCACGTTGTAGGCAAAGCCCAGTGAGGCCa                     >  8:172992/1‑149 (MQ=255)                                                                                                                         gATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGAAAGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCa                     >  8:40548/1‑149 (MQ=255)                                                                                                                                  attcattTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTa                                                                 >  2:41452/1‑96 (MQ=255)                                                                                                                                  attcattTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTa                                                                 <  1:41452/96‑1 (MQ=255)                                                                                                                                  attcattTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGTTCTTCACGTAGTAGGCAAAGCCCAGGGAg                         <  8:352626/136‑1 (MQ=255)                                                                                                                                  attcattTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAg                         >  7:352626/1‑136 (MQ=255)                                                                                                                                  attcattTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccaca            >  8:120627/1‑149 (MQ=255)                                                                                                                                   ttcattTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacac           >  4:263281/1‑149 (MQ=255)                                                                                                                                          cGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAAcc    >  5:119449/1‑149 (MQ=255)                                                                                                                                           ggC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTATTATTCAAAGCCCAGGTAGTCCACGATCTACACCACAAcc    <  2:279527/148‑1 (MQ=255)                                                                                                                                            gctgcCGTGCTGCTGCTCGGC‑GCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGc  >  4:176783/4‑149 (MQ=255)                                                                                                                                            gc‑gcCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGc  >  2:107832/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                               GACGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGC  >  NC_002947/5181001‑5181288

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