Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,017,836	A→T	6.2%	L9* (TTA→TAA)	PP_5461 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,017,836	0	A	T	6.2%	226.6 / 3.8	80	L9* (TTA→TAA)	PP_5461	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/35);  new base T (2/3);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

AGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCAAC  >  NC_002947/1017698‑1017982                                                                                                                                           |                                                                                                                                                    agcaagcaGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACTTTTAGGTGGaca                                                                                                                                           <  3:294369/149‑1 (MQ=255) agcaagcaGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGaca                                                                                                                                           <  5:292161/149‑1 (MQ=255)     agcaGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACTTTTAGGTGGACACCAt                                                                                                                                       <  7:104793/149‑1 (MQ=255)      gcaGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCAtt                                                                                                                                      >  5:302385/1‑149 (MQ=255)          gTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCAGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTg                                                                                                                                  >  1:278171/1‑149 (MQ=255)              cAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGcc                                                                                                                               >  1:186719/1‑148 (MQ=255)               aaaGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCCTTTCTGGCCCa                                                                                                                             >  8:108718/1‑149 (MQ=255)               aaaGCATAGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAAAACATCTTCCTGACGAAACCGGCGAAGGCCATCGTGTCCACCTAAAACTCGACACCAATTCTGGccc                                                                                                                              >  2:299431/1‑148 (MQ=255)                aaGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCAt                                                                                                                            >  7:300301/1‑149 (MQ=255)                 aGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATc                                                                                                                           <  2:331260/149‑1 (MQ=255)                     ttGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCACGc                                                                                                                       >  4:73686/1‑149 (MQ=255)                        aaaaTGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCAACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑a                                                                                                                   >  3:34174/1‑149 (MQ=255)                              ccTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTc                                                                                                              >  4:141430/1‑148 (MQ=255)                              ccTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTc                                                                                                              >  8:104112/1‑148 (MQ=255)                               cTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCa                                                                                                            <  2:117654/149‑1 (MQ=255)                                  gTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAg                                                                                                           <  4:164215/147‑1 (MQ=255)                                   tCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑Attt                                                                                                                >  8:303282/1‑141 (MQ=255)                                   tCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑Attt                                                                                                                <  7:303282/141‑1 (MQ=255)                                      cGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACTTTTAGGTGGac                                                                                                                                            >  4:65049/1‑111 (MQ=255)                                      cGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACTTTTAGGTGGac                                                                                                                                            <  3:65049/111‑1 (MQ=255)                                             aaaaaaaCTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGc                                                                                                         <  6:373834/138‑1 (MQ=255)                                             aaaaaaaCTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGc                                                                                                         >  5:373834/1‑138 (MQ=255)                                             aaaaaaaCTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAgcgc                                                                                              >  5:294047/1‑149 (MQ=255)                                                       acacTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTa                                                                                      <  2:186719/147‑1 (MQ=255)                                                         acTCCCCTGATCGTCGACTCGCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCTGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGT‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAAt                                                                                  <  1:228218/149‑1 (MQ=255)                                                             cccTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAAc                                                                              >  6:338792/1‑149 (MQ=255)                                                              ccTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAg                                                                                                           >  3:147763/1‑119 (MQ=255)                                                              ccTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAg                                                                                                           <  4:147763/119‑1 (MQ=255)                                                               cTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCt                                                                            >  4:133572/1‑149 (MQ=255)                                                                tGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCtt                                                                           >  1:43499/1‑149 (MQ=255)                                                                 gATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCTCCTTAAAGTGGACACCATTTCTGGCCCCGCCAGCGG‑ATTTTCCAGGCATCCTCCGCGCCGTACCTACCATAAACGttt                                                                          >  2:239397/1‑149 (MQ=255)                                                                 gATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCtt                                                                           >  7:287397/1‑148 (MQ=255)                                                                     gTCGAATCCCGGTTTCCGTCCGGCCAACTCATGTTCCAGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGGGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCaaa                                                                      <  8:336950/149‑1 (MQ=255)                                                                     gTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCa                                                                        >  5:365277/1‑147 (MQ=255)                                                                      tCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTAc                                                                                         <  7:170811/129‑1 (MQ=255)                                                                      tCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTAc                                                                                         >  8:170811/1‑129 (MQ=255)                                                                      tCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAg                                                                     >  7:177917/1‑149 (MQ=255)                                                                        gAATCCCGGTATCCGTCCGGCCAACACATCTTCCTTACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAATTTGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGcc                                                                   <  7:113163/149‑1 (MQ=255)                                                                        gAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGcc                                                                   >  4:146539/1‑149 (MQ=255)                                                                        gAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGcc                                                                   >  2:376992/1‑149 (MQ=255)                                                                        gAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATACTCAGCGCCGTACCTACAATAAACCTTTCAAAGc                                                                    >  5:69358/1‑148 (MQ=255)                                                                        gAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGAATCCTCAGCGCCGTACATACAATAAACCTTTCAAAgac                                                                   >  8:274517/1‑147 (MQ=255)                                                                          aTCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCa                                                                 <  2:8039/149‑1 (MQ=255)                                                                             ccGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCaga                                                               <  8:300301/148‑1 (MQ=255)                                                                             ccGGTATCCGTCCGGCAAACACATCTTCCTGACGAAACCGGCTAAGGCCATGGTGTCCACCTAAAAGTGTACACCATTTATGGCCCATAAATCGT‑TTTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCaga                                                               <  1:287391/148‑1 (MQ=255)                                                                                   tccgtccgGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGcgac                                                        >  8:235304/1‑149 (MQ=255)                                                                                            ccAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGgc                                               <  6:365277/149‑1 (MQ=255)                                                                                              aaCACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGgcgc                                             <  1:376992/149‑1 (MQ=255)                                                                                              aaCACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGgcgc                                             <  3:72406/149‑1 (MQ=255)                                                                                               acacATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGgcgc                                             >  1:164198/1‑148 (MQ=255)                                                                                               acacATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCAGCGcc                                            >  2:125917/1‑149 (MQ=255)                                                                                                   aTCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCa                                        <  3:141430/149‑1 (MQ=255)                                                                                                   aTCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCa                                        >  3:308612/1‑149 (MQ=255)                                                                                                     cTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATc                                      <  7:323733/149‑1 (MQ=255)                                                                                                        ccTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGT‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATc                                      <  8:177917/146‑1 (MQ=255)                                                                                                        ccTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGCCGGGCTCGtcgcactctg                                       <  2:278171/145‑7 (MQ=255)                                                                                                        ccTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGc                                    <  7:104112/148‑1 (MQ=255)                                                                                                          tGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCGTTTAAGGCCCATGAAGCGG‑ATTTTACAGGCCTACCCAGCGCCGTTCCTACAATAAACCATTCAAAGCCCAGAGTTCGACGAGCCCGGCGCCTCCATCCCaaa                                 >  1:274126/1‑149 (MQ=255)                                                                                                                  ccGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCt                         >  8:194368/1‑149 (MQ=255)                                                                                                                   cGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCtt                        >  5:144074/1‑149 (MQ=255)                                                                                                                     gCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCTACGAGCCCGGCGCCTCCATCGc                                    >  6:201963/1‑135 (MQ=255)                                                                                                                     gCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGc                                    <  5:201963/135‑1 (MQ=255)                                                                                                                     gCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTgg                      <  3:133572/149‑1 (MQ=255)                                                                                                                     gCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTgg                      <  2:43499/149‑1 (MQ=255)                                                                                                                     gCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTgg                      <  7:274517/149‑1 (MQ=255)                                                                                                                          ggCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCcac                 <  5:338792/149‑1 (MQ=255)                                                                                                                             cATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGGTATATTCCAGGCATCCTCAGCGCCGTACCTACAATAAAACTTTCAAAGCCCAGAGTGAGACGAGAACGGCGACTCCATCGCAAACATCGCCTTGGGCCACAg               >  1:83150/1‑149 (MQ=255)                                                                                                                              aTGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGcc             <  6:302385/149‑1 (MQ=255)                                                                                                                                  tgtCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCaa         >  3:236572/1‑149 (MQ=255)                                                                                                                                   gtCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAAt        >  6:357435/1‑149 (MQ=255)                                                                                                                                    tCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATg       <  8:287397/149‑1 (MQ=255)                                                                                                                                    tCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCAGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAAt        <  1:125917/148‑1 (MQ=255)                                                                                                                                       aCCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCcaca                >  3:128539/1‑137 (MQ=255)                                                                                                                                       aCCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCcaca                <  4:128539/137‑1 (MQ=255)                                                                                                                                        ccTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGgcgc                                             >  3:212265/1‑107 (MQ=255)                                                                                                                                        ccTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGgcgc                                             <  4:212265/107‑1 (MQ=255)                                                                                                                                        ccTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCt                                           >  7:286698/1‑109 (MQ=255)                                                                                                                                        ccTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCaa   <  2:164198/149‑1 (MQ=255)                                                                                                                                        ccTAAAAGTGGACAACATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCt                                           <  8:286698/109‑1 (MQ=255)                                                                                                                                          tAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCAAc  >  5:2848/1‑148 (MQ=255)                                                                                                                                           |                                                                                                                                                    AGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGCCCATCAAGCGG‑ATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCAAC  >  NC_002947/1017698‑1017982

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