Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,017,839	A→G	8.3%	L8P (CTT→CCT)	PP_5461 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,017,839	0	A	G	8.3%	201.1 / 7.6	72	L8P (CTT→CCT)	PP_5461	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/32);  new base G (3/3);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

TCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCAACCT  >  NC_002947/1017696‑1017984                                                                                                                                                |                                                                                                                                                      tCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGa                                                                                                                                                 >  3:359738/1‑149 (MQ=255) tCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGa                                                                                                                                                 >  3:365639/1‑149 (MQ=255) tCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGa                                                                                                                                                 <  6:153856/149‑1 (MQ=255)        gcaGGTATCAAAGCTTTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACTTTTAG‑GTGGACACCAtg                                                                                                                                          <  7:375096/149‑2 (MQ=255)          aGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTc                                                                                                                                        <  4:365639/149‑1 (MQ=255)          aGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTc                                                                                                                                        <  4:359738/149‑1 (MQ=255)             tATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTgg                                                                                                                                     <  5:228371/149‑1 (MQ=255)                      aTTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAg                                                                                                                            >  4:339070/1‑149 (MQ=255)                      aTTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGACCGGCCAACACATCTTCCTGACGAAACCGCCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAg                                                                                                                            >  6:366804/1‑149 (MQ=255)                       ttGAAAATGCCTTGTCCCGTACACAAAAAAACTTACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTTACGAAACCGGCGAAGGCCATTGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGc                                                                                                                           <  8:286355/149‑1 (MQ=255)                       ttGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTgg                                                                                                                                     >  1:333061/1‑139 (MQ=255)                       ttGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTgg                                                                                                                                     <  2:333061/139‑1 (MQ=255)                       ttGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGc                                                                                                                           >  4:215235/1‑149 (MQ=255)                          aaaaTGCCTTGTCCCGTACACAAAAAAACTGTCACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGa                                                                                                                        >  6:351022/1‑149 (MQ=255)                                  ttGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGATATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAg                                                                                                                >  6:311991/1‑149 (MQ=255)                                            cacaAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAg                                                                                                     <  5:22060/149‑1 (MQ=255)                                               aaaaaaaCTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACAATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAccgc                                                                                                  >  2:105691/1‑149 (MQ=255)                                                 aaaaaCTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCg                                                                                                <  8:37262/149‑1 (MQ=255)                                                 aaaaaCTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCg                                                                                                >  6:76019/1‑149 (MQ=255)                                                 aaaaaCTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCg                                                                                                >  2:5501/1‑149 (MQ=255)                                                    aaCTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGGCGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTa                                                                                              <  1:105691/149‑1 (MQ=255)                                                       tGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTa                                                                                          <  5:366804/149‑1 (MQ=255)                                                         acacTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTa                                                                                          <  5:339585/147‑1 (MQ=255)                                                         acacTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATAGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTAc                                                                                         >  3:369526/1‑148 (MQ=255)                                                          cacTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACaa                                                                                       <  5:69905/149‑1 (MQ=255)                                                           acTCCCCTGATAGTCGAATCCCGATATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAAt                                                                                      <  5:311991/149‑1 (MQ=255)                                                               cccTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATaaa                                                                                   <  2:33941/148‑1 (MQ=255)                                                                ccTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCACTCCTCAGCGCCGTACCTACAATAAAc                                                                                  <  1:5501/149‑1 (MQ=255)                                                                ccTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAAcc                                                                                 >  2:295642/1‑149 (MQ=255)                                                                    aTAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCttt                                                                              <  3:215235/148‑1 (MQ=255)                                                                    aTAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTc                                                                             <  6:191857/149‑1 (MQ=255)                                                                     tAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTc                                                                             <  3:339070/148‑1 (MQ=255)                                                                     tAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTc                                                                             >  1:399270/1‑148 (MQ=255)                                                                       gTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAg                                                                                                                <  6:125511/112‑1 (MQ=255)                                                                        tCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAg                                                                         >  2:374128/1‑149 (MQ=255)                                                                        tCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAg                                                                         >  1:208544/1‑149 (MQ=255)                                                                          gAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGcc                                                                       >  7:218316/1‑149 (MQ=255)                                                                          gAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATAGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGc                                                                        <  4:369526/148‑1 (MQ=255)                                                                                cGGTATCCGACCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGAA‑TCCTCAGCGACGTACCTACAATAAACCTTTCAAAGCCCAGAGt                                                                 >  8:167675/1‑149 (MQ=255)                                                                                        gtccgGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAg                                                         >  1:17223/1‑149 (MQ=255)                                                                                         tccgGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGc                                                        >  8:340247/1‑149 (MQ=255)                                                                                              ccAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGgc                                                   <  7:194488/149‑1 (MQ=255)                                                                                              ccAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGgc                                                   >  5:96519/1‑149 (MQ=255)                                                                                                aaCACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGgcgc                                                 >  5:210379/1‑149 (MQ=255)                                                                                                   acaTCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTc                                            >  5:29307/1‑149 (MQ=255)                                                                                                     aTCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCa                                          >  6:148053/1‑149 (MQ=255)                                                                                                      tCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑ACCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCaga                                                                   >  5:230506/1‑125 (MQ=255)                                                                                                      tCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑ACCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCaga                                                                   <  6:230506/125‑1 (MQ=255)                                                                                                          ccTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCAt                                         <  5:76019/145‑1 (MQ=255)                                                                                                           cTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCa                                     >  1:91761/1‑148 (MQ=255)                                                                                                            tGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCaaa                                   >  2:214910/1‑149 (MQ=255)                                                                                                            tGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCaaa                                   >  2:120537/1‑149 (MQ=255)                                                                                                             gACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCaaa                                   >  3:161755/1‑148 (MQ=255)                                                                                                              aCGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAAc                                  <  2:399270/148‑1 (MQ=255)                                                                                                                gAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCTGGCGC‑‑CTCCATCGCAAACATc                               <  1:295642/149‑1 (MQ=255)                                                                                                                 aaaCCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCg                              <  4:94995/149‑1 (MQ=255)                                                                                                                     cGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCtt                          >  6:288834/1‑149 (MQ=255)                                                                                                                       gagAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCATGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACTAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTgg                        <  7:167675/147‑1 (MQ=255)                                                                                                                       gCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTgg                        <  2:91761/149‑1 (MQ=255)                                                                                                                         gAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGc                      <  8:218316/149‑1 (MQ=255)                                                                                                                            agaaaTGGTGTCCACTTTTAG‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGc                      >  6:369425/5‑146 (MQ=255)                                                                                                                            agaaaTGGTGTCCACTTTTAG‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGc                      <  5:369425/142‑1 (MQ=255)                                                                                                                              aaaTGGTGTCCACTTTTAG‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGCCACAg                 >  6:143771/3‑149 (MQ=255)                                                                                                                               cATGGTGTCCACCTAAAA‑GTGGACACGATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGCCACAGc                >  1:218190/1‑149 (MQ=255)                                                                                                                               cATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTg                         >  8:23732/1‑140 (MQ=255)                                                                                                                               cATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTg                         <  7:23732/140‑1 (MQ=255)                                                                                                                                 tGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGCCACAGCCt              >  6:10039/1‑149 (MQ=255)                                                                                                                                   gtgtCCACTTTTAG‑GTGGACACCATTTCTGGTCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGgcgc                                                 >  1:145192/1‑114 (MQ=255)                                                                                                                                   gtgtCCACTTTTAG‑GTGGACACCATTTCTGGTCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGgcgc                                                 <  2:145192/114‑1 (MQ=255)                                                                                                                                      tCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAAt          <  5:148053/148‑1 (MQ=255)                                                                                                                                      tCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAAt          <  4:161755/148‑1 (MQ=255)                                                                                                                                        cACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGc       >  8:182117/1‑149 (MQ=255)                                                                                                                                             aaaa‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCAGCGGCATCGCAAACATCGCCTTGGGCCACAGCCTCAATCGCAAc    >  1:135610/1‑149 (MQ=255)                                                                                                                                              aaa‑GTGGACACCATTTCTGGCCCCTCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTggg                       >  2:372016/1‑127 (MQ=255)                                                                                                                                              aaa‑GTGGACACCATTTCTGGCCCCTCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTggg                       <  1:372016/127‑1 (MQ=255)                                                                                                                                              aaa‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCAACCt  <  2:218190/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                      TCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCAACCT  >  NC_002947/1017696‑1017984

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