Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,023,411	A→C	7.9%	Y148S (TAC→TCC)	rmd →	oxidoreductase Rmd

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,023,411	0	A	C	7.9%	196.7 / 6.2	63	Y148S (TAC→TCC)	rmd	oxidoreductase Rmd
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base C (3/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CATCGCGGCCTGTGGCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGCCGACACCATCGAGCTCGGCAACCT  >  NC_002947/2023265‑2023558                                                                                                                                                   |                                                                                                                                                    cATCGCGGCCTGTGGCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACa                                                                                                                                                   >  8:9489/1‑149 (MQ=255)      cgGCCTGTGGCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGgcg                                                                                                                                              <  6:69480/149‑1 (MQ=255)      cgGCCTGTGGCAAAACGCCTGAGTGCGTGTAGATCGCCAGCAGTGCCAACGTCTATAGCAACGTTTCCGAAGGGATTCTGGGCGAGCAGACACCGCCAGCGCCGGCCAATGCTTACGAGGTCAGCAAGTTGGCCAAGGAATACATGgcg                                                                                                                                              >  8:259727/1‑149 (MQ=255)               gCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGTATTCCATGgc                                                                                                                                               <  4:312924/139‑1 (MQ=255)               gCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGTATTCCATGgc                                                                                                                                               >  3:312924/1‑139 (MQ=255)               gCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCtgtg                                                                                                                                      <  1:2466/148‑1 (MQ=255)               gCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCtgtg                                                                                                                                      <  8:331275/148‑1 (MQ=255)                cAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTgtggt                                                                                                                                    >  4:161988/1‑149 (MQ=255)                cAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTgtggt                                                                                                                                    <  3:161988/149‑1 (MQ=255)                 aaaaCGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTgtggtg                                                                                                                                   >  6:229954/1‑149 (MQ=255)                   aaCGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTgtggt                                                                                                                                    >  3:373298/1‑146 (MQ=255)                   aaCGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTgtggt                                                                                                                                    <  4:373298/146‑1 (MQ=255)                     cGCCTGAGTGCGTGTTGATCGCCAGCAGAGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGAc                                                                                                                               >  1:213675/1‑149 (MQ=255)                           aGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGttt                                                                                                                          <  2:193080/148‑1 (MQ=255)                           aGTGCGGGTTGATCACCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGttt                                                                                                                          <  4:309762/148‑1 (MQ=255)                            gTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGc                                                                                                                        <  8:209124/149‑1 (MQ=255)                                 tgtTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATc                                                                                                                   <  1:275390/149‑1 (MQ=255)                                   ttGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGt                                                                                                                 <  2:213604/149‑1 (MQ=255)                                     gATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGtc                                                                                                                <  7:9489/148‑1 (MQ=255)                                        cGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGtcatca                                                                                                            >  3:372163/1‑149 (MQ=255)                                                   ccAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCttt                                                                                                  <  6:62195/148‑1 (MQ=255)                                                   ccAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCttt                                                                                                  >  8:263004/1‑148 (MQ=255)                                                   ccAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCttt                                                                                                  >  2:169570/1‑148 (MQ=255)                                                    cAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTc                                                                                                 >  3:227490/1‑148 (MQ=255)                                                    cAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCa                                                                                                <  4:227490/149‑1 (MQ=255)                                                     aaCGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCaa                                                                                               >  5:162615/1‑149 (MQ=255)                                                     aaCGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGACCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCaa                                                                                               >  4:248036/1‑149 (MQ=255)                                                     aaCGTCTATGGCAACATTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCtt                                                                                                   <  8:67624/145‑1 (MQ=255)                                                     aaCGTCTATGGCAACATTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCtt                                                                                                   >  7:67624/1‑145 (MQ=255)                                                           tATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTacac                                                                                         >  7:292852/1‑149 (MQ=255)                                                               gCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGgt                                                                                     <  5:229954/149‑1 (MQ=255)                                                                cAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGgtg                                                                                    >  2:357146/1‑149 (MQ=255)                                                                cAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGgtg                                                                                    >  8:145971/1‑149 (MQ=255)                                                                 aaCGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTgtgg                                                                                                                                     <  7:134828/99‑1 (MQ=255)                                                                 aaCGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGAATACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTgtgg                                                                                                                                     >  8:134828/1‑99 (MQ=255)                                                                  aCGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGtt                                                                                  <  7:263004/149‑1 (MQ=255)                                                                    gTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAAt                                                                                              <  1:218478/135‑1 (MQ=255)                                                                    gTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAAt                                                                                              >  2:218478/1‑135 (MQ=255)                                                                         cGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGggcagg                                                                           <  8:292852/149‑1 (MQ=255)                                                                              ggATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCagag                                                                      <  8:74703/149‑1 (MQ=255)                                                                                      ggCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTcc                                                               <  1:169570/148‑1 (MQ=255)                                                                                       gCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAAGACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCctg                                                             >  5:271426/1‑149 (MQ=255)                                                                                         gAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGtt                                                                                  >  3:10214/1‑126 (MQ=255)                                                                                         gAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGtt                                                                                  <  4:10214/126‑1 (MQ=255)                                                                                               aCCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGccc                                                       >  6:241083/1‑147 (MQ=255)                                                                                                         cgcCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCg                                           >  7:205162/1‑149 (MQ=255)                                                                                                             ggCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCAtt                                       >  4:250308/1‑149 (MQ=255)                                                                                                                cAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTc                                    <  6:162615/149‑1 (MQ=255)                                                                                                                cAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGGTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTc                                    <  7:145971/149‑1 (MQ=255)                                                                                                                cAATGATTACGCGGTCAGCAAGTTGCCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTc                                    <  3:248036/149‑1 (MQ=255)                                                                                                                   tGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTcgcc                                 <  8:205162/149‑1 (MQ=255)                                                                                                                            ggTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAg                            <  7:137988/145‑1 (MQ=255)                                                                                                                            ggTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAg                            >  8:137988/1‑145 (MQ=255)                                                                                                                             gTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGCCGa                       <  5:170878/149‑1 (MQ=255)                                                                                                                                aGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGCCGacac                    >  2:309491/1‑149 (MQ=255)                                                                                                                                   agccgcgCCATGTATTCCATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGc                          >  6:207684/7‑140 (MQ=255)                                                                                                                                   agccgcgCCATGTATTCCATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGc                          <  5:207684/134‑1 (MQ=255)                                                                                                                                    aGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGCCGACACCATc                <  1:357146/149‑1 (MQ=255)                                                                                                                                     gTTGGTCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTc                                            >  1:151136/1‑120 (MQ=255)                                                                                                                                     ccgcgCCATGTATTCCATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGCCGACACCATc                >  4:342256/5‑148 (MQ=255)                                                                                                                                          ccATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGCCGACACCATCGAGCTc          >  1:75720/1‑149 (MQ=255)                                                                                                                                                 aTACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGCCGACACCATCGAGCTCGGCAAc    >  7:287431/1‑148 (MQ=255)                                                                                                                                                  tACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATGACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGCCGACACCATCGAGCTCGGCAACCt  <  5:328745/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CATCGCGGCCTGTGGCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGCCGACACCATCGAGCTCGGCAACCT  >  NC_002947/2023265‑2023558

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