Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,388,295	A→T	7.1%	M338L (ATG→TTG)	pyrD →	dihydroorotate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,388,295	0	A	T	7.1%	210.7 / 4.6	68	M338L (ATG→TTG)	pyrD	dihydroorotate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (25/38);  new base T (3/2);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

TTGCCGATCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTC  >  NC_002947/2388148‑2388444                                                                                                                                                    |                                                                                                                                                      ttGCCGATCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATg                                                                                                                                                     <  2:573647/150‑1 (MQ=255) ttGCCGATCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATg                                                                                                                                                     <  1:296133/150‑1 (MQ=255)   gCCGATCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGcc                                                                                                                                                   >  2:146608/1‑150 (MQ=255)   gCCGATCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGcc                                                                                                                                                   >  2:1258543/1‑150 (MQ=255)     cGATCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCcgc                                                                                                                                                 >  1:841988/1‑150 (MQ=255)     cGATCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTAATTCGCGAAGCAGTGGATGCTATCGCGGCAATGCcgc                                                                                                                                                 >  1:407094/1‑150 (MQ=255)      gATCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCcgcg                                                                                                                                                <  2:835198/150‑1 (MQ=255)      gATCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCcgcg                                                                                                                                                >  2:293527/1‑150 (MQ=255)                  cgGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAgcccg                                                                                                                                    <  1:1338281/150‑1 (MQ=255)                         aTTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCggt                                                                                                                             >  1:1008757/1‑150 (MQ=255)                           tACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTgg                                                                                                                           <  2:159486/150‑1 (MQ=255)                             ccGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGcccccccacag                                                                                                                         >  2:181401/1‑143 (MQ=255)                             ccGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTgggg                                                                                                                         <  2:1091778/150‑1 (MQ=255)                              cGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACGAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTggggg                                                                                                                        >  2:983404/1‑150 (MQ=255)                              cGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTggggg                                                                                                                        >  2:1211231/1‑150 (MQ=255)                              cGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTggggg                                                                                                                        >  1:245930/1‑150 (MQ=255)                              cGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTggggg                                                                                                                        <  1:779567/150‑1 (MQ=255)                              cGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTACCGCGGCATTGCcgcgcta                                                                                                                                             <  1:1421588/129‑4 (MQ=255)                              cGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTACCGCGGCATTGCcgcgata                                                                                                                                             >  2:1421588/1‑126 (MQ=255)                                 aGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGc                                                                                                                     <  2:841988/150‑1 (MQ=255)                                   gccgccATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGTTGGATGCTGTCGCGGCATTTCTGCGATGACCAGCCCGCCCCTGTGGGGGGGCgg                                                                                                                   >  2:1007733/1‑150 (MQ=255)                                     cgccATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCgggt                                                                                                                 <  2:452985/150‑1 (MQ=255)                                       ccATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCgggtgg                                                                                                               <  2:245930/150‑1 (MQ=255)                                       ccATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCgggtgg                                                                                                               <  1:117810/150‑1 (MQ=255)                                        cATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCgggtggg                                                                                                              <  1:127815/150‑1 (MQ=255)                                                    aagaTCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCtgatga                                                                                                  <  1:137068/150‑1 (MQ=255)                                                    aagaTCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCtgatga                                                                                                  <  2:860164/150‑1 (MQ=255)                                                     agaTCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGTACTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAg                                                                                                 <  1:1007733/150‑1 (MQ=255)                                                     agaTCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAg                                                                                                 <  2:710854/150‑1 (MQ=255)                                                      gaTCGCGGCCGGGGAGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAgg                                                                                                <  1:293527/150‑1 (MQ=255)                                                                       gCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATaaaa                                                                               <  1:982605/150‑1 (MQ=255)                                                                         cTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATaaaaaa                                                                             <  2:462829/150‑1 (MQ=255)                                                                                aGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGATGTTTGCCAAGGCATAAAAAAGGGcccc                                                                      <  1:1041059/150‑1 (MQ=255)                                                                                       cTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAgg                                                               >  1:241887/1‑150 (MQ=255)                                                                                          ggggTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCaa                                                                                        <  2:114619/122‑1 (MQ=255)                                                                                          gggCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCaa                                                                                        >  1:114619/1‑122 (MQ=255)                                                                                          gggCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGc                                                            <  1:146608/150‑1 (MQ=255)                                                                                             cTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGcccc                                                         >  2:967319/1‑150 (MQ=255)                                                                                                cATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                      <  1:967319/150‑1 (MQ=255)                                                                                                cATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                      <  1:416056/150‑1 (MQ=255)                                                                                                 aTCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTggg                                                     >  1:618292/1‑150 (MQ=255)                                                                                                  tCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGc                                                    >  1:292051/1‑150 (MQ=255)                                                                                                    tACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCt                                                  >  1:1302013/1‑150 (MQ=255)                                                                                                      cAACGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCtt                                                                    <  1:360916/130‑1 (MQ=255)                                                                                                      cAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCtt                                                                    >  2:360916/1‑130 (MQ=255)                                                                                                       aaaGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcg                                               <  1:1211231/150‑1 (MQ=255)                                                                                                         aGGGCCGGCGCTGCTTCGCGAAGCGCTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcgcc                                             >  2:375474/1‑150 (MQ=255)                                                                                                         aGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGGGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcgcc                                             >  2:291662/1‑150 (MQ=255)                                                                                                             ccGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGAAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGCGCCCCTGGGCCTTcgcccgcc                                         <  1:377682/150‑1 (MQ=255)                                                                                                                     gATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGAt                                 <  1:1258543/150‑1 (MQ=255)                                                                                                                      aTTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATg                                <  1:291662/150‑1 (MQ=255)                                                                                                                        tCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATggg                              <  2:1008757/150‑1 (MQ=255)                                                                                                                              aGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGc                        >  2:1092673/1‑150 (MQ=255)                                                                                                                                  gTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATgg                    <  2:618292/150‑1 (MQ=255)                                                                                                                                   tGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGt                   >  2:505063/1‑150 (MQ=255)                                                                                                                                     gATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATaaaaaa                                                                             <  2:1046052/90‑1 (MQ=255)                                                                                                                                     gATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTaa                 >  1:1131651/1‑150 (MQ=255)                                                                                                                                     gATGCTATCGCGGCAATGCCGCGGTGACCAGCACGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATaaaaaa                                                                             >  1:1046052/1‑90 (MQ=255)                                                                                                                                       tGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTGCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATTTGGCGTGCATGGTAAAt               <  1:254273/150‑1 (MQ=255)                                                                                                                                        gCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATc              <  1:375474/150‑1 (MQ=255)                                                                                                                                         tcaccgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCTCGCCGCCCGGATGGGGCGTGCATg                     >  2:337332/5‑142 (MQ=255)                                                                                                                                         tcaccgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCTCGCCGCCCGGATGGGGCGTGCATg                     <  1:337332/138‑1 (MQ=255)                                                                                                                                            tCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAg                                                                                       <  1:980328/73‑1 (MQ=255)                                                                                                                                            tCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAg                                                                                       >  2:980328/1‑73 (MQ=255)                                                                                                                                                ggCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcg                                               >  2:734258/1‑109 (MQ=255)                                                                                                                                                ggCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcg                                               <  1:734258/109‑1 (MQ=255)                                                                                                                                                ggCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAg      <  2:241887/150‑1 (MQ=255)                                                                                                                                                 gCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAgt     <  1:1092673/150‑1 (MQ=255)                                                                                                                                                    aTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTc  <  2:1131651/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      TTGCCGATCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTC  >  NC_002947/2388148‑2388444

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