Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,113,568	C→T	6.4%	T71I (ACC→ATC)	PP_3619 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,113,568	0	C	T	6.4%	178.7 / 2.9	62	T71I (ACC→ATC)	PP_3619	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/32);  new base T (2/2);  total (28/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 = 8.97e-01

AACCAAAAATCCACGGACAAATCAGCAAAATAACCACAACATTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAACTCTCTGACATATACCTCAATAGAGCCAC  >  NC_002947/4113421‑4113709                                                                                                                                                    |                                                                                                                                              aaCCAAAAATCCACGGACAAATCAGCAAAATAACCACAACATTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCAcc                                                                                                                                              >  5:298598/1‑149 (MQ=255)  aCCAAAAATCCACGGACAAATCAGCAAAATAACCACAACATTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCAccc                                                                                                                                             <  5:136541/149‑1 (MQ=255)   ccAAAAATCCACGGACAAATCAGCAAAATAACCACAACATTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCa                                                                                                                                            >  7:4311/1‑149 (MQ=255)   ccAAAAATCCACGGACAAATCAGCAAAATAACCACAACATTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACAATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCa                                                                                                                                            <  4:181374/149‑1 (MQ=255)   ccAAAAATCCACGGACAAATCAGCAAAATAACCACAACATTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCAAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCa                                                                                                                                            >  7:151578/1‑149 (MQ=255)    cAAAAATCCACGGACAAATCAGCAAAATAACCACAACATTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCATCCAc                                                                                                                                           >  5:351724/1‑149 (MQ=255)    cAAAAATCCACGGACAAATCAGCAAAATAACCACAACATTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCATCCAc                                                                                                                                           >  5:269048/1‑149 (MQ=255)     aaaaaTCCACGGACAAATCAGCAAAATAACCACAACATTAGCGCAGCTTAGAAATCCATGGGTTTCCGCCTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACt                                                                                                                                          >  5:234865/1‑149 (MQ=255)     aaaaaTCCACGGACAAATCAGCAAAATAACCACAACATTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACt                                                                                                                                          <  8:90810/149‑1 (MQ=255)               ggACAAATCAGCAAAATAACCACAACATTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTa                                                                                                                                <  7:220031/149‑1 (MQ=255)                gACAAATCAGCAAAATAACCACAACATTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTa                                                                                                                                >  7:298161/1‑148 (MQ=255)                     aTCAGCAAAATAACCACAACATTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACa                                                                                                                          <  6:298598/149‑1 (MQ=255)                             aaTAACCACAACATTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAgggg                                                                                                                  <  8:4311/149‑1 (MQ=255)                                     caacaTTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTAc                                                                                                                               >  5:52762/1‑128 (MQ=255)                                     caacaTTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTAc                                                                                                                               <  6:52762/128‑1 (MQ=255)                                      aacaTTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGtat                                                                                                         >  5:253939/1‑149 (MQ=255)                                        caTTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGtattt                                                                                                       >  6:55603/1‑149 (MQ=255)                                             gcgcAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTg                                                                                                   <  2:19545/148‑1 (MQ=255)                                             gcgcAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTg                                                                                                   <  3:19057/148‑1 (MQ=255)                                              cgcAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCa                                                                                                 <  1:3211/149‑1 (MQ=255)                                              cgcAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCAAAGATTTATTCCCCCACGATTTTAGCCATCGCTATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCa                                                                                                 <  4:102660/149‑1 (MQ=255)                                                    ttAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCatat                                                                                                             <  7:241507/131‑1 (MQ=255)                                                    ttAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCatat                                                                                                             >  8:241507/1‑131 (MQ=255)                                                      aGAAATCCTTGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACa                                                                                            <  8:151578/146‑1 (MQ=255)                                                       gAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGCTTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCACTCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGctc                                                                                        <  3:246441/149‑1 (MQ=255)                                                          atcctggggTTCCCGCCTCCGCATTCAGGGCCACCA‑GTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCg                                                                                    <  6:350661/143‑1 (MQ=255)                                                                        gACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATa                                                                                              >  2:76407/1‑126 (MQ=255)                                                                        gACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATa                                                                                              <  1:76407/126‑1 (MQ=255)                                                                        gACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTAtttt                                                                       <  1:170019/149‑1 (MQ=255)                                                                        gACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTAtttt                                                                       >  6:369996/1‑149 (MQ=255)                                                                        gACTCAGCCTGCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTAtttt                                                                       <  4:201549/149‑1 (MQ=255)                                                                         aCTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTAtttt                                                                       <  8:259143/148‑1 (MQ=255)                                                                         aCTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTAtttt                                                                       >  7:259143/1‑148 (MQ=255)                                                                                     tGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTa                                                          <  6:253939/149‑1 (MQ=255)                                                                                      ggCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCt                                                            >  3:261927/1‑146 (MQ=255)                                                                                       gCCAAGATGTATACCGCCACGATTTAAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGCGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAg                                                         <  8:291808/148‑1 (MQ=255)                                                                                        ccAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTTCAATACAGCTCCTCGGACCATTTATTTTTACTCAATTATTAgac                                                       >  4:640/1‑147 (MQ=255)                                                                                        ccAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGcc                                                       >  6:367504/1‑149 (MQ=255)                                                                                         cAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTAAATCAATTACCAACAGGAAGGGGCATATCTAGTTATTTCAATAAAGCTCCTCGGACCATTTAGTTTTACTCGATTCTTAGcc                                                       >  6:414162/1‑148 (MQ=255)                                                                                           aGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTAc                                                    >  5:220106/1‑149 (MQ=255)                                                                                                tataCCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTACTCGGACCATTTATTTTTAATCCATTCTTAGACTACCgtgt                                               >  4:393923/1‑149 (MQ=255)                                                                                                    ccGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTc                                           <  3:640/149‑1 (MQ=255)                                                                                                     ctccACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCATCCACTCAATCAATTACCAACATGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCt                                          <  6:351724/147‑1 (MQ=255)                                                                                                     cGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCATCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCt                                          <  6:269048/149‑1 (MQ=255)                                                                                                       ccACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCtt                                         >  6:834/1‑148 (MQ=255)                                                                                                       ccACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTta                                        <  1:242838/149‑1 (MQ=255)                                                                                                           gATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTtaataa                                    <  7:63443/149‑1 (MQ=255)                                                                                                            aTTTTAGCCCTCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCACTTACCACCAGGAAGGGGCATCTGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATaaa                                   <  6:234865/149‑1 (MQ=255)                                                                                                                    cATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAActct                           <  5:55603/149‑1 (MQ=255)                                                                                                                    cATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAActc                            <  4:261927/148‑1 (MQ=255)                                                                                                                     aTCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAActctc                          >  2:334344/1‑149 (MQ=255)                                                                                                                     aTCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAActctc                          >  3:371003/1‑149 (MQ=255)                                                                                                                        gcgATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGGATATGTATTTATTGCAATACAGCTCATCGGAACATTTATTTTTACTCAATTCTTAGACTACCGTGTTTTCTTAATAAAAAAACTATCTGa                       >  6:218064/1‑149 (MQ=255)                                                                                                                         cgATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAACTCTCTGa                       >  7:318249/1‑148 (MQ=255)                                                                                                                              gTATCCTACACACGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATAGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAACTCTCTGACatat                  <  2:175691/148‑1 (MQ=255)                                                                                                                                aTCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAACTCTCTGACATATAc                <  8:318249/148‑1 (MQ=255)                                                                                                                                 tCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAACTCTCTGACATATACCt              <  5:367504/149‑1 (MQ=255)                                                                                                                                   cTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAACTCTCTGACATATACCTc             >  3:118587/1‑148 (MQ=255)                                                                                                                                       acaGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAACTCTCTGACATATACCTCAATag        <  6:220106/149‑1 (MQ=255)                                                                                                                                         aGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAACTCTCTGACATATACATCAATa         >  5:374983/1‑146 (MQ=255)                                                                                                                                          ggCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAACTCTCTGACATATACCTCAATAGAGc     <  5:369996/149‑1 (MQ=255)                                                                                                                                           gCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGATCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAACTCTCTGACATATACCTCAATAGAGc     >  4:149753/1‑148 (MQ=255)                                                                                                                                             tttttCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAACTCTCTGACATATACCTCAATAGAGcca   >  5:138322/1‑148 (MQ=255)                                                                                                                                               tttCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAACTCTCTGACATATACCTCAATAGAGccac  <  8:298161/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                              AACCAAAAATCCACGGACAAATCAGCAAAATAACCACAACATTAGCGCAGCTTAGAAATCCATGGGTTTCCGACTCAGCCTTCATGGCCAAGATGTATACCGCCACGATTTTAGCCATCGCGATTGTATCCTACACAGGCTTTTTCACCCACTCAATCAATTACCAACAGGAAGGGGCATATGTATTTATTGCAATACAGCTCCTCGGACCATTTATTTTTACTCCATTCTTAGCCTACCGTGTTTTCTTAATAAAAAAACTCTCTGACATATACCTCAATAGAGCCAC  >  NC_002947/4113421‑4113709

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑