Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,388,303	G→A	6.4%	R340R (CGG→CGA)	pyrD →	dihydroorotate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,388,303	0	G	A	6.4%	206.1 / 2.8	72	R340R (CGG→CGA)	pyrD	dihydroorotate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (28/39);  new base A (3/2);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

TCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTC  >  NC_002947/2388155‑2388447                                                                                                                                                     |                                                                                                                                                 tCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGAGGTTGCTATCGTGGCAATGCCGCGGt                                                                                                                                                 <  2:1607338/150‑1 (MQ=255)   aTTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGa                                                                                                                                               <  1:1446765/150‑1 (MQ=255)   aTTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGAAGAGAAGATCGCGGCCGGGGCGAGACTGGTACAGATATACTCGGGATTCATCTACAAAGGGCCGGCACTGATTCGCGAAGCGGTTGATGCTATCGCGGCAATGCCGCGGTGa                                                                                                                                               >  2:391629/1‑150 (MQ=255)    ttGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGAc                                                                                                                                              >  2:1633124/1‑150 (MQ=255)                  aTTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCtgt                                                                                                                                >  1:532168/1‑150 (MQ=255)                   tcnccGAAGGCCGCCATGCTGCCGAGAAGAGCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGNGGTGGACGCNATCGCGGCNATGCCGCGGTGACCAGCCCGCCCCtgtg                                                                                                                               <  1:380999/147‑1 (MQ=255)                    tACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTgg                                                                                                                              <  1:1767406/150‑1 (MQ=255)                     acagAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTTGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTggg                                                                                                                             >  1:619099/4‑150 (MQ=255)                       cGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTggggg                                                                                                                           >  1:1782459/1‑150 (MQ=255)                                  aTGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCgggtgggt                                                                                                                <  2:1566128/150‑1 (MQ=255)                                  aTGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCgggtgggt                                                                                                                <  1:1616161/150‑1 (MQ=255)                                     ctgcCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCa                                                                                                             <  2:930190/150‑1 (MQ=255)                                      tgcCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGGGGGTGCAg                                                                                                            >  1:1739318/1‑150 (MQ=255)                                         cGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGATGACCAGCCCGCGCCTGTGGGCGGGCGGGTGGGTGCAGCtg                                                                                                         >  2:1451342/1‑150 (MQ=255)                                            gaagaTCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCtgatg                                                                                                      >  2:1293749/1‑150 (MQ=255)                                                 tCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGAGGGGGCAGCTGCTGAGGtt                                                                                                 >  2:233027/1‑150 (MQ=255)                                                        cGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAg                                                                                          <  2:368840/150‑1 (MQ=255)                                                         ggggCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAgg                                                                                         <  2:532168/150‑1 (MQ=255)                                                          gggCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGc                                                                                        <  1:1143869/150‑1 (MQ=255)                                                                   gggTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAg                                                                               <  1:46136/149‑1 (MQ=255)                                                                    ggTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAgg                                                                              <  2:1922234/150‑1 (MQ=255)                                                                     gTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAggg                                                                             >  1:598663/1‑150 (MQ=255)                                                                              taCTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCaa                                                                    <  2:456205/150‑1 (MQ=255)                                                                                cTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAgg                                                                  <  1:556408/150‑1 (MQ=255)                                                                                 tCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAggg                                                                 >  1:1265065/1‑150 (MQ=255)                                                                                  cGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCttgacgaag                                                                <  2:935751/150‑8 (MQ=255)                                                                                    ggCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGcc                                                              >  1:1508599/1‑150 (MQ=255)                                                                                      cTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAGGGGccct                                                            <  2:619099/150‑2 (MQ=255)                                                                                       ttCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCAGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCCTAAAAAAGGGCCCATTGAAGGGGCCCCt                                                           >  1:1744393/1‑150 (MQ=255)                                                                                         cATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                         <  2:598663/150‑1 (MQ=255)                                                                                         cATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                         <  1:1084020/150‑1 (MQ=255)                                                                                          aTCTACAAAGGGACGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTggg                                                        >  2:1553282/1‑150 (MQ=255)                                                                                              aCAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCtgat                                                                                                       <  2:1760512/99‑1 (MQ=255)                                                                                              aCAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCtgat                                                                                                       >  1:1760512/1‑99 (MQ=255)                                                                                               cAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTc                                                   <  2:507419/150‑1 (MQ=255)                                                                                                aaaGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcg                                                  <  1:22715/150‑1 (MQ=255)                                                                                                  aGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcgcc                                                <  1:531907/150‑1 (MQ=255)                                                                                                  aGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGAAAGAGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcgcc                                                <  2:1808580/150‑1 (MQ=255)                                                                                                    ggCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAg                                                                   >  2:806358/1‑129 (MQ=255)                                                                                                    ggCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAg                                                                   <  1:806358/129‑1 (MQ=255)                                                                                                       cGGCGCTGATTCGCGAAGCGGTGGATGCTATGGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCACGGGGCCCCTGGGCCGTCGCccgccg                                           >  1:166665/1‑150 (MQ=255)                                                                                                         gcgcTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAg                                                                                                    <  2:132742/91‑1 (MQ=255)                                                                                                         gcgcTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAg                                                                                                    >  1:132742/1‑91 (MQ=255)                                                                                                         gcgcTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTACAGCTGATGAGGTTTGCCAAGTCATAAAAAAGGGCCCCTTCAAGGGGCCCCTTGGCCTTCGCCCgccgcc                                         >  2:214440/1‑150 (MQ=255)                                                                                                             tGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGa                                     >  1:270299/1‑150 (MQ=255)                                                                                                              gATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGAt                                    >  1:1675187/1‑150 (MQ=255)                                                                                                               aTTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGGGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATg                                   <  1:870412/150‑1 (MQ=255)                                                                                                               aTTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAgg                                                                                         <  1:1610247/96‑1 (MQ=255)                                                                                                               aTTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAgg                                                                                         >  2:1610247/1‑96 (MQ=255)                                                                                                               aTTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATg                                   <  1:1293749/150‑1 (MQ=255)                                                                                                                ttCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATgg                                  <  2:1675187/150‑1 (MQ=255)                                                                                                                     gAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCtt                                                    <  1:692087/127‑1 (MQ=255)                                                                                                                     gAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCtt                                                    >  2:692087/1‑127 (MQ=255)                                                                                                                                    accgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa              >  2:985969/3‑150 (MQ=255)                                                                                                                                    aTCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCtga                                                                                                        <  2:611574/60‑1 (MQ=255)                                                                                                                                    aTCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCtga                                                                                                        >  1:611574/1‑60 (MQ=255)                                                                                                                                    aTCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa              <  2:507352/150‑1 (MQ=255)                                                                                                                                     ccgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                         >  2:236152/2‑106 (MQ=255)                                                                                                                                     ccgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                         <  1:236152/105‑1 (MQ=255)                                                                                                                                        cgGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCa          >  1:710710/1‑150 (MQ=255)                                                                                                                                         ggCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAg         <  1:1355595/150‑1 (MQ=255)                                                                                                                                         ggCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGATAACCCAAAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAg         <  1:1553282/150‑1 (MQ=255)                                                                                                                                          gCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAGGGGccctttt                                                         <  2:587466/101‑5 (MQ=255)                                                                                                                                          gCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAGGGGccctttt                                                         >  1:587466/1‑97 (MQ=255)                                                                                                                                          gCAATGCCGCGGTGACCAGCCCGCCCCTGTCGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGTCCCTGGGCCTTCGCCCGCCTCCCGGATGGGGCGTGCATCGCAATTGGTATTCAgt        <  2:1019615/150‑1 (MQ=255)                                                                                                                                           cAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTaaa                   >  2:401449/1‑138 (MQ=255)                                                                                                                                           cAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTaaa                   <  1:401449/138‑1 (MQ=255)                                                                                                                                                ccGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAg                                                                                                            >  1:1418436/1‑44 (MQ=255)                                                                                                                                                ccGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAg                                                                                                            <  2:1418436/44‑1 (MQ=255)                                                                                                                                                ccGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTc  <  1:985969/150‑1 (MQ=255)                                                                                                                                                 cgcgGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGc                           >  1:942249/1‑124 (MQ=255)                                                                                                                                                 cgcgGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGc                           <  2:942249/124‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 TCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTC  >  NC_002947/2388155‑2388447

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑