Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,023,409	A→T	8.2%	E147D (GAA→GAT) ‡	rmd →	oxidoreductase Rmd

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,023,409	0	A	T	8.2%	172.0 / 5.8	61	E147D (GAA→GAT) ‡	rmd	oxidoreductase Rmd
Reads supporting (aligned to +/- strand):  ref base A (27/29);  new base T (3/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

CATCGCGGCCTGTGGCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGCCGACACCATCGAGCTCGGCAAC  >  NC_002947/2023265‑2023556                                                                                                                                                 |                                                                                                                                                    cATCGCGGCCTGTGGCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACa                                                                                                                                                 >  8:9489/1‑149 (MQ=255)      cgGCCTGTGGCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGgcg                                                                                                                                            <  6:69480/149‑1 (MQ=255)      cgGCCTGTGGCAAAACGCCTGAGTGCGTGTAGATCGCCAGCAGTGCCAACGTCTATAGCAACGTTTCCGAAGGGATTCTGGGCGAGCAGACACCGCCAGCGCCGGCCAATGCTTACGAGGTCAGCAAGTTGGCCAAGGAATACATGgcg                                                                                                                                            >  8:259727/1‑149 (MQ=255)               gCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGTATTCCATGgc                                                                                                                                             >  3:312924/1‑139 (MQ=255)               gCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGTATTCCATGgc                                                                                                                                             <  4:312924/139‑1 (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                                                                                                >  8:263004/1‑148 (MQ=255)                                                   ccAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCttt                                                                                                >  2:169570/1‑148 (MQ=255)                                                   ccAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCttt                                                                                                <  6:62195/148‑1 (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                                                                                                 >  7:67624/1‑145 (MQ=255)                                                     aaCGTCTATGGCAACATTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCtt                                                                                                 <  8:67624/145‑1 (MQ=255)                                                           tATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTacac                                                                                       >  7:292852/1‑149 (MQ=255)                                                               gCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGgt                                                                                   <  5:229954/149‑1 (MQ=255)                                                                cAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGgtg                                                                                  >  8:145971/1‑149 (MQ=255)                                                                cAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGgtg                                                                                  >  2:357146/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                        <  5:207684/134‑1 (MQ=255)                                                                                                                                   agccgcgCCATGTATTCCATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGc                        >  6:207684/7‑140 (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)                                                                                                                                                 |                                                                                                                                                    CATCGCGGCCTGTGGCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGCCGACACCATCGAGCTCGGCAAC  >  NC_002947/2023265‑2023556

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