Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,375,488	A→T	10.0%	V31D (GTC→GAC)	PP_5637 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,375,488	0	A	T	10.0%	194.2 / 12.0	70	V31D (GTC→GAC)	PP_5637	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/28);  new base T (5/2);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

ACCACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGTATA  >  NC_002947/4375339‑4375634                                                                                                                                                      |                                                                                                                                                   accatGATTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGa                                                                                                                                                    >  2:1064364/1‑150 (MQ=255) accacGACTGAGCATATAGCACGTCTTCCACCTCTTCTACACCGTTAATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGCGCCCAGGCCTAGGGTTAAACACCCAGGCCTGGGTCGCGa                                                                                                                                                    <  2:1180795/150‑1 (MQ=255)    acGACTGCGCATATAGCACGCCTTCGATCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTc                                                                                                                                                 >  1:244904/1‑150 (MQ=255)          gcgcATATAGCACGCCTTCGATCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCa                                                                                                                                           <  2:244904/150‑1 (MQ=255)           cgcATATAGCACGCCTTCTACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCag                                                                                                                                          <  2:280663/150‑1 (MQ=255)           cgcATATAGCACGCCTTCTACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCag                                                                                                                                          <  2:281685/150‑1 (MQ=255)                ataGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCt                                                                                                                                     >  2:834210/1‑150 (MQ=255)                         cTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGCGCTGGGCCATCACGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGAt                                                                                                                            >  1:1287079/1‑150 (MQ=255)                         cTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGAt                                                                                                                            >  1:335397/1‑150 (MQ=255)                          ttCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATg                                                                                                                           >  1:721618/1‑150 (MQ=255)                            cGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGt                                                                                                                         <  1:284694/150‑1 (MQ=255)                                  cttctACCCCGTTGATGTTGATGCCAAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCaa                                                                                                                   >  1:194064/1‑150 (MQ=255)                                   ttctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAt                                                                                                                  <  2:4618/150‑1 (MQ=255)                                   ttctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATATCGGTCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAt                                                                                                                  >  1:1028995/1‑150 (MQ=255)                                      tACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAtgtg                                                                                                               >  1:1401169/1‑150 (MQ=255)                                      tACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAtgtg                                                                                                               >  1:1209053/1‑150 (MQ=255)                                        ccccGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGa                                                                                                             >  2:389740/1‑150 (MQ=255)                                        ccccGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGa                                                                                                             <  2:335397/150‑1 (MQ=255)                                         cccGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGAc                                                                                                            >  1:553502/1‑150 (MQ=255)                                                     gatgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAg                                                                                                >  2:1030736/1‑150 (MQ=255)                                                         ccGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGc                                                                                            <  1:389740/150‑1 (MQ=255)                                                         ccGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGc                                                                                            <  2:553502/150‑1 (MQ=255)                                                           gAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTg                                                                                          >  2:68415/1‑150 (MQ=255)                                                           gAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTg                                                                                          <  2:721618/150‑1 (MQ=255)                                                                  gCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTc                                                                                   <  1:834210/150‑1 (MQ=255)                                                                  gCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCCGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTc                                                                                   <  2:1028995/150‑1 (MQ=255)                                                                   ccATGCCGTTGCCCTTCGCATCGTGCAGCCGGGGGATAGTGTCGGACCCAGGAGTAGGGTTAAACACCCAGGCCACCGTCGCGACTGGACCCAGAGGCTCGCTGTGAGGGTCGCCAATGGGGAGGTCCGCGCGGAGGGGCTGGACTGTcc                                                                                  >  2:118348/1‑150 (MQ=255)                                                                   ccATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTcc                                                                                  <  2:1401169/150‑1 (MQ=255)                                                                   ccATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTcc                                                                                  >  2:1288238/1‑150 (MQ=255)                                                                   ccAGGCCGTTGACCTTCGAATCGTGAAGCCGGGGGATAATGTCGGGCCCAGGAGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTAGACCAAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCTGAGGGGCTGGACTTTcc                                                                                  <  2:246077/150‑1 (MQ=255)                                                                         cGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGc                                                                            >  2:685520/1‑150 (MQ=255)                                                                         cGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGc                                                                            >  2:280432/1‑150 (MQ=255)                                                                               ccTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTg                                                                      >  1:666040/1‑150 (MQ=255)                                                                                      aTCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCaa                                                               <  2:1287079/150‑1 (MQ=255)                                                                                      aTCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCaa                                                               <  1:1064364/150‑1 (MQ=255)                                                                                       tCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCt                                                                                           <  2:694397/121‑1 (MQ=255)                                                                                       tCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCt                                                                                           >  1:694397/1‑121 (MQ=255)                                                                                              gCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTCAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGc                                                       <  2:666040/150‑1 (MQ=255)                                                                                                 gggggATAATGTCGGGCCCAGGCGTAGGGTTAAACACGCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGTGGAGGGGCTGGACTTTCCGTAGGCGGTGTGAGGGCAAGGCAACGCCAt                                                    >  2:1385303/1‑150 (MQ=255)                                                                                                  ggggATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATg                                                   >  1:570547/1‑150 (MQ=255)                                                                                                         aTGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCgcgg                                            <  2:570547/150‑1 (MQ=255)                                                                                                            tCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTTCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAg                                                                    <  1:1151550/123‑1 (MQ=255)                                                                                                            tCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTa                                                                               <  1:583815/112‑1 (MQ=255)                                                                                                            tCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTa                                                                               >  2:583815/1‑112 (MQ=255)                                                                                                            tCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAg                                                                    >  2:1151550/1‑123 (MQ=255)                                                                                                            tCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGt                                         >  2:1030455/1‑150 (MQ=255)                                                                                                              gggCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCg                                       >  2:1216149/1‑150 (MQ=255)                                                                                                              gggCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCg                                       >  2:219688/1‑150 (MQ=255)                                                                                                                gCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGaa                                     <  2:1209053/150‑1 (MQ=255)                                                                                                                   cAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTAGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACga                                  <  2:194064/150‑1 (MQ=255)                                                                                                                      gCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCCCGCCATGCTCGCGGCGTCGAACGAgaa                               <  1:1288238/150‑1 (MQ=255)                                                                                                                          aGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaaggaa                           <  1:280432/150‑1 (MQ=255)                                                                                                                             gTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTa                        <  1:685520/150‑1 (MQ=255)                                                                                                                             gTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGACGACCGAgaaggaagga                        >  2:812582/1‑148 (MQ=255)                                                                                                                              ttAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTAc                       >  1:543987/1‑150 (MQ=255)                                                                                                                                aaaCACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACAt                     >  1:649877/1‑150 (MQ=255)                                                                                                                                    acCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATgggg                 >  1:1108835/1‑150 (MQ=255)                                                                                                                                    acCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCgg                 >  2:494932/1‑150 (MQ=255)                                                                                                                                     cccAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGc                >  1:1354002/1‑150 (MQ=255)                                                                                                                                      ccAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGc                                                  <  2:1135092/115‑1 (MQ=255)                                                                                                                                      ccAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGc                                                  >  1:1135092/1‑115 (MQ=255)                                                                                                                                      ccAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCa               >  1:12925/1‑150 (MQ=255)                                                                                                                                          ctcTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaag                              <  1:1415652/129‑1 (MQ=255)                                                                                                                                          ctcTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaag                              >  2:1415652/3‑131 (MQ=255)                                                                                                                                           tcTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATg                    <  2:402257/139‑1 (MQ=255)                                                                                                                                           tcTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATg                    >  1:402257/2‑140 (MQ=255)                                                                                                                                           tcTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAAt          >  2:525495/2‑150 (MQ=255)                                                                                                                                           tcTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAAt          >  2:925976/2‑150 (MQ=255)                                                                                                                                             tGGGTCGCGACTCGACCCAGAGGCTTGATGTGATAGTAAGAGGCATTTTCGTCCTTGCTGAGGGGCTGGACTTTTCTTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATAc        <  1:1385303/150‑1 (MQ=255)                                                                                                                                              gggTCGAGTCGCGACCCAGAGGCTCGCCGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACt       >  2:483464/1‑150 (MQ=255)                                                                                                                                                   gcgACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGtata  >  2:477960/1‑150 (MQ=255)                                                                                                                                                   gcgACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGtata  <  2:12925/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   ACCACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGTATA  >  NC_002947/4375339‑4375634

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