Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,023,407	G→T	8.2%	E147* (GAA→TAA) ‡	rmd →	oxidoreductase Rmd

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,023,407	0	G	T	8.2%	187.1 / 6.8	61	E147* (GAA→TAA) ‡	rmd	oxidoreductase Rmd
Reads supporting (aligned to +/- strand):  ref base G (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 = 3.48e-01

GAAGCCATCGCGGCCTGTGGCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGCCGACACCATCGAGCTC  >  NC_002947/2023260‑2023550                                                                                                                                                    |                                                                                                                                                gAAGCCATCGCGGCCTGTGGCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATgg                                                                                                                                                 >  5:62195/1‑148 (MQ=255)      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                                                                                                                            <  4:373298/146‑1 (MQ=255)                        aaCGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTgtggt                                                                                                                            >  3:373298/1‑146 (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                                                                                           <  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                                                                            >  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                                                                                      >  2:218478/1‑135 (MQ=255)                                                                         gTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAAt                                                                                      <  1:218478/135‑1 (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)                                                                                                                                                    |                                                                                                                                                GAAGCCATCGCGGCCTGTGGCAAAACGCCTGAGTGCGTGTTGATCGCCAGCAGTGCCAACGTCTATGGCAACGTTTCCGAAGGGATGCTGGGCGAGCAGACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTGCTGCCCAAGATTGTTTCGCATTTTCGCCGCAAGGCCGACACCATCGAGCTC  >  NC_002947/2023260‑2023550

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