Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,079,884	T→C	6.2%	L360P (CTG→CCG)	rng →	endoribonuclease G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,079,884	0	T	C	6.2%	188.5 / 2.6	64	L360P (CTG→CCG)	rng	endoribonuclease G
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base C (2/2);  total (33/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.12e-01

GAAGAAACCATCTTCAAGACCAACCTCGAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGAGCCGTGCCTGGCCTGTCAGGGCCGCGGCAAG  >  NC_002947/1079742‑1080032                                                                                                                                                |                                                                                                                                                      gaagaaACCATCTTCAAGACCAACCTCGAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGa‑aa                                                                                                                                                 >  8:19813/1‑148 (MQ=255)  aagaaACCATCTTCAAGACCAACCTCGAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGa‑aaaa                                                                                                                                               <  7:314819/149‑1 (MQ=255)         cATCTTCAAGACCAACCTCGAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTg                                                                                                                                         >  3:213421/1‑148 (MQ=255)         cATCTTCAAGACCAACCTCGAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCATGGA‑AAAACAGCTg                                                                                                                                         >  7:306356/1‑148 (MQ=255)         aaTCTTCAAGACCAACATCGAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTgg                                                                                                                                        >  1:184755/2‑149 (MQ=255)          aTCTTCAAGACCAACCTCGAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGa                                                                                                                                       <  5:55465/149‑1 (MQ=255)             ttCAAGACCAACCTCGAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAgcg                                                                                                                                    >  7:296949/1‑149 (MQ=255)              tCAAGACCAACCTCGAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAgcgc                                                                                                                                   >  4:120484/1‑149 (MQ=255)               cAAGACCAACCTCGAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAgcgc                                                                                                                                   >  7:83097/1‑148 (MQ=255)                     cAACCTCGAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACg                                                                                                                            >  4:23109/1‑149 (MQ=255)                     cAACCTCGAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACg                                                                                                                            >  2:118099/1‑149 (MQ=255)                      aaCCTCGAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGc                                                                                                                           <  7:19813/149‑1 (MQ=255)                          tCGAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGcc                                                                                                                          >  7:268107/1‑146 (MQ=255)                            gAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGa                                                                                                                      <  4:60009/148‑1 (MQ=255)                                      ccGCCATCGCCCGACAACTTCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACatcatc                                                                                                           <  3:194323/149‑1 (MQ=255)                                      ccGCCATCGCCCGACAACTGCGGCTGCGCATCATTGGCGGCATCATCATCATCGCCTTAATCGACATGGAAGAAGAAGATCACCAGCGCCAGGTACTGTCGTACCCTGGA‑AATACAGGTGGAGCGCGATCACGCCAAGACCAAcatc                                                                                                              <  1:250497/147‑1 (MQ=255)                                              gCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCCGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCAc                                                                                                   >  7:84681/1‑149 (MQ=255)                                              gCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCCGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCAc                                                                                                   >  7:89327/1‑149 (MQ=255)                                               cccGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCAcc                                                                                                  >  2:223805/1‑149 (MQ=255)                                                 cGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGGCTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGa                                                                                                <  7:113371/149‑1 (MQ=255)                                                      aCTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCCGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTgg                                                                                           <  8:89327/149‑1 (MQ=255)                                                      aCTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCCGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTgg                                                                                           <  8:84681/149‑1 (MQ=255)                                                        tgcgcctgcgcAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGc                                                                                         >  5:23224/1‑149 (MQ=255)                                                           gcctgcgcAACATTGGGGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTg                                                                                      <  8:83097/149‑1 (MQ=255)                                                                cgcAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGc                                                                                  <  3:120484/148‑1 (MQ=255)                                                                  cAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGAACACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGa                                                                               <  4:54039/149‑1 (MQ=255)                                                                     cATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CCTACCCTGGA‑AAAACAGCAGGAGCGCGATCACGCCAAGACCAACATCATCCGCATCAACGAGCTGGGCATGGTGCACATGa                                                                            >  6:42086/1‑149 (MQ=255)                                                                       ttGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGcc                                                                                                                          >  6:212102/1‑101 (MQ=255)                                                                       ttGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGcc                                                                                                                          <  5:212102/101‑1 (MQ=255)                                                                         ggcggcTTCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGGTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTTGGCCTGGTGCAGATGAccc                                                                         <  5:42086/148‑1 (MQ=255)                                                                               atcatcatcatcGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAgcg                                                                  >  7:194395/1‑149 (MQ=255)                                                                                  atcatcatcGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGACCTg                                                                                      <  6:178208/126‑1 (MQ=255)                                                                                  atcatcatcGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGACCTg                                                                                      >  5:178208/1‑126 (MQ=255)                                                                                     atcatcGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCAcgc                                                             <  8:18276/148‑1 (MQ=255)                                                                                        atcGACTTACTCGACATGTAAGAAGAAGCTCACCAGCGCCAGGTACTT‑CGTACCCTGGA‑AAAACAGCTGGATCGCGATCACGCCAAGACCAACATCATCTGCATCCCCGAGCTGGGCCTGGTGCAGATGACCCGAAAGCGCACGCGGGa                                                         <  2:184755/149‑1 (MQ=255)                                                                                           gACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGaaa                                                       <  1:223805/148‑1 (MQ=255)                                                                                            aCTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGc                                                     <  8:268107/149‑1 (MQ=255)                                                                                              ttCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGt                                                                                    <  2:82286/116‑1 (MQ=255)                                                                                              ttCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGt                                                                                    >  1:82286/1‑116 (MQ=255)                                                                                                   cGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGaa                                               >  5:7361/1‑148 (MQ=255)                                                                                                            aGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCttt                                     <  3:23109/149‑1 (MQ=255)                                                                                                             gACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCtt                                      >  1:99308/1‑147 (MQ=255)                                                                                                             gACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCtt                                      <  2:99308/147‑1 (MQ=255)                                                                                                             gACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTg                                    <  8:194395/149‑1 (MQ=255)                                                                                                                gAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGc                                   >  7:77062/1‑147 (MQ=255)                                                                                                                     tCACCAGCGCCAGGTACTT‑CGTACCCTGGA‑ATAACAGCTGGAGCGCGATCACTCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGAGCCGt                            <  8:306356/149‑1 (MQ=255)                                                                                                                     tCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGAGCCGt                            <  6:23224/149‑1 (MQ=255)                                                                                                                     tCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGAGCCGt                            <  4:213421/149‑1 (MQ=255)                                                                                                                        ccAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGAGCCGTgcc                         >  4:320993/1‑149 (MQ=255)                                                                                                                        ccAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGAGCCGTgcc                         >  1:169565/1‑149 (MQ=255)                                                                                                                         cAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGa                                 <  2:212328/140‑1 (MQ=255)                                                                                                                         cAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGa                                 >  1:212328/1‑140 (MQ=255)                                                                                                                            cgcCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCt                                                   <  2:195787/119‑1 (MQ=255)                                                                                                                            cgcCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCt                                                   >  1:195787/1‑119 (MQ=255)                                                                                                                            cgcCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGAGCCGTgcctggc                     <  6:7361/149‑1 (MQ=255)                                                                                                                              ccAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTg                                                                                      >  4:104782/1‑82 (MQ=255)                                                                                                                               cAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGAGCCGTgcctggcctg                  >  1:26536/1‑149 (MQ=255)                                                                                                                                     cTG‑CGTACCCTTCATAAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCGGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGAGCCGTGCCTGGCCTGTCAg              <  1:118099/148‑1 (MQ=255)                                                                                                                                      tG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCg                                                                        >  2:333819/1‑88 (MQ=255)                                                                                                                                      tG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCg                                                                        <  1:333819/88‑1 (MQ=255)                                                                                                                                      tG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGAGCCGTGCCTGGCCTGTCAGGGc           >  4:137227/1‑149 (MQ=255)                                                                                                                                              ccTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGACCAGGTGCTTTGCGAGCCGTGCCTGGCCTGTCAGGGCCGCGGc     >  7:179556/1‑148 (MQ=255)                                                                                                                                              ccTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGAGCCGTGCCTGGCCTGTCAGGGCCGCGGc     >  6:288501/1‑148 (MQ=255)                                                                                                                                              ccTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGAGCCGTGCCTGGCCTGTCAGGGCCGCGGc     >  8:283516/1‑148 (MQ=255)                                                                                                                                                tGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGAGCCGTGCCTGGCCTGTCAGGGCCGCGGCAAg  >  4:256844/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                      GAAGAAACCATCTTCAAGACCAACCTCGAGGCTGCTACCGCCATCGCCCGACAACTGCGCCTGCGCAACATTGGCGGCATCATCATCATCGACTTCATCGACATGGAAGACGAAGATCACCAGCGCCAGGTACTG‑CGTACCCTGGA‑AAAACAGCTGGAGCGCGATCACGCCAAGACCAACATCATCGGCATCACCGAGCTGGGCCTGGTGCAGATGACCCGCAAGCGCACGCGGGAAAGCCTGGAACAGGTGCTTTGCGAGCCGTGCCTGGCCTGTCAGGGCCGCGGCAAG  >  NC_002947/1079742‑1080032

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