Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,338,493	A→T	11.3%	intergenic (+113/+29)	PP_3812 → / ← PP_5632	lipase/allo‑isoleucine metabolism protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,338,493	0	A	T	11.3%	116.8 / 3.4	55	intergenic (+113/+29)	PP_3812/PP_5632	lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base T (3/3);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGAA  >  NC_002947/4338347‑4338641                                                                                                                                                    |                                                                                                                                                        ggCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCATTAGCGGGTa‑‑a‑a                                                                                                                                                    >  3:100129/1‑149 (MQ=255)   cGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑Acc                                                                                                                                                  <  2:199904/149‑1 (MQ=255)           gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTT‑‑T‑ACCCGCGAAAg                                                                                                                                          <  6:153362/149‑1 (MQ=255)           gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAg                                                                                                                                          <  3:295261/149‑1 (MQ=255)             aCTTAAAGTGCATCTGCAGTGAGGCCATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGTGCAGGTCCGGCCCTTTCGCGGGTT‑‑T‑ACCCGCGAAAgg                                                                                                                                         <  1:67460/149‑1 (MQ=255)              cTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCATTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTT‑‑T‑ACCCGCGAAAGGGc                                                                                                                                       <  1:306585/149‑1 (MQ=255)                tAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCATTAGCCGGGA‑‑A‑AACCCCGAAAAGGCCg                                                                                                                                     >  5:1833/1‑149 (MQ=255)                 aaaGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCgg                                                                                                                                    >  8:259726/1‑149 (MQ=255)                  aaGTGCATCTGCAGTGAGGCGATCAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATCAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCgg                                                                                                                                    <  3:246799/148‑1 (MQ=255)                               gTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                                                                                                                      <  5:69974/149‑1 (MQ=255)                               gTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                                                                                                                      <  2:269440/149‑1 (MQ=255)                                     atagaaGCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAgg                                                                                                                                         >  4:46121/6‑124 (MQ=38)                                     atagaaGCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAgg                                                                                                                                         <  3:46121/119‑1 (MQ=38)                                           gCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCt                                                                                                          >  2:95748/1‑149 (MQ=255)                                            cAAAAGACCCGCCCAGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCAATAGCGGGTA‑‑A‑ACCCCCGAAAGGGCCCGGAAATACTGCACAGGGCTCAAACCCta                                                                                                         >  4:349535/1‑148 (MQ=255)                                            cAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCtt                                                                                                         >  3:121189/1‑149 (MQ=255)                                             aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                                                                                                        >  5:368304/1‑149 (MQ=255)                                             aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                                                                                                        >  2:50111/1‑149 (MQ=255)                                             aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                                                                                                        <  1:95748/149‑1 (MQ=255)                                              aTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGTTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                                                                                                       <  1:344111/149‑1 (MQ=255)                                                  aCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTAGCGGGTT‑‑A‑ACCCGCCAAACGGCCCGCAAAACCGGCTCAGGGCTCAAACCCATCcaca                                                                                                    >  2:327237/1‑148 (MQ=255)                                                    ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAt                                                                                                  >  1:243463/1‑148 (MQ=255)                                                    ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCATTCGCGGGGG‑‑A‑AACCGCGAAACGGCCGGCACATACAGGTCAGGGCGGAAACCCGGCcaca                                                                                                    >  1:157397/1‑146 (MQ=255)                                                                   ttttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATACTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTc                                                                                  >  3:302985/1‑149 (MQ=255)                                                                    tttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTAGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCg                                                                                 >  5:84817/1‑149 (MQ=255)                                                                         ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCGTTCGCGGGGA‑‑A‑ACCCGCGAAAGGGCCGGGCAATCCTGGGCAGGGCGCAAAACCTTGCACAAAGATCACCACAGCCTTAGCCAcc                                                                            >  3:43655/1‑149 (MQ=255)                                                                              tttGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAGTAGACCCGCGAAAGGGCCGGCAAATCCTGCAAAGGGCGGAAACCCTTCCAAAATGAAAACCTCAGCCTACGCCAcccc                                                                          >  1:303079/1‑149 (MQ=255)                                                                              tttGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCGTTAGCGGGGA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCGCAAACCCGTCCACAATGATCACCACAGCCATAGCCACCCCGTg                                                                       >  2:301288/1‑149 (MQ=255)                                                                                 gCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCg                                                                    >  7:258786/1‑149 (MQ=255)                                                                                   ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGAt                                                                  <  7:259726/149‑1 (MQ=255)                                                                                    gAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATg                                                                 >  5:362439/1‑149 (MQ=255)                                                                                         aGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATgg                                                                <  2:59206/145‑1 (MQ=255)                                                                                         aGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGGA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATgg                                                                >  1:59206/1‑145 (MQ=255)                                                                                          gCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCa                                                           <  1:301288/149‑1 (MQ=255)                                                                                           cccGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCAt                                                          <  6:368304/149‑1 (MQ=255)                                                                                           cccGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCAt                                                          <  6:362439/149‑1 (MQ=255)                                                                                           cccGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCAt                                                          <  6:348096/149‑1 (MQ=255)                                                                                             cGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCATTAGCGGGGT‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCGGGTCAGGGCGCAAACCCGTGCACCAAGGTGACCACCGCCGTCGCCACCCCGTGGCGATGGGTGGATgg                                                        >  8:223085/1‑148 (MQ=255)                                                                                                    ttCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCATTAGCGGGGA‑‑A‑AACCGCGCAAGGGGCGGCAAATACTGCTCAGGGGGCAAACCCGTGCACAAAGATAACCACAGCCGTAGCCACCCCGTGGCGATGGGAGCATGGCTCTTgg                                                 >  6:131349/1‑149 (MQ=255)                                                                                                                 ccAGCGGTGCG‑GCAGGGCCGGCCGTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGcgcg                                        <  7:296458/145‑1 (MQ=255)                                                                                                                  cAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCg                                   <  5:307932/149‑1 (MQ=255)                                                                                                                   aGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCTAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCgg                                  <  1:327237/149‑1 (MQ=255)                                                                                                                   aGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCgg                                  <  1:159872/149‑1 (MQ=255)                                                                                                                    gCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCgg                                  >  1:238894/1‑148 (MQ=255)                                                                                                                        tGCG‑GCAGGTCGGGCCCTTTCCCGGCAA‑‑A‑ACCCGCTAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACa                             <  3:349535/149‑1 (MQ=255)                                                                                                                              cAGGTCCGGCCCTTGCGCGGGTA‑‑A‑ACCCGCTAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGa                        <  2:157397/149‑1 (MQ=255)                                                                                                                                     ggCCCTTTCGCGGGTT‑‑T‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTg                 >  6:157035/1‑149 (MQ=255)                                                                                                                                      gCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTATAGCGCGGAGCGGTAACAGGCGACGGCCTgc                >  4:75313/1‑149 (MQ=255)                                                                                                                                      gCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgc                >  2:38589/1‑149 (MQ=255)                                                                                                                                      gCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgc                <  8:258786/149‑1 (MQ=255)                                                                                                                                        ccTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcgc              <  5:131349/149‑1 (MQ=255)                                                                                                                                        ccTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcgc              >  2:213156/1‑149 (MQ=255)                                                                                                                                          tttCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGcgcg             >  3:66239/1‑148 (MQ=255)                                                                                                                                            tCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa           <  7:158136/148‑1 (MQ=255)                                                                                                                                             cgcgGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTAc          <  2:243463/148‑1 (MQ=255)                                                                                                                                                  gTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATc     >  8:290954/1‑148 (MQ=255)                                                                                                                                                    a‑‑a‑aCCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGaa  <  4:43655/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                        GGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTA‑‑A‑ACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGAA  >  NC_002947/4338347‑4338641

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