Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,505,443	C→G	6.3%	G67R (GGC→CGC)	PP_3103 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,505,443	0	C	G	6.3%	232.8 / 2.1	72	G67R (GGC→CGC)	PP_3103	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base G (2/3);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGGCGATGGCGGACAGCCCCAGGCCCATGAGGATCAGGTT  >  NC_002947/3505301‑3505590                                                                                                                                                |                                                                                                                                                    ctCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCgg                                                                                                                                               <  3:279768/149‑1 (MQ=255)  tCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCTGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCggg                                                                                                                                              >  6:68480/1‑149 (MQ=255)  tCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCggg                                                                                                                                              >  1:232074/1‑149 (MQ=255)   cATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGt                                                                                                                                             >  4:365838/1‑149 (MQ=255)    aTGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTc                                                                                                                                            >  4:84937/1‑149 (MQ=255)           cGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGATGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGACCCGAGCGGCGGTTCCGCGGCGCAGGCCCAggagc                                                                                                                                     <  6:71419/149‑4 (MQ=37)              ccTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACGCGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCg                                                                                                                                  >  2:227532/1‑149 (MQ=255)                  cATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTATCAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGc                                                                                                                               <  1:285621/148‑1 (MQ=255)                   aTCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCtt                                                                                                                             <  7:200416/149‑1 (MQ=255)                         ccccGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGc                                                                                                                        >  6:317623/1‑148 (MQ=255)                            cGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATc                                                                                                                    <  4:255092/149‑1 (MQ=255)                                    gtaggtGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCg                                                                                                             >  4:268222/1‑148 (MQ=255)                                       ggtGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAAc                                                                                                          <  1:114521/148‑1 (MQ=255)                                          gACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCgg                                                                                                               <  6:382276/140‑1 (MQ=255)                                          gACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCgg                                                                                                               >  5:382276/1‑140 (MQ=255)                                          gACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCAGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCCGCGAACCAgg                                                                                                      >  5:272560/1‑149 (MQ=255)                                           aCGATGGCTGTCAGTAACAGCAGCACCACCTGCTCCTGCCCGCGAGCAAGCTGGCGGGCAGCGCGCTCGACACGCTC‑GTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCAATCCTCGGCGAACCAgg                                                                                                      <  2:241893/148‑1 (MQ=255)                                              aTGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGt                                                                                                  <  2:276770/149‑1 (MQ=255)                                                 gCTTTCAGTAACAGCAGCAACAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCGGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTgg                                                                                                <  6:272560/148‑1 (MQ=255)                                                 gCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTgg                                                                                                <  5:326673/148‑1 (MQ=255)                                                  cTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGc                                                                                               >  1:174473/1‑148 (MQ=255)                                                   tGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCa                                                                                             >  2:235176/1‑149 (MQ=255)                                                     tCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATg                                                                                           >  5:389755/1‑149 (MQ=255)                                                        gTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCt                                                                                        >  6:382653/1‑149 (MQ=255)                                                           aCAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACGCGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGc                                                                                     <  1:227532/149‑1 (MQ=255)                                                           aCAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGc                                                                                     >  5:230424/1‑149 (MQ=255)                                                             agcagcACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTg                                                                                   <  5:68480/149‑1 (MQ=255)                                                               cagcaCCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGaa                                                                                 <  2:232074/149‑1 (MQ=255)                                                               cagcaCCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGaa                                                                                 <  1:283347/149‑1 (MQ=255)                                                                agcaCCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAg                                                                                <  8:71324/149‑1 (MQ=255)                                                                 gcaCCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAgg                                                                               <  6:182837/149‑1 (MQ=255)                                                                 gcaCCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAgg                                                                               <  7:370522/149‑1 (MQ=255)                                                                 gcaCCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAgg                                                                               <  1:152048/149‑1 (MQ=255)                                                                 gcaCCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAgg                                                                               <  3:365838/149‑1 (MQ=255)                                                                  caCCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGAGCGATCGACACTCTA‑CTGGGCCTTCGCCGCGGAACCGCCGCTCGGGTCCAGGACAGCGGGCTTGAGGCCATCCACGGAGAAACACGCCGTGGCCATGCCTTGCTCAAggg                                                                              >  8:217556/1‑149 (MQ=255)                                                                    ccAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGACCCGAGCGGCGGTTCCGCGGCGCAGGCCCAg                                                                                                                                         <  8:360739/88‑1 (MQ=18)                                                                    ccAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGACCCGAGCGGCGGTTCCGCGGCGCAGGCCCAg                                                                                                                                         >  7:360739/1‑88 (MQ=18)                                                                    ccAACTGCTCCTGCCAGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTc                                                                            >  3:67205/1‑149 (MQ=255)                                                                       aCTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTCGATGGCCATCCGCCGCGGACCCGCGCGACGGGTCCGGGCCAGCGGGGTTGCGGCCAGCAACGCCCAACACGGCCGTGGCCATGCGTTGGTGAAgggggga                                                                          >  1:91379/1‑145 (MQ=255)                                                                       aCTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGa                                                                          >  3:193355/1‑148 (MQ=255)                                                                        cTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGAca                                                                        <  5:317623/149‑1 (MQ=255)                                                                        cTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGAca                                                                        <  5:382653/149‑1 (MQ=255)                                                                        cTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGAca                                                                        <  3:389353/149‑1 (MQ=255)                                                                         tGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGAcag                                                                       >  6:99996/1‑149 (MQ=255)                                                                            tCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGAcaggca                                                                    <  7:155901/149‑1 (MQ=255)                                                                              cTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGAcaggcagg                                                                  >  6:191973/1‑149 (MQ=255)                                                                                 cccGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAggcaggccg                                                               >  6:49372/1‑149 (MQ=255)                                                                                       aGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAggcaggccggcaggc                                                         <  7:193450/149‑1 (MQ=255)                                                                                         cAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGccc                                                       <  2:174473/149‑1 (MQ=255)                                                                                          aaGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCt                                                      >  6:25418/1‑149 (MQ=255)                                                                                                gCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATaa                                                 <  5:191973/148‑1 (MQ=255)                                                                                                gCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATaa                                                 >  5:106868/1‑148 (MQ=255)                                                                                                        cgcgcTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAg                                                                                <  8:365592/109‑1 (MQ=255)                                                                                                        cgcgcTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAg                                                                                >  7:365592/1‑109 (MQ=255)                                                                                                             tCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATaa                                                 >  4:332798/1‑135 (MQ=255)                                                                                                             tCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATaa                                                 <  3:332798/135‑1 (MQ=255)                                                                                                             tCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGGCGAt                                   >  2:266431/1‑149 (MQ=255)                                                                                                               gACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGGCGAt                                   <  1:266431/147‑1 (MQ=255)                                                                                                                   cGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGGCGATGGCGGa                             <  7:227887/149‑1 (MQ=255)                                                                                                                      tC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGGCGATGGCGGACAg                          >  7:341246/1‑149 (MQ=255)                                                                                                                         cTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGGCGATGGCGGACAGc                         <  3:84937/148‑1 (MQ=255)                                                                                                                         cTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGGCGATGGCGGACAGc                         <  6:106868/148‑1 (MQ=255)                                                                                                                         cTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCAACCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGGCGATGGCGGACAg                          <  7:107579/147‑1 (MQ=255)                                                                                                                         cTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCAACCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGGCGATGGCGGACAg                          >  8:107579/1‑147 (MQ=255)                                                                                                                           gggCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGGCGATGGCGGACAGccc                       >  5:98767/1‑148 (MQ=255)                                                                                                                              ccTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGCCCAAGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGTCGATCTCGTACACccccccg                   >  6:271305/1‑146 (MQ=255)                                                                                                                               cTGCGCCGCTGAACCGGCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATCCTCGGCGATGGCGGACAGCCCCAGGc                  <  7:217556/149‑1 (MQ=255)                                                                                                                                  cgccgcGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGGCGATGGCGGACAGCCCCAGGccc                <  6:389755/148‑1 (MQ=255)                                                                                                                                      gcgGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGGCGATGGCGGACAGCCCCAGGCCCATGAg           >  3:306447/1‑149 (MQ=255)                                                                                                                                       cgGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATaaaa                                               <  2:155171/112‑1 (MQ=255)                                                                                                                                       cgGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATaaaa                                               >  1:155171/1‑112 (MQ=255)                                                                                                                                           aCCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAACATACTCGGCGATGGCGGACAGCCCCAGGCCCATGAGGAt        >  6:347161/1‑147 (MQ=255)                                                                                                                                             cgccgcTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGGCGATGGCGGACAGCCCCAGGCCCATGAGGATCAgg    <  4:67205/149‑1 (MQ=255)                                                                                                                                               ccgcTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGGCGATGGCGGACAGCCCCAGGCCCATGAGGATCAGGtt  >  2:136712/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    CTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCGCTCGACACGCTC‑CTGGGCCTGCGCCGCGGAACCGCCGCTCGGGTCCAGGCCAGCGGGCTTGAGGCCATCCTCGGCGAACCAGGCCGTGGCCATGCCTTGCTGAAGGGTCGACAGGCAGGCCGGCAGGCCCTGATAAAAATACTCGGCGATGGCGGACAGCCCCAGGCCCATGAGGATCAGGTT  >  NC_002947/3505301‑3505590

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