Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	479,499	A→T	7.1%	intergenic (‑31/+25)	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,499	0	A	T	7.1%	200.7 / 3.9	70	intergenic (‑31/+25)	cca/ycgB	multifunctional tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase/phosphatase/type IV piliation protein
Reads supporting (aligned to +/- strand):  ref base A (29/36);  new base T (3/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGTGG  >  NC_002947/479351‑479644                                                                                                                                                      |                                                                                                                                                   gATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATaa                                                                                                                                                    >  4:264745/1‑149 (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‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTc                                                                                                                                       <  2:312431/147‑1 (MQ=255)                gCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTc                                                                                                                                       <  6:296671/147‑1 (MQ=255)                gCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCAt                                                                                                                                     >  6:99582/1‑149 (MQ=255)                  aTTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATgg                                                                                                                                   <  4:69715/149‑1 (MQ=255)                  aTTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATgg                                                                                                                                   <  4:292101/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                                                                                                                 <  1:69710/149‑1 (MQ=255)                                    cccACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAggtggt                                                                                                                 <  2:85094/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‑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)                                                                                                                                                  gATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACtc                                   >  2:81159/1‑117 (MQ=255)                                                                                                                                                  gATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACtc                                   <  3:317875/117‑1 (MQ=255)                                                                                                                                                  gATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACtc                                   <  1:81159/117‑1 (MQ=255)                                                                                                                                                  gATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACtc                                   >  4:317875/1‑117 (MQ=255)                                                                                                                                                  gATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGTg   >  6:6233/1‑149 (MQ=255)                                                                                                                                                   aTAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGTgg  <  4:131797/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   GATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCA‑GCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCG‑GCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGTGG  >  NC_002947/479351‑479644

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