Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,017,834	C→T	8.3%	G10S (GGT→AGT)	PP_5461 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,017,834	0	C	T	8.3%	203.5 / 6.4	72	G10S (GGT→AGT)	PP_5461	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (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 = 1.03e-01

GCGTCTGCAGTCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGC  >  NC_002947/1017686‑1017979                                                                                                                                                     |                                                                                                                                                    gCGTCTGCAGTCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCAcc                                                                                                                                                     <  1:216684/149‑1 (MQ=255)           tCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGa                                                                                                                                          <  6:153856/149‑1 (MQ=255)           tCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGa                                                                                                                                          >  3:359738/1‑149 (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                                                                                         <  8:37262/149‑1 (MQ=255)                                                           aaaaaCTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCg                                                                                         >  2:5501/1‑149 (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                                                                      <  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                                            >  5:96519/1‑149 (MQ=255)                                                                                                        ccAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGgc                                            <  7:194488/149‑1 (MQ=255)                                                                                                          aaCACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGgcgc                                          >  5:210379/1‑149 (MQ=255)                                                                                                             acaTCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTc                                       >  5:29307/1‑149 (MQ=255)                                                                                                               aTCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCa                                     >  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‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCAt                                    <  5:76019/145‑1 (MQ=255)                                                                                                                     cTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCa                                >  1:91761/1‑148 (MQ=255)                                                                                                                      tGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCaaa                              >  2:214910/1‑149 (MQ=255)                                                                                                                      tGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCaaa                              >  2:120537/1‑149 (MQ=255)                                                                                                                       gACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCaaa                              >  3:161755/1‑148 (MQ=255)                                                                                                                        aCGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAAc                             <  2:399270/148‑1 (MQ=255)                                                                                                                          gAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCTGGCGCCTCCATCGCAAACATc                          <  1:295642/149‑1 (MQ=255)                                                                                                                           aaaCCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCg                         <  4:94995/149‑1 (MQ=255)                                                                                                                               cGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCtt                     >  6:288834/1‑149 (MQ=255)                                                                                                                                 gagAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCATGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACTAGCCCGGCGCCTCCATCGCAAACATCGCCTTgg                   <  7:167675/147‑1 (MQ=255)                                                                                                                                 gCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTgg                   <  2:91761/149‑1 (MQ=255)                                                                                                                                   gAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGc                 <  8:218316/149‑1 (MQ=255)                                                                                                                                      agaaaTGGTGTCCACTTTTAG‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGc                 >  6:369425/5‑146 (MQ=255)                                                                                                                                      agaaaTGGTGTCCACTTTTAG‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGc                 <  5:369425/142‑1 (MQ=255)                                                                                                                                        aaaTGGTGTCCACTTTTAG‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAg            >  6:143771/3‑149 (MQ=255)                                                                                                                                         cATGGTGTCCACCTAAAA‑GTGGACACGATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGc           >  1:218190/1‑149 (MQ=255)                                                                                                                                         cATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTg                    <  7:23732/140‑1 (MQ=255)                                                                                                                                         cATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTg                    >  8:23732/1‑140 (MQ=255)                                                                                                                                           tGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCt         >  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‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAAt     <  4:161755/148‑1 (MQ=255)                                                                                                                                                tCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAAt     <  5:148053/148‑1 (MQ=255)                                                                                                                                                  cACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGc  >  8:182117/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCGTCTGCAGTCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAA‑GTGGACACCATTTCTGGCCCATCAAGCGGATTTTCCAGGCA‑TCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGC  >  NC_002947/1017686‑1017979

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