Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,812,821	C→G	9.9%	intergenic (+72/‑38)	rplT → / → pheS	50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,812,821	0	C	G	9.9%	179.4 / 11.2	62	intergenic (+72/‑38)	rplT/pheS	50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base G (3/3);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

GTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGCGA  >  NC_002947/2812672‑2812968                                                                                                                                                      |                                                                                                                                                     gTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTATGGCTTAAGTACCCGACCCAGTGACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTc                                                                                                                                                      <  2:731420/150‑1 (MQ=255)   tCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCaa                                                                                                                                                    <  2:1317487/150‑1 (MQ=255)      ggCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGc                                                                                                                                                 >  2:634634/1‑150 (MQ=255)          gATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGACTTCCTTGTTCGCTACAATCTCTCGGCGGTGCTAGCGCGGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCtt                                                                                                                                             <  1:118454/150‑1 (MQ=255)             cTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCTTTGAAGGCTCTTccc                                                                                                                                          <  2:1238574/150‑1 (MQ=255)               ggCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCt                                                                                                                                        >  2:1793768/1‑150 (MQ=255)                         gAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTAt                                                                                                                              <  2:803953/150‑1 (MQ=255)                           aaaaGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCt                                                                                                                            >  1:252977/1‑150 (MQ=255)                               gcggcgTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAg                                                                                                                        >  1:617044/1‑150 (MQ=255)                                cggcgTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGNAGt                                                                                                                       <  2:1187019/150‑1 (MQ=255)                                     tttGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTa                                                                                                                  <  2:617044/150‑1 (MQ=255)                                          tgcGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGa                                                                                                                         <  1:1137236/138‑1 (MQ=255)                                          tgcGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGa                                                                                                                         >  2:1137236/1‑138 (MQ=255)                                             gATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGACAGTCATCGATAGGGGAAGAGCCTTCACAGCTCGTCCCCAATTTTCGTGTCTGGAAGCTGCACATGGcaa                                                                                                          >  1:205270/1‑150 (MQ=255)                                                    gagaAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGa                                                                                                  >  1:333132/1‑150 (MQ=255)                                                    gagaAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGa                                                                                                  >  2:854638/1‑150 (MQ=255)                                                         aGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATgcgc                                                                                             >  1:884529/1‑150 (MQ=255)                                                             aaaGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGt                                                                                         <  2:1079817/150‑1 (MQ=255)                                                               aGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGtct                                                                                       >  1:1808714/1‑150 (MQ=255)                                                                 ccGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTTAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTcc                                                                                     <  2:919651/150‑1 (MQ=255)                                                                  cGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTccc                                                                                    <  1:854638/150‑1 (MQ=255)                                                                         gCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCt                                                                             >  2:1508427/1‑150 (MQ=255)                                                                            tAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTaga                                                                          <  2:1767596/150‑1 (MQ=255)                                                                             aaGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTagag                                                                         <  2:1294993/150‑1 (MQ=255)                                                                               gTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGc                                                                       <  1:634634/150‑1 (MQ=255)                                                                                  cccGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGt                                                                    <  1:1793768/150‑1 (MQ=255)                                                                                   ccGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTcccc                                                                                                                                         <  1:969716/81‑1 (MQ=255)                                                                                   ccGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTcccc                                                                                                                                         >  2:969716/1‑81 (MQ=255)                                                                                     gcgACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGa                                                                 >  1:1248263/1‑150 (MQ=255)                                                                                          aaTCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAAcgcg                                                            >  2:599187/1‑150 (MQ=255)                                                                                          aaTCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAAcgcg                                                            >  2:1377628/1‑150 (MQ=255)                                                                                              aCCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTCCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCCAGGTGGCGCATG‑GTGCTAATGCCTAGAGGCCGTGGAACGCGCTGa                                                        >  1:1866924/1‑150 (MQ=255)                                                                                               cccGGCGGTGCTCGCGCTGCCGGGCATTGAACGTCATCGATAGGGTAAAAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCGAAGCCCTAGAGGCCGTGGAACGCGCTGaa                                                       <  1:108456/150‑1 (MQ=255)                                                                                               cccGGCAGCGCTAGCACCGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGaa                                                       >  1:1574206/1‑150 (MQ=255)                                                                                               cccGGCAGCGCTAGCACCGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGaa                                                       <  2:1574206/150‑1 (MQ=255)                                                                                               cccGGCAGCGCTAGCACCGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGaa                                                       <  2:1574040/150‑1 (MQ=255)                                                                                               cccGGCAGCGCTAGCACCGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGaa                                                       >  1:1574040/1‑150 (MQ=255)                                                                                                           aGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGcc                                           <  1:1377628/150‑1 (MQ=255)                                                                                                            gcgcTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGccc                                          <  1:642893/150‑1 (MQ=255)                                                                                                                 gCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACCTGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGaa                                     <  2:1808714/150‑1 (MQ=255)                                                                                                                   cGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACa                                   <  2:333132/150‑1 (MQ=255)                                                                                                                   cGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACa                                   >  2:1209495/1‑150 (MQ=255)                                                                                                                        aTTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTTTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAAAGCGCTGAAGACATCAACGCCCTGGAACAGAtac                              >  1:418947/1‑148 (MQ=255)                                                                                                                          tGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCgg                            <  2:1248263/150‑1 (MQ=255)                                                                                                                           gAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCggg                           <  1:599187/150‑1 (MQ=255)                                                                                                                           gAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCAACTTTTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCggg                           <  1:1260682/150‑1 (MQ=255)                                                                                                                                tCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAAt                      >  1:686262/1‑150 (MQ=255)                                                                                                                                tCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAAt                      >  2:315526/1‑150 (MQ=255)                                                                                                                                 cATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAAtt                     <  2:252977/150‑1 (MQ=255)                                                                                                                                  aTCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTa                    >  1:1850218/1‑150 (MQ=255)                                                                                                                                   aaaaTAGGGGAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGaa                                                       >  1:1396454/4‑114 (MQ=38)                                                                                                                                   aaaaTAGGGGAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGaa                                                       <  2:1396454/111‑1 (MQ=38)                                                                                                                                      aTAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTAtctc                <  1:315526/150‑1 (MQ=255)                                                                                                                                          gggAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTg                                                                                           >  1:1652791/1‑71 (MQ=255)                                                                                                                                          gggAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTg                                                                                           <  2:1652791/71‑1 (MQ=255)                                                                                                                                           ggAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCg               >  2:761078/1‑146 (MQ=255)                                                                                                                                           ggAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCg               <  1:761078/146‑1 (MQ=255)                                                                                                                                           ggAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCaa           >  2:1904803/1‑150 (MQ=255)                                                                                                                                            gAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCaag          >  2:1566521/1‑150 (MQ=255)                                                                                                                                                agCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTACTTATCTCGGCAAGAAgg      >  1:760676/1‑150 (MQ=255)                                                                                                                                                   cTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGCg   <  1:1511386/150‑1 (MQ=255)                                                                                                                                                    ttCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGCGa  >  2:540667/1‑150 (MQ=255)                                                                                                                                                    ttCAAAGCTCTTCCCCTATTTTCGTATATGGAGGCTGTACATGGGAAG‑CCTGGATGCGCTGGGTTCCCAGGCCCTAGAGGCCGTGGAACGAGCACAAGACAACAACGCCCTGGAATACATCCGTGTTAATCATCTCGGCATGAAGGGCGa  >  2:592296/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGCGA  >  NC_002947/2812672‑2812968

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