Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,143,615	C→A	7.1%	A85E (GCA→GAA)	PP_3647 →	oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,143,615	0	C	A	7.1%	223.8 / 6.1	70	A85E (GCA→GAA)	PP_3647	oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base A (3/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-01

CATTCAGCGCGCTGTCGATCGACCCGGCGCTGGTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGATCATTGTCGGCGCGGTG  >  NC_002947/4143471‑4143763                                                                                                                                                 |                                                                                                                                                     cATTCAGCGCGCTGTCGATCGACCCGGCGCTGGTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAgg                                                                                                                                                   <  6:284330/148‑1 (MQ=255)   ttCAGCGCGCTGTCGATCGACCCGGCGCTGGTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAggcaa                                                                                                                                                <  8:105607/149‑1 (MQ=255)       gcgcgcTGTCGATCGACCCGGCGCTGGTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAggcaaggg                                                                                                                                             >  8:163966/1‑148 (MQ=255)           gcTGTCGATCGACCCGGCGCTGGTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAggcaagggcaag                                                                                                                                         >  6:5169/1‑148 (MQ=255)               tcgatcgaCCCGGCGCTGGTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAg                                                                                                                                                    >  7:137895/1‑133 (MQ=255)               tcgatcgaCCCGGCGCTGGTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAg                                                                                                                                                    <  8:137895/133‑1 (MQ=255)               tcgatcgaCCCGGCGCTGGTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGAc                                                                                                                                      <  6:153693/147‑1 (MQ=255)               tcgatcgaCCCGGCGCTGGTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACaa                                                                                                                                    <  4:270775/149‑1 (MQ=255)                      aCCCGGCGCTGGTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAg                                                                                                                             >  7:8032/1‑149 (MQ=255)                      aCCCGGCGCTGGTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAg                                                                                                                             <  5:189436/149‑1 (MQ=255)                       cccGGCGCTGGTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAgg                                                                                                                            <  5:375588/149‑1 (MQ=255)                         cGGCGCTGGTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCTATCCTGCGCTACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAATGCCAAGGGc                                                                                                                          <  5:191852/149‑1 (MQ=255)                            cgcTGGTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTc                                                                                                                       <  6:179804/149‑1 (MQ=255)                                 gTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGt                                                                                                                   <  7:163966/148‑1 (MQ=255)                                   gCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGa                                                                                                                                       <  1:256251/126‑1 (MQ=255)                                   gCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCAGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGa                                                                                                                                       >  2:256251/1‑126 (MQ=255)                                         cTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTg                                                                                                           <  5:374143/148‑1 (MQ=255)                                            ccccAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTgagc                                                                                                        <  8:8032/148‑1 (MQ=255)                                             cccAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTgagc                                                                                                        >  1:229/1‑147 (MQ=255)                                                      gCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCaa                                                                                             <  7:370360/149‑1 (MQ=255)                                                        cTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAAcc                                                                                           >  4:232783/1‑149 (MQ=255)                                                        cTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAAcc                                                                                           >  3:54314/1‑149 (MQ=255)                                                         tCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAAccc                                                                                          >  8:4838/1‑149 (MQ=255)                                                         tCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAAccc                                                                                          >  4:216677/1‑149 (MQ=255)                                                                  tACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGc                                                                                 >  7:358264/1‑149 (MQ=255)                                                                     ccGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGcc                                                                                >  3:73469/1‑147 (MQ=255)                                                                      cGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAg                                                                             <  6:67840/149‑1 (MQ=255)                                                                       gATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAgg                                                                            >  2:198225/1‑149 (MQ=255)                                                                             gcgcGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACgg                                                                      >  1:187622/1‑149 (MQ=255)                                                                               gcgACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGc                                                                     >  2:49438/1‑148 (MQ=255)                                                                                cgACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGa                                                                   >  4:174904/1‑149 (MQ=255)                                                                                     gcaagcaGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGatcatc                                                              <  8:98745/149‑1 (MQ=255)                                                                                       aagcaGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGatcatc                                                              >  4:297416/1‑147 (MQ=255)                                                                                          caGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTg                                                         <  3:216677/149‑1 (MQ=255)                                                                                             ttCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCg                                                       >  3:162544/1‑148 (MQ=255)                                                                                             ttCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCg                                                       >  2:76132/1‑148 (MQ=255)                                                                                                gCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAg                                                                             >  5:306388/1‑123 (MQ=255)                                                                                                 ccATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCATTCGCAGGCAAGGGAAAGGACAAGGCCAAGGGCGTCGAGTGGCCCCTGAGCGCGCTGGGAACCCCCCTGCTGGCCAAGGCCACGGCGATCAGCGAGGGCGAATtg                                                  >  6:240976/1‑149 (MQ=255)                                                                                                   aTCCCCCTGCTGGCCCCCGTACATCCCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATtgtg                                                <  5:240976/149‑1 (MQ=255)                                                                                                   aTCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATtgtg                                                <  7:4838/149‑1 (MQ=255)                                                                                                   aTCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATtgtg                                                <  1:76132/149‑1 (MQ=255)                                                                                                     ccACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATtgtg                                                >  1:348337/1‑147 (MQ=255)                                                                                                      cACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGgc                                              <  1:49438/148‑1 (MQ=255)                                                                                                       aCCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATtgt                                                 >  3:251645/1‑144 (MQ=255)                                                                                                       aCCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATtgt                                                 <  4:251645/144‑1 (MQ=255)                                                                                                       aCCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGgcgc                                            >  4:324455/1‑149 (MQ=255)                                                                                                         ctgctgTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGa                                          <  6:49639/149‑1 (MQ=255)                                                                                                          tgctgTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCTACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAg                                         >  1:331653/1‑149 (MQ=255)                                                                                                                        cAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGa                           <  3:324455/149‑1 (MQ=255)                                                                                                                         aGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGAc                          <  5:5169/149‑1 (MQ=255)                                                                                                                         aGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGAc                          <  3:174904/149‑1 (MQ=255)                                                                                                                         aGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGAc                          >  2:179101/1‑149 (MQ=255)                                                                                                                                 gAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGAt                  >  8:204228/1‑149 (MQ=255)                                                                                                                                  ttgcccttGCCTGCGAAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGAt                  >  2:107600/8‑148 (MQ=255)                                                                                                                                    gCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCCACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGATCAt               >  2:339999/1‑149 (MQ=255)                                                                                                                                     cccttGCCTGCGAAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTacga                                    >  6:338948/5‑127 (MQ=255)                                                                                                                                     cccttGCCTGCGAAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTacga                                    <  5:338948/123‑1 (MQ=255)                                                                                                                                     cGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGATCAtt              <  3:232783/149‑1 (MQ=255)                                                                                                                                     cGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGATCAtt              <  8:358264/149‑1 (MQ=255)                                                                                                                                      gTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGATCATTg             >  2:394524/1‑149 (MQ=255)                                                                                                                                       cttGCCTGCGAAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGATCATTGt            >  8:386894/3‑149 (MQ=255)                                                                                                                                          gCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGa                   >  3:69002/1‑139 (MQ=255)                                                                                                                                          gCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGa                   <  4:69002/139‑1 (MQ=255)                                                                                                                                          gCCTGCGAAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGATCATTGTCg          <  1:107600/148‑1 (MQ=255)                                                                                                                                            cTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGATCATTGTCGgcg       >  8:84053/1‑149 (MQ=255)                                                                                                                                            cTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGATCATTGTCGgcg       >  2:53508/1‑149 (MQ=255)                                                                                                                                             ttCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGATCATTGTCGgcgc      >  3:99718/1‑149 (MQ=255)                                                                                                                                              tCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGATCATTGTCGGcgcg     >  2:375841/1‑149 (MQ=255)                                                                                                                                               cGCTGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGTGTGGCAGCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGCCCAGGGCGATCATCGTGTGCGATTTGTGGCGCGAGTACGACGGCGGTGACCACGCGATCATTGTCGGcgcg     <  2:187622/148‑1 (MQ=255)                                                                                                                                               cGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGATCATTGTCGGcgcg     <  4:162544/148‑1 (MQ=255)                                                                                                                                                 cAGGCAAGGGCAAGGACAAGGCCAAGGGAGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGATCATTGTCGGCGCGGTg  <  3:331321/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CATTCAGCGCGCTGTCGATCGACCCGGCGCTGGTGCTGTTCTGCCCCAACTACGCCTCCGACACCTACCCGATCCTGCGCGACAGCAAGCAGTTCGCCATCCACCTGCTGTCCGCCGACCAGACCGCCGAGGCGTATGCCTTCGCAGGCAAGGGCAAGGACAAGGCCAAGGGCGTCGAGTGGCACCTGAGCGAGCTGGGCAACCCCCTGCTGGCCAAGGCCACGGCGATCATCGAGTGCGAATTGTGGCGCGAGTACGACGGCGGTGACCACGCGATCATTGTCGGCGCGGTG  >  NC_002947/4143471‑4143763

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