Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	479,493	A→G	7.0%	intergenic (‑25/+31)	cca ← / ← ycgB	multifunctional tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase/phosphatase/type IV piliation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	479,493	0	A	G	7.0%	206.8 / 3.5	71	intergenic (‑25/+31)	cca/ycgB	multifunctional tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase/phosphatase/type IV piliation protein
Reads supporting (aligned to +/- strand):  ref base A (30/36);  new base G (3/2);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CCGGCCGATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTC  >  NC_002947/479345‑479635                                                                                                                                                      |                                                                                                                                                ccGGCCGATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACa                                                                                                                                                 <  3:195145/149‑1 (MQ=255)       gATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATaa                                                                                                                                           >  4:264745/1‑149 (MQ=255)            ccATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAg                                                                                                                                                <  1:140079/139‑1 (MQ=255)            ccATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAg                                                                                                                                                >  2:140079/1‑139 (MQ=255)            ccATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGt                                                                                                                                      >  7:303605/1‑149 (MQ=255)            ccATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGt                                                                                                                                      >  1:312431/1‑149 (MQ=255)                ttGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTcggc                                                                                                                                  >  1:244227/1‑149 (MQ=255)                  gCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTcggcg                                                                                                                                 <  7:73383/148‑1 (MQ=255)                   cTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCACGCTCTCCATAACTGCCGACAGGATAACCCGTcggcg                                                                                                                                 <  1:241370/148‑1 (MQ=255)                    tAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGt                                                                                                                               >  1:171905/1‑148 (MQ=255)                      gCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGa                                                                                                                                              >  1:175780/1‑131 (MQ=255)                      gCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGa                                                                                                                                              <  2:175780/131‑1 (MQ=255)                      gCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTc                                                                                                                              <  6:296671/147‑1 (MQ=255)                      gCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTc                                                                                                                              <  2:312431/147‑1 (MQ=255)                      gCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCAt                                                                                                                            >  6:99582/1‑149 (MQ=255)                        aTTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATgg                                                                                                                          <  4:292101/149‑1 (MQ=255)                        aTTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATgg                                                                                                                          <  4:69715/149‑1 (MQ=255)                         tttCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACGGGTTATCCCGTCGGCAGTTATg                                                                                                                           >  2:352915/1‑147 (MQ=255)                         tttCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATggg                                                                                                                         <  3:186082/149‑1 (MQ=255)                          ttcttcGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACGGGTTATCCCGTCGGCAGTTATgga                                                                                                                         <  1:352915/148‑2 (MQ=255)                           tcttcGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGt                                                                                                                       >  5:44602/1‑149 (MQ=255)                            cttcGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGtt                                                                                                                      >  8:176155/1‑149 (MQ=255)                            cttcGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGt                                                                                                                       <  3:273827/148‑1 (MQ=255)                             ttcGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTc                                                                                                                                     >  3:39370/1‑133 (MQ=255)                             ttcGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTc                                                                                                                                     <  4:39370/133‑1 (MQ=255)                             ttcGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTg                                                                                                                     >  6:57358/1‑149 (MQ=255)                             ttcGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTg                                                                                                                     >  7:6008/1‑149 (MQ=255)                                     tGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCagag                                                                                                             >  1:85094/1‑149 (MQ=255)                                         gcCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAggtgg                                                                                                         >  2:69710/1‑149 (MQ=255)                                          cccACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAggtggt                                                                                                        <  2:85094/149‑1 (MQ=255)                                          cccACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAggtggt                                                                                                        <  1:69710/149‑1 (MQ=255)                                               caccaGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCaga                                                                                                              <  2:117743/138‑1 (MQ=255)                                               caccaGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCaga                                                                                                              >  1:117743/1‑138 (MQ=255)                                               caccaGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGAt                                                                                                    >  6:176808/1‑148 (MQ=255)                                                  cagccagTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGTTTGCATCAGAGGTGGTGGATaa                                                                                                  <  8:303605/147‑1 (MQ=255)                                                  cagccagTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGtt                                                                                                                      <  8:342558/127‑1 (MQ=255)                                                  cagccagTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACa                                                                                                <  5:181197/149‑1 (MQ=255)                                                     ccagTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAggt                                                                                                           >  2:165737/1‑135 (MQ=255)                                                     ccagTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAggt                                                                                                           <  1:165737/135‑1 (MQ=255)                                                        gTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGa                                                                                                                                              >  2:140466/1‑97 (MQ=255)                                                        gTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGa                                                                                                                                              <  1:140466/97‑1 (MQ=255)                                                        gTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAgg                                                                                          >  3:250119/1‑149 (MQ=255)                                                         tCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGt                                                                                         <  2:244227/149‑1 (MQ=255)                                                           gATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTc                                                                                        <  4:250119/148‑1 (MQ=255)                                                                cGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCTTGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCg                                                                                  >  7:324782/1‑149 (MQ=255)                                                                cGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGc                                                                                                                    >  2:336384/1‑115 (MQ=255)                                                                cGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGc                                                                                                                    <  1:336384/115‑1 (MQ=255)                                                                      cAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATa                                                                            <  5:57358/149‑1 (MQ=255)                                                                          gcggcCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCg                                                                        >  3:131797/1‑149 (MQ=255)                                                                           cggcCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCgg                                                                       <  3:264745/149‑1 (MQ=255)                                                                                gAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCtt                                                                  <  3:334440/149‑1 (MQ=255)                                                                                   cAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCg                                                               >  2:134680/1‑149 (MQ=255)                                                                                             cgcACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGc                                                                                <  2:238770/122‑1 (MQ=255)                                                                                                      ccgccAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCTGCGCTTTCGCTATGCTCGCCGCgcggcg                                            <  1:203931/149‑1 (MQ=255)                                                                                                                tGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATgtggt                                  <  6:44602/148‑1 (MQ=255)                                                                                                                tGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATgtgg                                   <  7:176155/147‑1 (MQ=255)                                                                                                                  tAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGa                                <  2:171905/148‑1 (MQ=255)                                                                                                                          catgtgca‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACtcat                        <  5:99582/148‑1 (MQ=255)                                                                                                                             gtgca‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACt                    >  6:287027/1‑149 (MQ=255)                                                                                                                             gtgca‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACt                    >  3:292983/1‑149 (MQ=255)                                                                                                                               gca‑GCTCTCCATAACTGCCGACAGGATAACCAGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACTTg                  >  5:230808/1‑149 (MQ=255)                                                                                                                                ca‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACTTGc                 >  1:319554/1‑149 (MQ=255)                                                                                                                                    ctctCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACTTGCTCg              <  5:176808/149‑1 (MQ=255)                                                                                                                                        ccATGACCGCCGACGGGTTATCCTGTCGGCAGTTATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACTTGCTCGCCtt          >  5:240/1‑149 (MQ=38)                                                                                                                                        ccATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGGTATGCTCGCCGCGCGGCG‑GCATGTGGTGAg                               <  1:29968/128‑1 (MQ=255)                                                                                                                                        ccATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAg                               >  2:29968/1‑128 (MQ=255)                                                                                                                                         cATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGAGCAAGTGATGAGTACTCACCACATGccgcc             <  6:134609/146‑5 (MQ=39)                                                                                                                                         cATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACTTGCTCGCCTTg         <  1:134680/149‑1 (MQ=255)                                                                                                                                             aCTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCAc      <  5:287027/148‑1 (MQ=255)                                                                                                                                             aCTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACg     >  3:74460/1‑149 (MQ=255)                                                                                                                                                 ccGACGGGTTATCCTGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGtc  <  6:19121/148‑1 (MQ=255)                                                                                                                                                 ccGACGGGTTATCCTGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGtc  >  5:19121/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                CCGGCCGATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTC  >  NC_002947/479345‑479635

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