Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,643,261	G→A	8.0%	intergenic (+46/‑77)	dbpA → / → yhiN	23S rRNA specific ATP‑dependent RNA helicase/FAD/NAD(P)‑binding domain‑containing oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,643,261	0	G	A	8.0%	180.5 / 3.8	63	intergenic (+46/‑77)	dbpA/yhiN	23S rRNA specific ATP‑dependent RNA helicase/FAD/NAD(P)‑binding domain‑containing oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base A (3/2);  total (32/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 = 8.84e-01

CAGGCGTTTGTGGCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCT  >  NC_002947/5643114‑5643399                                                                                                                                                     |                                                                                                                                            cAGGCGTTTGTGCTCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCgg                                                                                                                                             <  1:120501/148‑1 (MQ=255)          gtgGCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCATGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAAt                                                                                                                                     >  1:69210/1‑148 (MQ=255)             gCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCAcc                                                                                                                                <  2:219851/149‑1 (MQ=255)              ccGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCgg                                                                                                                                             >  2:346648/1‑135 (MQ=255)              ccGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCgg                                                                                                                                             <  1:346648/135‑1 (MQ=255)                          cgcTGGCCAAGCAGGCCATGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCgg                                                                                                                   <  6:179750/149‑1 (MQ=255)                          cgcTGGCCAAGCAGGCCATGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCgg                                                                                                                   <  2:108383/149‑1 (MQ=255)                            cTGGCCAAGCAGGCCATGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGt                                                                                                                 <  7:26477/149‑1 (MQ=255)                                 cAAGCAGGCCATGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCa                                                                                                                                  >  4:267613/1‑127 (MQ=255)                                 cAAGCAGGCCATGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCa                                                                                                                                  <  3:267613/127‑1 (MQ=255)                                 cAAGCAGGCCATGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAAcc                                                                                                            >  6:78107/1‑149 (MQ=255)                                 cAAGCAGGCCATGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAAcc                                                                                                            >  3:246180/1‑149 (MQ=255)                                    gCAGGCCATGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGc                                                                                                         <  8:231042/149‑1 (MQ=255)                                      aGGCCATGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATg                                                                                                                                    >  6:108597/1‑120 (MQ=255)                                            tGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAg                                                                                                 >  6:319430/1‑149 (MQ=255)                                               aGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGt                                                                                              >  7:259958/1‑149 (MQ=255)                                                      tGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACACGCGAATGCGGCAGTGAATTCACCATCGCATTCGCGGGTGAACACTCTCCCACAGGGTTctcca                                                                                        >  8:297665/1‑145 (MQ=38)                                                       gAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCAt                                                                                       >  4:306615/1‑149 (MQ=255)                                                         aCAGCGGCAAGATCATGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCAtt                                                                                      >  2:115513/1‑148 (MQ=255)                                                         aCAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCa                                                                                           <  8:105732/143‑1 (MQ=255)                                                         aCAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCa                                                                                           >  7:105732/1‑143 (MQ=255)                                                         aCAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCAttt                                                                                     <  5:78107/149‑1 (MQ=255)                                                          cAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTc                                                                                    <  8:259958/149‑1 (MQ=255)                                                              ggCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑t                                                                               <  2:69210/149‑1 (MQ=255)                                                               gCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑Tc                                                                              >  5:364270/1‑149 (MQ=255)                                                                   gATCAAGGGTCGCGCCCTGAAAGTACTCATTATCTATAGCCAGGCAAAGCCCCCTGTATTAGCGGGTTTACCCGCTAATGCTTCAGTGAATTCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGaa                                                                          <  7:297665/149‑1 (MQ=37)                                                                    aTCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAggg                                                                                               >  1:63724/1‑128 (MQ=255)                                                                    aTCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAggg                                                                                               <  2:63724/128‑1 (MQ=255)                                                                     tCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCaaa                                                                                 <  8:100238/141‑1 (MQ=255)                                                                     tCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCaaa                                                                                 >  7:100238/1‑141 (MQ=255)                                                                     tCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAgg                                                                        <  5:319430/149‑1 (MQ=255)                                                                       aaGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGtt                                                                      <  4:246180/149‑1 (MQ=255)                                                                              gcgcCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCaaa                                                                                 >  6:280698/1‑132 (MQ=255)                                                                              gcgcCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCaaa                                                                                 <  5:280698/132‑1 (MQ=255)                                                                                  ccTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCa                                                           >  4:234452/1‑149 (MQ=255)                                                                                    tGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACt                                                         >  7:268689/1‑149 (MQ=255)                                                                                    tGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACt                                                         >  1:314406/1‑149 (MQ=255)                                                                                       aaGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCAc                                                          >  8:214613/1‑145 (MQ=255)                                                                                       aaGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCAc                                                          <  7:214613/145‑1 (MQ=255)                                                                                            cGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGAc                                                 <  3:306615/149‑1 (MQ=255)                                                                                              cATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGt                                               >  3:189802/1‑149 (MQ=255)                                                                                                     tAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGatcatc                                        >  1:264107/1‑149 (MQ=255)                                                                                                         gCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTg                                    <  8:268689/149‑1 (MQ=255)                                                                                                         gCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTg                                    <  2:264107/149‑1 (MQ=255)                                                                                                          ccAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGaa                                                                          >  4:117546/1‑110 (MQ=255)                                                                                                          ccAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGaa                                                                          <  3:117546/110‑1 (MQ=255)                                                                                                            aGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCaaa                                                                                 >  4:325747/1‑102 (MQ=255)                                                                                                            aGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCaaa                                                                                 <  3:325747/102‑1 (MQ=255)                                                                                                              gCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcgcc                               <  1:115513/149‑1 (MQ=255)                                                                                                                    cccccTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcgccggcgcc                         <  6:364270/149‑1 (MQ=255)                                                                                                                    cccccTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcgccggcgcc                         >  3:124207/1‑149 (MQ=255)                                                                                                                      cccTGTGGGAGCGGGTTCACCCGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGcgccgc                       >  1:9380/1‑149 (MQ=18)                                                                                                                       ccTGTGGGAGCGGGTTCACCCGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAACTCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTg                                    <  1:72824/136‑1 (MQ=14)                                                                                                                       ccTGTGGGAGCGGGTTCACCCGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTg                                    >  2:72824/1‑135 (MQ=17)                                                                                                                       ccTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCgccgcc                      >  2:296790/1‑149 (MQ=255)                                                                                                                         tGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACtt                                                        >  2:91719/1‑112 (MQ=255)                                                                                                                             ggcgCGGGTTAACCCGCTAATGCGGAAGTGAATGCACCATCGCTTTCGCGGGTTAACCGGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCt                 <  2:314406/145‑1 (MQ=255)                                                                                                                              gAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGa               >  7:30824/1‑149 (MQ=255)                                                                                                                               aGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGAt              <  3:234452/149‑1 (MQ=255)                                                                                                                                 cGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGAtgt            <  4:189802/149‑1 (MQ=255)                                                                                                                                 cGGGTTCACCCGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGAtgt            <  2:9380/149‑1 (MQ=21)                                                                                                                                   ggTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTgc          <  1:296790/149‑1 (MQ=255)                                                                                                                                          ccGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTc                                                                                    <  5:271300/69‑1 (MQ=255)                                                                                                                                          ccGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTc                                                                                    >  6:271300/1‑69 (MQ=255)                                                                                                                                           cgcgAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt  >  5:2797/1‑149 (MQ=25)                                                                                                                                            gcgAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt  >  5:130527/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                            CAGGCGTTTGTGGCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGC‑TGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAA‑TCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCT  >  NC_002947/5643114‑5643399

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