Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	479,486	T→C	6.0%	intergenic (‑18/+38)	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,486	0	T	C	6.0%	199.1 / 2.5	66	intergenic (‑18/+38)	cca/ycgB	multifunctional tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase/phosphatase/type IV piliation protein
Reads supporting (aligned to +/- strand):  ref base T (30/32);  new base C (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CCACCGGCCGATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCAC  >  NC_002947/479342‑479631                                                                                                                                                  |                                                                                                                                                  ccACCGGCCGATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTTCAGCTCTCCATAACTGcc                                                                                                                                                >  5:314878/1‑148 (MQ=255)           aTAGCCATTTGCTAGCATTTATTCCACTCTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAAc                                                                                                                                     <  6:120970/149‑1 (MQ=255)              gCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAAccc                                                                                                                                   <  5:80040/148‑1 (MQ=255)               ccATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATaa                                                                                                                                      >  3:154751/1‑144 (MQ=255)               ccATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATaa                                                                                                                                      <  4:154751/144‑1 (MQ=255)                 aTTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCAAAACTGCCGACAGGATAACCCGTcg                                                                                                                               <  2:6233/149‑1 (MQ=255)                     tcTAGGATTTCTTCGACGGTGGCGCCCTCCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTTTAGATGTGCATGTTCAGCTCTCCATAACTTCCGACGGGTTATCCTGTcggca                                                                                                                            <  2:326105/147‑2 (MQ=255)                         gCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTc                                                                                                                         <  6:13502/147‑1 (MQ=255)                            tttCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATggg                                                                                                                    <  6:264954/149‑1 (MQ=255)                            tttCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATggg                                                                                                                    >  5:82806/1‑149 (MQ=255)                              tcttcGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATgggg                                                                                                                   <  1:347437/148‑1 (MQ=255)                                ttcGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTg                                                                                                                >  4:16300/1‑149 (MQ=255)                                  cGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGtt                                                                                                                 <  6:283193/146‑1 (MQ=255)                                  cGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGc                                                                                                               <  6:256333/148‑1 (MQ=255)                                    aCGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAgg                                                                                                                                          >  6:33369/1‑119 (MQ=255)                                      ggTGGCGCCCACCACCAGCCAGTCGATATCGCTGACTAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCAt                                                                                                             >  1:62802/1‑147 (MQ=255)                                           cgcCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAggtg                                                                                                     <  6:82806/149‑1 (MQ=255)                                           cgcCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAggtg                                                                                                     >  1:188364/1‑149 (MQ=255)                                            gcCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGTTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAggtgg                                                                                                    <  8:259273/149‑1 (MQ=255)                                                    ccagccagTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCag                                                                                                          <  3:205519/135‑1 (MQ=255)                                                    ccagccagTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCag                                                                                                          >  4:205519/1‑135 (MQ=255)                                                        ccagTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAggtggt                                                                                                   >  1:200787/1‑138 (MQ=255)                                                        ccagTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAggtggt                                                                                                   <  2:200787/138‑1 (MQ=255)                                                           gTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACa                                                                                           <  6:171431/143‑1 (MQ=255)                                                           gTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACa                                                                                           >  5:171431/1‑143 (MQ=255)                                                               atatCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGgtggtg                                                                                                  >  1:183958/1‑132 (MQ=255)                                                               atatCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGgtggtg                                                                                                  <  2:183958/132‑1 (MQ=255)                                                               atatCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCa                                                                                 >  4:11467/1‑149 (MQ=255)                                                                        aC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCTGAGGTGGTGGATAACAGCCAGGTCCATGCGGcc                                                                          <  2:104844/147‑1 (MQ=255)                                                                                      gCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATGACCGCCGACGGGTTATCCTGTCGGCAGTTATg                                                                                                                      <  3:165828/90‑1 (MQ=21)                                                                                      gCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATGACCGCCGACGGGTTATCCTGTCGGCAGTTATg                                                                                                                      >  4:165828/1‑90 (MQ=21)                                                                                         aaCGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGc                                                         >  7:166228/1‑148 (MQ=255)                                                                                            ggTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCt                                                        <  5:16531/146‑1 (MQ=255)                                                                                            ggTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCt                                                        >  6:16531/1‑146 (MQ=255)                                                                                             gTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGc                                                    <  2:188364/149‑1 (MQ=255)                                                                                                gcgcACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTATCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTcg                                                 <  4:325078/149‑1 (MQ=255)                                                                                                           cgccAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGcggcgg                                      >  3:213299/1‑149 (MQ=255)                                                                                                           cgccAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGcggcgg                                      >  1:193936/1‑149 (MQ=255)                                                                                                            gccAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGgtggtg                                                                                                  >  3:157034/1‑88 (MQ=255)                                                                                                            gccAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGgtggtg                                                                                                  <  4:157034/88‑1 (MQ=255)                                                                                                             ccAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCa                                    <  3:11467/149‑1 (MQ=255)                                                                                                                    tGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATgtggt                              <  4:329667/148‑1 (MQ=255)                                                                                                                    tGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATgtggt                              <  2:193936/148‑1 (MQ=255)                                                                                                                     gTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATgtggtg                             <  8:166228/148‑1 (MQ=255)                                                                                                                      tAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGa                            <  5:175486/148‑1 (MQ=255)                                                                                                                      tAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGa                            <  2:62802/148‑1 (MQ=255)                                                                                                                       aGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCcgcgc                                           >  7:205360/1‑132 (MQ=255)                                                                                                                       aGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCcgcgc                                           <  8:205360/132‑1 (MQ=255)                                                                                                                       aGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGATGAg                           >  3:292770/1‑148 (MQ=255)                                                                                                                       aGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCCCTATGCTCGCCGCGCGGCGCCATGTGGTGAg                           >  7:108914/1‑148 (MQ=255)                                                                                                                           gtgcatgtgcaGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACt                       >  1:80315/1‑148 (MQ=255)                                                                                                                            tgcatgtgcaGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACtca                     <  4:141555/149‑1 (MQ=255)                                                                                                                             gcatgtgcaGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACtcat                    >  5:332350/1‑149 (MQ=255)                                                                                                                               atgtgcaGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACtcatca                  >  7:18778/1‑149 (MQ=255)                                                                                                                                    caGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGATGAGTACTCACCACATGc             <  4:292770/149‑1 (MQ=255)                                                                                                                                     aGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGa                            <  4:8504/133‑1 (MQ=255)                                                                                                                                     aGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGa                            >  3:8504/1‑133 (MQ=255)                                                                                                                                       aaccccATGACCGCCGACGGGTTATCCTGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCt                                                   >  1:80634/5‑108 (MQ=37)                                                                                                                                       aaccccATGACCGCCGACGGGTTATCCTGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGGCGGCGCTTTCGCTATGCt                                                   <  2:80634/104‑1 (MQ=25)                                                                                                                                          tcCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCt       >  6:342187/1‑149 (MQ=255)                                                                                                                                          tcCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCt       >  7:242777/1‑149 (MQ=255)                                                                                                                                           ccATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCtt      >  1:233580/1‑149 (MQ=255)                                                                                                                                            cATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGa                            <  6:103152/126‑1 (MQ=255)                                                                                                                                            cATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGa                            >  5:103152/1‑126 (MQ=255)                                                                                                                                            cATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCACCTTg     >  7:161519/1‑149 (MQ=255)                                                                                                                                               aaCTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCAc  >  2:363602/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CCACCGGCCGATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGAC‑AGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCAC  >  NC_002947/479342‑479631

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