Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,444,709	G→C	5.4%	A56G (GCC→GGC)	nicF ←	maleamate amidohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,444,709	0	G	C	5.4%	234.4 / 3.2	74	A56G (GCC→GGC)	nicF	maleamate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base G (36/34);  new base C (2/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CGGGTTGCCTTCGACCATGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGG  >  NC_002947/4444570‑4444855                                                                                                                                            |                                                                                                                                                    cGGGTTGCCTTCGACCATGTCCTTCATCACCGGCGCCTTGGGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTAGCCGAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTc                                                                                                                                           >  1:99889/1‑150 (MQ=255)      tGCCTTCGACCATGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTg                                                                                                                                      >  2:753912/1‑150 (MQ=255)       gCCTTCGACCATGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGCCGGCCTGCtggcgctg                                                                                                                                     <  1:1612056/150‑8 (MQ=255)        ccTTCGACCATGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGcc                                                                                                                                    <  2:1503997/150‑1 (MQ=255)        ccTTCGACCATGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGcc                                                                                                                                    >  1:1580370/1‑150 (MQ=255)         cTTCGACCATGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCcg                                                                                                                                   >  2:522989/1‑150 (MQ=255)         cTTCGACCATGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCcg                                                                                                                                   >  1:1707611/1‑150 (MQ=255)         cTTAGACCATGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCcg                                                                                                                                   >  2:138632/1‑150 (MQ=255)          ttCGACCATGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCcgc                                                                                                                                  >  2:812054/1‑150 (MQ=255)          ttCGACCATGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCcgc                                                                                                                                  <  1:753912/150‑1 (MQ=255)             gACCATGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑Ac                                                                                                                              <  2:320313/150‑1 (MQ=255)              aCCATGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑Acc                                                                                                                             >  2:525878/1‑150 (MQ=255)               ccATGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCa                                                                                                                            >  2:607924/1‑150 (MQ=255)                  tGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGc                                                                                                                         >  1:1394194/1‑150 (MQ=255)                  tGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGCCGGCCTGCTCGACTGcc                                                                                                                                    >  1:8284/1‑140 (MQ=255)                  tGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGCCGGCCTGCTCGACTGcc                                                                                                                                    <  2:8284/140‑1 (MQ=255)                     ccTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTAGATGGTCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCg                                                                                                                      >  1:1552969/1‑150 (MQ=255)                       ttCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCggg                                                                                                                    <  1:1686810/150‑1 (MQ=255)                                ggCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCCGGGCgcaccgc                                                                                                           >  1:984207/1‑144 (MQ=255)                                    ccTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTg                                                                                                       >  1:204653/1‑150 (MQ=255)                                          gcACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCg                                                                                                 >  1:161770/1‑150 (MQ=255)                                            aCCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGt                                                                                                        >  1:610219/1‑141 (MQ=255)                                            aCCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGt                                                                                                        <  2:610219/141‑1 (MQ=255)                                            aCCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCggg                                                                                               >  2:617956/1‑150 (MQ=255)                                              ccACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGggt                                                                                             >  2:1133483/1‑150 (MQ=255)                                               cacaCGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGggtg                                                                                            <  1:1824138/150‑1 (MQ=255)                                               cacaCGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGggtg                                                                                            >  2:2280/1‑150 (MQ=255)                                                acacGCCGCCGTCAGCAAAGTGCGGAGGCTGGTAGCGGATGTTGGTATGCACCCCCAAGGTACCGCAGTCGCGGGGCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑GCCACGCCCGGGGCGTACAGCGGTGCGCCCGGggtgg                                                                                           <  1:689732/150‑1 (MQ=255)                                                 cacGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGggtggt                                                                                          <  2:1580370/150‑1 (MQ=255)                                                 cacGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACGGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGgggggt                                                                                          >  1:1153106/1‑150 (MQ=255)                                                  acGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAGGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGgtggtg                                                                                         <  1:32691/150‑1 (MQ=255)                                                    gccgccGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCGGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTa                                                                                       <  2:809878/150‑1 (MQ=255)                                                      cgccgTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGc                                                                                     >  2:1169870/1‑150 (MQ=255)                                                           tCAGCAAAATGCGGAGGCTGCTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGc                                                                                >  2:169096/1‑150 (MQ=255)                                                             aGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCAt                                                                              >  2:1415587/1‑150 (MQ=255)                                                                 aaaTGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAg                                                                          >  2:1807406/1‑150 (MQ=255)                                                                 aaaTGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAg                                                                          >  1:106874/1‑150 (MQ=255)                                                                     gCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGa                                                                      <  1:2280/150‑1 (MQ=255)                                                                     gCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGa                                                                      <  1:1912147/150‑1 (MQ=255)                                                                     gCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGa                                                                      <  1:1709396/150‑1 (MQ=255)                                                                       ggAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCGCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGatc                                                                    <  1:607924/150‑1 (MQ=255)                                                                          ggCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGatcatc                                                                 <  2:204653/150‑1 (MQ=255)                                                                            cTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGGCAGCAGGCCGGCGGCCTGCTCGATTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATcag                                                               <  1:1070204/150‑1 (MQ=255)                                                                               gTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATcagcag                                                            >  2:673044/1‑150 (MQ=255)                                                                                   cGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGTGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGcc                                                        <  2:986569/150‑1 (MQ=255)                                                                                   cGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGcc                                                        <  1:522989/150‑1 (MQ=255)                                                                                    ggATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGACCAGCGCCAGCAGGCCGGCGGTCTGCTCGACTGCCGGGTTTCTCG‑ACGGGCGATTGAGCGTTACGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCg                                                       >  1:897000/1‑150 (MQ=255)                                                                                         ttGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTg                                                  <  1:525878/150‑1 (MQ=255)                                                                                          tGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGc                                                 <  1:812054/150‑1 (MQ=255)                                                                                           ggTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTAGACTGCGTCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCg                                                <  2:58736/150‑1 (MQ=255)                                                                                             tATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGaa                                                                           <  2:805858/121‑1 (MQ=255)                                                                                             tATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGaa                                                                           >  1:805858/1‑121 (MQ=255)                                                                                                 caccacCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑GCCACTCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCaaa                                          <  2:355854/150‑1 (MQ=255)                                                                                                  accaccAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGccc                                                                                   <  2:745939/108‑1 (MQ=255)                                                                                                  accaccAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGccc                                                                                   >  1:745939/1‑108 (MQ=255)                                                                                                  accaccAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAg                                         <  1:138632/150‑1 (MQ=255)                                                                                                   caaccaAGGTACCGGAGTCGCGGGCCAGCGCTAGCAGGCCGGCGGGCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGGCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGTCGGCTTGCGGCCAAAGc                                        <  2:610591/148‑1 (MQ=255)                                                                                                    caccaAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGggtggt                                                                                          >  1:949162/1‑99 (MQ=255)                                                                                                    caccaAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCCTACAGCGGTGCGCCCGGggtggt                                                                                          <  2:949162/99‑1 (MQ=255)                                                                                                     accaAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCg                                      >  1:237284/1‑150 (MQ=255)                                                                                                       caAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGAt                                    >  2:695971/1‑150 (MQ=255)                                                                                                         aGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATgc                                  <  1:695971/150‑1 (MQ=255)                                                                                                              ccGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGatcatc                                                                 >  1:1843612/1‑114 (MQ=255)                                                                                                              ccGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGatcatc                                                                 <  2:1843612/114‑1 (MQ=255)                                                                                                                gCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGcc                           >  2:1076474/1‑150 (MQ=255)                                                                                                                   gTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCca                        >  1:1530789/1‑150 (MQ=255)                                                                                                                    tCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCcac                       >  1:1225913/1‑150 (MQ=255)                                                                                                                    ccgcgGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCcac                       <  2:1904097/149‑1 (MQ=255)                                                                                                                          gccagcgccagcAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACAccc                 >  1:1827318/1‑150 (MQ=255)                                                                                                                            cagcgccagcAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTg               <  2:508068/150‑1 (MQ=255)                                                                                                                             agcgccagcAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTgg              <  1:1133483/150‑1 (MQ=255)                                                                                                                                   agcagGCCGCCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTg        >  1:359498/1‑150 (MQ=255)                                                                                                                                      agGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCa     <  2:1707611/150‑1 (MQ=255)                                                                                                                                        gccgccGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACg   >  2:720269/1‑150 (MQ=255)                                                                                                                                         ccGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACgg  <  1:673044/150‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                    CGGGTTGCCTTCGACCATGTCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCG‑ACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGG  >  NC_002947/4444570‑4444855

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