Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,017,836	A→T	8.3%	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	8.3%	204.2 / 7.5	72	L9* (TTA→TAA) ‡	PP_5461	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/32);  new base T (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 = 8.84e-01

TCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCAAC  >  NC_002947/1017696‑1017982                                                                                                                                             |                                                                                                                                                       tCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGa                                                                                                                                               >  3:359738/1‑149 (MQ=255) tCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGa                                                                                                                                               <  6:153856/149‑1 (MQ=255) tCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGa                                                                                                                                               >  3:365639/1‑149 (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                                                                                              >  2:5501/1‑149 (MQ=255)                                                 aaaaaCTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCg                                                                                              <  8:37262/149‑1 (MQ=255)                                                 aaaaaCTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCg                                                                                              >  6:76019/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                                                                           >  1:399270/1‑148 (MQ=255)                                                                     tAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTc                                                                           <  3:339070/148‑1 (MQ=255)                                                                       gTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAg                                                                                                              <  6:125511/112‑1 (MQ=255)                                                                        tCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAg                                                                       >  1:208544/1‑149 (MQ=255)                                                                        tCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAg                                                                       >  2:374128/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:120537/1‑149 (MQ=255)                                                                                                            tGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCaaa                                 >  2:214910/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                    <  5:369425/142‑1 (MQ=255)                                                                                                                            agaaaTGGTGTCCACTTTTAG‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGc                    >  6:369425/5‑146 (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                       <  7:23732/140‑1 (MQ=255)                                                                                                                               cATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTg                       >  8:23732/1‑140 (MQ=255)                                                                                                                                 tGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGCCACAGCCt            >  6:10039/1‑149 (MQ=255)                                                                                                                                   gtgtCCACTTTTAG‑GTGGACACCATTTCTGGTCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGgcgc                                               <  2:145192/114‑1 (MQ=255)                                                                                                                                   gtgtCCACTTTTAG‑GTGGACACCATTTCTGGTCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGgcgc                                               >  1:145192/1‑114 (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)                                                                                                                                             |                                                                                                                                                       TCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGC‑‑CTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCAAC  >  NC_002947/1017696‑1017982

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