Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	74,378	A→T	6.1%	K285M (AAG→ATG) ‡	PP_0063 →	lipid A biosynthesis lauroyl acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	74,378	0	A	T	6.1%	200.9 / 4.1	66	K285M (AAG→ATG) ‡	PP_0063	lipid A biosynthesis lauroyl acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base T (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

AAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTCCTACAGGGACCGCGTCAAGCCGATGTCTTG  >  NC_002947/74242‑74525                                                                                                                                         |                                                                                                                                                      aaGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAGTGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTGCAAGAAGAGACCGGCt                                                                                                                                          >  1:417198/1‑150 (MQ=255) aaGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGACTAGCTCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCt                                                                                                                                          <  1:328908/150‑1 (MQ=255)  aGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTAGAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGACAGTACATGTGGAGCATCAAGCGCTTCAAGAAGCGCCCGGCTg                                                                                                                                         >  2:291099/1‑150 (MQ=255)          cacGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTgg                                                                                                                                        <  1:1129723/141‑1 (MQ=255)          cTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCTTGAAgcgcttc                                                                                                                                             <  2:149260/138‑4 (MQ=255)          cTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCTTGAAgcgcttc                                                                                                                                             >  1:149260/1‑135 (MQ=255)          cTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCTTGAAgcgct                                                                                                                                               >  2:314372/1‑135 (MQ=255)          cTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCTTGAAgcgct                                                                                                                                               <  1:314372/136‑2 (MQ=255)          cTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTgg                                                                                                                                        >  2:1129723/1‑143 (MQ=255)           tCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAACGCTTCAAGAAGCGCCCGGCTGACGAGCcgg                                                                                                                                >  2:504047/1‑149 (MQ=255)             gAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGTCGAGCCgcgc                                                                                                                              >  1:828266/1‑150 (MQ=255)               atgtGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACACTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTg                                                                                                                            <  2:587741/147‑1 (MQ=255)                gCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTgg                                                                                                                           <  2:12300/150‑1 (MQ=255)                      ccGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTaa                                                                                                                     <  2:491082/150‑1 (MQ=255)                      ccGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTaa                                                                                                                     <  2:315377/150‑1 (MQ=255)                      ccGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTaa                                                                                                                     <  2:84687/150‑1 (MQ=255)                         gAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGccc                                                                                                                                              >  2:631386/1‑122 (MQ=255)                         gAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGccc                                                                                                                                              <  1:631386/122‑1 (MQ=255)                               aTGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCAc                                                                                                            >  1:201673/1‑150 (MQ=255)                                 gTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCActg                                                                                                          >  2:526242/1‑148 (MQ=255)                                      gCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCcgc                                                                                                                                <  1:528279/123‑1 (MQ=255)                                      gCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCcgc                                                                                                                                >  2:528279/1‑123 (MQ=255)                                       cACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCctct                                                                                                    >  1:1253354/1‑150 (MQ=255)                                       cACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCctct                                                                                                    >  2:739541/1‑150 (MQ=255)                                               tgAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGTCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑c                                                                                          <  2:828266/150‑1 (MQ=255)                                               tgAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑c                                                                                          <  2:918294/150‑1 (MQ=255)                                                gAATGTGTCGGCGGCGGCGATGAGCAAGGTGGACGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCGGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑Cg                                                                                         >  2:868689/1‑150 (MQ=255)                                                  aTGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTATAGGAA‑‑Cggg                                                                                       >  2:1028463/1‑150 (MQ=255)                                                                 cgATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGCAGACAACGGATGGCGAGCCATATAGCCATTAACTCGTTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGc                                                                        <  1:739541/150‑1 (MQ=255)                                                                 cgATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGc                                                                        >  2:504648/1‑150 (MQ=255)                                                                         aaGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTc                                                                                                              >  1:928328/1‑106 (MQ=255)                                                                         aaGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTc                                                                                                              <  2:928328/106‑1 (MQ=255)                                                                               gtCGAACGCGATGTGGGCGAGTACCCCAGCCAGGACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCAc                                                          <  2:439572/150‑1 (MQ=255)                                                                                  gAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGAtt                                                             <  1:224062/144‑1 (MQ=255)                                                                                  gAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGAtt                                                             >  2:224062/1‑144 (MQ=255)                                                                                            tgCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCaagaa                                                                                                                                                    <  1:408907/49‑1 (MQ=255)                                                                                            tgCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCaagaa                                                                                                                                                    >  2:408907/1‑49 (MQ=255)                                                                                                      ccccAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGa                                                               <  2:803509/122‑1 (MQ=255)                                                                                                      ccccAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGa                                                               >  1:803509/1‑122 (MQ=255)                                                                                                      ccccAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTAGCGGGTAAACCCg                                   >  1:1297171/1‑150 (MQ=255)                                                                                                       cccAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGt                                                                                                         <  2:298211/81‑1 (MQ=255)                                                                                                       cccAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGt                                                                                                         >  1:298211/1‑81 (MQ=255)                                                                                                       cccAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGc                                  <  1:605746/150‑1 (MQ=255)                                                                                                        ccagccagTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTAGTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCt                                 <  2:1261349/150‑1 (MQ=255)                                                                                                         cagccagTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTc                                >  1:889954/1‑150 (MQ=255)                                                                                                         cagccagTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTc                                >  1:962452/1‑150 (MQ=255)                                                                                                           gccagTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTCCt                              >  2:772297/1‑150 (MQ=255)                                                                                                               gTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTCCTACAg                          >  2:499852/1‑150 (MQ=255)                                                                                                                    tgtgGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTCCTACAGGGAcc                     >  2:1419580/1‑150 (MQ=255)                                                                                                                    tgtgGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTCCTACAGGGAcc                     >  2:1189093/1‑150 (MQ=255)                                                                                                                     gtgGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTCCTACAg                          <  2:451035/144‑1 (MQ=255)                                                                                                                     gtgGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTCCTACAg                          >  1:451035/1‑144 (MQ=255)                                                                                                                      tgGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑Cggg                                                                                       <  1:1328861/82‑1 (MQ=255)                                                                                                                      tgGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑Cggg                                                                                       >  2:1328861/1‑82 (MQ=255)                                                                                                                      atgAGGATGAAGCTCTTCAAGAAGCGTGCGGCTGTCGCGCAGGGCTAGTATTAACCCATTCACCGTCTTCTGTACTAA‑‑CGGGTTTACCCGCGAATGCAATTTTTGATGCACCGACGTATTCGCGGGTAACCCCGCTCCTACAGGCACcgc                   <  2:417198/148‑1 (MQ=255)                                                                                                                              gAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCCCCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGCGTAACCCCGCTCCTACAGGGACCGCGTCAAGcc           <  1:504047/150‑1 (MQ=255)                                                                                                                               aaGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTAccc                                                                                >  2:985152/1‑80 (MQ=255)                                                                                                                               aaGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTAccc                                                                                <  1:985152/80‑1 (MQ=255)                                                                                                                               aaGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTCCTACAGGGACCGCGTCAAGCCg          >  2:560652/1‑150 (MQ=255)                                                                                                                                 gcgcTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTCCTACAg                          >  1:906731/1‑132 (MQ=255)                                                                                                                                 gcgcTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTCCTACAg                          <  2:906731/132‑1 (MQ=255)                                                                                                                                  cgcTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTCCTACAGGGACCGCGTCAAGCCGATg       <  1:1189093/150‑1 (MQ=255)                                                                                                                                  cgcTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGAAATTCCCCGTC‑‑CTGTAGGAATTCACTGGTGCCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTCAACCCGCTCCTACAGGGACCGCGTCAAGCCGATg       <  1:526242/150‑1 (MQ=255)                                                                                                                                   gcTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACg                                                     >  1:1185267/1‑103 (MQ=255)                                                                                                                                   gcTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACg                                                     <  2:1185267/103‑1 (MQ=255)                                                                                                                                     ttCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTCTTCAGCCTCCTGCCTGTCGCTATAGGAA‑‑CGGGTTTACACAAGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTCCTACAGGGACCGCGTCAAGCCGATGTCt    <  1:1028463/150‑1 (MQ=255)                                                                                                                                     ttCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACGTACGTATTCGCGGGTAAACCCGCTCCTACAGGGACCGCGTCAAGCCGAGGTCt    >  2:1046445/1‑150 (MQ=255)                                                                                                                                       cAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTACTCGCGGGTAAACCCGCTCCTACAGGGACCGCGTCAAGCCGATGTCTTg  <  2:201673/150‑1 (MQ=255)                                                                                                                                       aaagaagCGCCCGGCTGGCGAGCAGAGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTCCTACAGGGACCGCGTCAAGCCGATGTCTTg  <  2:889954/149‑1 (MQ=255)                                                                                                                                         |                                                                                                                                                      AAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAGCAAGGTGGTCGAACGCTATGTGCGCGAGTACCCCAGCCAGTACATGTGGAGCATGAAGCGCTTCAAGAAGCGCCCGGCTGGCGAGCCGCGCTGGTATTAAGCCATTCACCGTCCTCTGTAGGAA‑‑CGGGTTTACCCGCGAATGCGATGTTGGATTCACCGACGTATTCGCGGGTAAACCCGCTCCTACAGGGACCGCGTCAAGCCGATGTCTTG  >  NC_002947/74242‑74525

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑