Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,388,818	A→G	8.2%	H287H (CAT→CAC) ‡	PP_2991 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,388,818	0	A	G	8.2%	208.6 / 7.1	73	H287H (CAT→CAC) ‡	PP_2991	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base G (2/4);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GAGTGGTTTGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGAAGGACACGTCGACTTTTCCGGTCTGG  >  NC_002947/3388671‑3388965                                                                                                                                                    |                                                                                                                                                    gAGTGGTTTGGGTGCTTAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAAATCGACAGAGCGCCGTCGCCTAGTGTGAGGNTGGCTGAATCTGCGTCGCGGCTCGAGCCAt                                                                                                                                                    <  1:34970/149‑1 (MQ=255)    tGGTTTGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGAc                                                                                                                                                 <  2:367472/149‑1 (MQ=255)      gTTTGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCGCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACg                                                                                                                                                >  5:280545/1‑148 (MQ=255)       tttGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTAGTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACg                                                                                                                                                >  1:259193/1‑147 (MQ=255)       tttGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACg                                                                                                                                                >  3:154451/1‑147 (MQ=255)       tctGGGGGCAGAGAGCGCACTCTGAACCGCTGGCACGTCTTCGTACTTCCTCATTTTCCCCAACTTCAGGGGGCGGCGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCGGAATCTGCGGCATGGCTCGAGCCGCGACGCa                                                                                                                                              <  4:298529/147‑1 (MQ=255)          gggTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCCACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCATGGCTCGAGCCGCGACGCAga                                                                                                                                            <  1:14717/148‑2 (MQ=255)             tGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGc                                                                                                                                        <  4:167106/149‑1 (MQ=255)             tGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGc                                                                                                                                        <  7:40556/149‑1 (MQ=255)                aTAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCATGGCTCGAGCCGCGACGCa                                                                                                                                              <  7:16523/140‑1 (MQ=255)                aTAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCATGGCTCGAGCCGCGACGCa                                                                                                                                              >  8:16523/1‑140 (MQ=255)                 tAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCggg                                                                                                                                     >  3:109147/1‑148 (MQ=255)                    agCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAAc                                                                                                                                 >  7:144585/1‑149 (MQ=255)                    agCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAAc                                                                                                                                 >  7:31389/1‑149 (MQ=255)                      cgcACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCg                                                                                                                               <  4:109147/149‑1 (MQ=255)                       gcACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGc                                                                                                                              <  8:31389/149‑1 (MQ=255)                         aCTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCAt                                                                                                                            <  8:144585/149‑1 (MQ=255)                           tctGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATaa                                                                                                                          >  3:332552/1‑149 (MQ=255)                           tctGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCATGGCTCGAGCCGCGACGCa                                                                                                                                              >  7:228380/1‑129 (MQ=255)                           tctGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCATGGCTCGAGCCGCGACGCa                                                                                                                                              <  8:228380/129‑1 (MQ=255)                              gAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCg                                                                                                                       >  4:176503/1‑149 (MQ=255)                                aaCGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTa                                                                                                                     <  8:238637/149‑1 (MQ=255)                                aaCGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCACGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTa                                                                                                                     <  6:280545/149‑1 (MQ=255)                                                     ttCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCTTATGCCAACTGGCCCTGAGGATc                                                                                                >  8:248216/1‑149 (MQ=255)                                                       catcatTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATc                                                                                                >  4:230736/1‑147 (MQ=255)                                                        atcatTTTCCCGAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGATCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCga                                                                                              <  3:302139/148‑1 (MQ=255)                                                        atcatTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGAGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGccc                                                                                                        <  4:19512/138‑1 (MQ=255)                                                        atcatTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGAGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGccc                                                                                                        >  3:19512/1‑138 (MQ=255)                                                            ttttCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGAt                                                                                           <  7:340195/147‑1 (MQ=255)                                                            ttttACCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGTGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAAATGGCCCTGATTATCgagaa                                                                                           >  8:340195/1‑146 (MQ=255)                                                            attaccccCACTTCAGGGTCGGGAGGTAAGTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCGATGACGCAGTTGGCGGGAAACCGCATAACCGTATCCCAACTGGCCCTGATGATCGAGATTg                                                                                         <  2:315261/145‑1 (MQ=255)                                                                  ccAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAg                                                                                   >  2:298593/1‑149 (MQ=255)                                                                      cTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAgctgc                                                                               >  3:30638/1‑149 (MQ=255)                                                                      cTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAgctgc                                                                               <  2:279579/149‑1 (MQ=255)                                                                         cAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCt                                                                                                       <  4:153351/122‑1 (MQ=255)                                                                         cAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCt                                                                                                       >  3:153351/1‑122 (MQ=255)                                                                           gggTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGt                                                                           >  5:333625/1‑148 (MQ=255)                                                                            ggTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTTGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGa                                                                         <  8:84256/149‑1 (MQ=255)                                                                            ggTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCtg                                                                                                      >  1:288293/1‑120 (MQ=255)                                                                            ggTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCtg                                                                                                      <  2:288293/120‑1 (MQ=255)                                                                             gTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGAc                                                                        >  3:318737/1‑149 (MQ=255)                                                                                cGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGa                                                                      >  5:345301/1‑148 (MQ=255)                                                                                 ggAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGata                                                                    >  6:225316/1‑149 (MQ=255)                                                                                   aGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATc                                                                  <  6:345301/149‑1 (MQ=255)                                                                                     gTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGa                                                                <  2:259193/149‑1 (MQ=255)                                                                                         ttCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTc                                                            >  5:12223/1‑149 (MQ=255)                                                                                         ttCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTc                                                            <  3:176503/149‑1 (MQ=255)                                                                                            gACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTcgc                                                          >  1:269576/1‑148 (MQ=255)                                                                                            gACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTcgc                                                          <  3:230736/148‑1 (MQ=255)                                                                                               agagCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCgaga                                                      <  6:136985/149‑1 (MQ=255)                                                                                               agagCGCCGGCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCgaga                                                      <  4:71472/149‑1 (MQ=255)                                                                                                gagCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCga                                                        >  3:199936/1‑146 (MQ=255)                                                                                                      cGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTg                                               <  6:12223/149‑1 (MQ=255)                                                                                                            cTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGc                                         >  7:259591/1‑149 (MQ=255)                                                                                                             tAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGTGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGct                                        >  4:180111/1‑149 (MQ=255)                                                                                                               gtgtGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGct                                        <  4:318737/147‑1 (MQ=255)                                                                                                                   gAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCg                                  >  3:90610/1‑149 (MQ=255)                                                                                                                     gggTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCt                                <  3:163180/149‑1 (MQ=255)                                                                                                                      ggTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCtt                               <  5:225316/149‑1 (MQ=255)                                                                                                                      ggTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGGCGCGCGTCTTCGCtt                               >  8:305287/1‑149 (MQ=255)                                                                                                                        tGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCttt                              >  4:6262/1‑148 (MQ=255)                                                                                                                          gCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGaa                           <  5:47091/149‑1 (MQ=255)                                                                                                                          gCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGaa                           <  3:114731/149‑1 (MQ=255)                                                                                                                           cTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGAAg                          <  8:259591/149‑1 (MQ=255)                                                                                                                              aaTCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGAAGGa                        <  4:332552/148‑1 (MQ=255)                                                                                                                                   gCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGAAGGACACGTc                  >  4:40317/1‑149 (MQ=255)                                                                                                                                    cGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGCTTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGAAGGACACGTc                  >  7:265320/1‑148 (MQ=255)                                                                                                                                    cGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGAAGGacac                     <  3:6262/145‑1 (MQ=255)                                                                                                                                      tCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGAAGGACACGTCGAc               <  7:305287/149‑1 (MQ=255)                                                                                                                                           gCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGAAGGACACGTCGACTTTTc          <  4:90610/149‑1 (MQ=255)                                                                                                                                                aGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGAAGGACACGTCGACTTTTCCGGTc     >  7:203204/1‑149 (MQ=255)                                                                                                                                                   cATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGAAGGACACGTCGACTTTTCCGGTCTgg  >  1:354746/1‑149 (MQ=255)                                                                                                                                                   cATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGAAGGACACGTCGACTTTTCCGGTCTgg  <  1:298593/149‑1 (MQ=255)                                                                                                                                                    aTGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGt                                              <  7:289078/104‑1 (MQ=255)                                                                                                                                                    aTGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGt                                              >  8:289078/1‑104 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GAGTGGTTTGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGAAGGACACGTCGACTTTTCCGGTCTGG  >  NC_002947/3388671‑3388965

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