Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	806,091	G→C	6.0%	D413E (GAC→GAG) ‡	PP_0693 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	806,091	0	G	C	6.0%	220.5 / 4.5	67	D413E (GAC→GAG) ‡	PP_0693	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/26);  new base C (2/2);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

CTGGGGCCTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCC  >  NC_002947/805943‑806235                                                                                                                                                     |                                                                                                                                                 cTGGGGCCTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACg                                                                                                                                                  >  6:54115/1‑149 (MQ=255) cTGGGGCCTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACACACAGGTCGCTGTAGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACg                                                                                                                                                  >  7:186994/1‑149 (MQ=255)       ccTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCctgc                                                                                                                                            >  8:253482/1‑149 (MQ=255)          cGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGGCGTTGAACAGACGCATCAGGAGGTCGACGTCCtgctgc                                                                                                                                         >  8:2716/1‑149 (MQ=255)            aaGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCAACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCAGCAGGACGTCGACCGCCTGATGCgt                                                                                                                                       <  4:79181/149‑2 (MQ=255)               cGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCATCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGaa                                                                                                                                    <  8:141537/149‑1 (MQ=255)               cGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGaa                                                                                                                                    <  5:229267/149‑1 (MQ=255)                 ttCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAAct                                                                                                                                  >  8:238484/1‑149 (MQ=255)                      ccaagAAGCCGTCGAGCACCGCATCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGAGAACTCCCAGCGGCAGGCTGAACAGCGAGCCGTAGAACAGACGCATCAGGCGGTCGACGTCCAGCTGCTAGAACTCTTcg                                                                                                                             >  2:331766/5‑149 (MQ=255)                      ccTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTcg                                                                                                                             <  5:54115/149‑1 (MQ=255)                      ccTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTcg                                                                                                                             >  6:351769/1‑149 (MQ=255)                        ttGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTc                                                                                                                              >  4:315103/1‑146 (MQ=255)                        ttGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTc                                                                                                                              <  3:315103/146‑1 (MQ=255)                            aGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAgg                                                                                                                       >  8:151553/1‑149 (MQ=255)                               cGTCGAGCACCGCTTCCAGCGATTTTACGCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCgg                                                                                                                    >  1:207745/1‑149 (MQ=255)                                gTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCggg                                                                                                                   <  8:186994/149‑1 (MQ=255)                                gTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCggg                                                                                                                   >  8:267343/1‑149 (MQ=255)                                    aGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATa                                                                                                               <  5:183711/149‑1 (MQ=255)                                           cTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgttg                                                                                                        >  5:118348/1‑149 (MQ=255)                                              ccAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgttgtt                                                                                                      <  4:276946/148‑1 (MQ=255)                                              ccAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgttgt                                                                                                       <  7:151553/147‑1 (MQ=255)                                               cAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgttgttgt                                                                                                    >  6:282451/1‑149 (MQ=255)                                               cAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgttgttgt                                                                                                    <  6:118348/149‑1 (MQ=255)                                               cAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgttgttgt                                                                                                    >  1:294586/1‑149 (MQ=255)                                                aGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGGtgttgttgt                                                                                                    >  4:22866/1‑148 (MQ=255)                                                 gCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgttgttgt                                                                                                    >  8:320333/1‑147 (MQ=255)                                                          cccACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATcc                                                                                         <  2:294586/149‑1 (MQ=255)                                                                    gCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAgcggggcc                                                                                <  2:207745/148‑7 (MQ=255)                                                                      tGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTg                                                                                             <  8:132624/133‑1 (MQ=255)                                                                      tGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTg                                                                                             >  7:132624/1‑133 (MQ=255)                                                                                cTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGAGGTCCTGCTGCGAGAACGCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCttaa                                                                   <  7:238484/149‑3 (MQ=255)                                                                                   gATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAg                                                                  <  4:268941/147‑1 (MQ=255)                                                                                   gATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAg                                                                  <  5:282451/147‑1 (MQ=255)                                                                                        cgcgATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgt                                                                                                          >  5:102113/1‑102 (MQ=255)                                                                                        cgcgATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgt                                                                                                          <  6:102113/102‑1 (MQ=255)                                                                                         gcgATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGttttt                                                          <  5:351769/149‑1 (MQ=255)                                                                                           gATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGAGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTc                                                         >  2:302911/1‑148 (MQ=255)                                                                                             tctcCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTc                                                         >  7:305951/1‑146 (MQ=255)                                                                                             tctcCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCa                                                       <  8:305951/148‑1 (MQ=255)                                                                                              ctcCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAgg                                                     >  3:67419/1‑149 (MQ=255)                                                                                               tcCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTc                                                         <  8:53121/144‑1 (MQ=255)                                                                                               tcCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTc                                                         >  7:53121/1‑144 (MQ=255)                                                                                                cccAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCGGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTc                                                   >  8:161873/1‑149 (MQ=255)                                                                                                 ccAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCCGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCAAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTcg                                                  >  1:262038/1‑149 (MQ=255)                                                                                                  cAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTcgc                                                 <  7:320333/149‑1 (MQ=255)                                                                                                  cAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTcgc                                                 >  3:175341/1‑149 (MQ=255)                                                                                                                aCAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCGGGTc                                                   >  8:45070/1‑133 (MQ=255)                                                                                                                aCAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCGGGTc                                                   <  7:45070/133‑1 (MQ=255)                                                                                                                 cAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCtcgtc                                    >  3:361123/1‑147 (MQ=255)                                                                                                                           ttGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCgg                         >  5:360494/1‑148 (MQ=255)                                                                                                                             gAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCAc                      >  8:345014/1‑149 (MQ=255)                                                                                                                              aaCAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCAt                                                                                                                <  7:260198/58‑1 (MQ=255)                                                                                                                              aaCAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCAt                                                                                                                >  8:260198/1‑58 (MQ=255)                                                                                                                              aaCAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCAc                      <  7:267343/148‑1 (MQ=255)                                                                                                                                 aGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGaa                  >  5:110565/1‑149 (MQ=255)                                                                                                                                 aGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGaa                  <  6:360494/149‑1 (MQ=255)                                                                                                                                   tcGCAGCAGGACGTCGACCGCCTGATGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTAc                                             <  4:127418/119‑1 (MQ=255)                                                                                                                                   tcGCAGCAGGACGTCGACCGCCTGATGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTAc                                             >  3:127418/2‑120 (MQ=255)                                                                                                                                   tcGCAGCAGGACGTCGACCGCCTGATGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACt                >  7:222999/2‑149 (MQ=255)                                                                                                                                   aCGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACt                >  8:162741/1‑149 (MQ=255)                                                                                                                                   aCGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACt                >  1:185351/1‑149 (MQ=255)                                                                                                                                       aTCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCg            >  3:230009/1‑149 (MQ=255)                                                                                                                                        tCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCCTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGt           >  3:235690/1‑149 (MQ=255)                                                                                                                                            gCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACt       >  4:353486/1‑149 (MQ=255)                                                                                                                                             cGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTg      >  8:319826/1‑149 (MQ=255)                                                                                                                                              ggTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTgg     <  2:69225/149‑1 (MQ=255)                                                                                                                                               gTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGATGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTggg    >  7:324525/1‑149 (MQ=255)                                                                                                                                                tCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGc   <  2:185351/149‑1 (MQ=255)                                                                                                                                                  gACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGcc  >  8:287333/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CTGGGGCCTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCC  >  NC_002947/805943‑806235

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