Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,091,449	A→C	5.8%	N20H (AAC→CAC)	garD →	galactarate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,091,449	0	A	C	5.8%	272.4 / 6.5	86	N20H (AAC→CAC)	garD	galactarate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (44/37);  new base C (2/3);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

GCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTGCGCCGCTACGGGCAAAT  >  NC_002947/4091301‑4091597                                                                                                                                                     |                                                                                                                                                     gCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACa                                                                                                                                                      <  5:366705/149‑1 (MQ=255)     ccGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAAcgt                                                                                                                                                  >  5:193789/1‑149 (MQ=255)     ccGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAAcgt                                                                                                                                                  >  7:43922/1‑149 (MQ=255)          aCTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTgg                                                                                                                                             >  5:294610/1‑149 (MQ=255)          aCTCAAAAAACACGCCCGCGCCGGGTCGCCCGGAATGGGCGGATGTACAAACGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACACTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTgg                                                                                                                                             <  1:84414/149‑1 (MQ=255)            tGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGgtc                                                                                                                                           >  3:22478/1‑149 (MQ=255)            tGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGgtc                                                                                                                                           >  3:335151/1‑149 (MQ=255)              aaaaaaCACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACATAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGgtcgt                                                                                                                                         >  1:252554/1‑149 (MQ=255)               aaacccccGCCCGCGCCGGTGCGCCCGGCATGGGCGGATGTACAAACGCCTGAGAACCTTCCAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCCCTACGTCCGCCTGCACGACGACGACAACGTCGTGgtcgtc                                                                                                                                        <  8:317451/142‑1 (MQ=255)                 aaaCACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACCACGACGTTgtcgtcgt                                                                                                                                      <  5:313691/149‑1 (MQ=255)                       gCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAacga                                                                                                                                >  8:242718/1‑149 (MQ=255)                             gcCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGAcggcgg                                                                                                                          <  4:22478/149‑1 (MQ=255)                                   tCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTggg                                                                                                                    <  2:121659/149‑1 (MQ=255)                                   tCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTggg                                                                                                                    <  7:242718/149‑1 (MQ=255)                                       ccGGCAGGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGaa                                                                                                                <  6:193789/149‑1 (MQ=255)                                        cGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAg                                                                                                               <  2:252554/149‑1 (MQ=255)                                         ggCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAgg                                                                                                              <  8:43922/149‑1 (MQ=255)                                            aTGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGgcgc                                                                                                           >  6:366047/1‑149 (MQ=255)                                            aTGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTACGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGAcggc                                                                                                                            >  7:139712/1‑132 (MQ=255)                                            aTGGGCGGATGTACAAAAGCATGAGAACCTGACATGTAGGGCCCAAAAATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGAcggc                                                                                                                            <  8:139712/132‑1 (MQ=255)                                               ggCGGATGTACAAAAGCCAGAGAACCTGACAAGTAGGGCCCAAACATGGAGTTGATCGAACATTCCGACTCGCCCAGCTACGTCCGCCTGCACGACGACGACAACGACGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGacct                                                                                                        >  6:175669/1‑148 (MQ=255)                                                       tACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCg                                                                                                >  8:259338/1‑149 (MQ=255)                                                         cAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGAc                                                                                              <  5:366047/149‑1 (MQ=255)                                                         cAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGAc                                                                                              >  8:156670/1‑149 (MQ=255)                                                         cAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGGGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGAc                                                                                              <  1:192632/149‑1 (MQ=255)                                                          aaaaGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACg                                                                                             >  1:303139/1‑149 (MQ=255)                                                          aaaaGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACg                                                                                             >  1:348982/1‑149 (MQ=255)                                                          aaaaGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACg                                                                                             >  6:193297/1‑149 (MQ=255)                                                          aaaaGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACg                                                                                             >  6:172203/1‑149 (MQ=255)                                                          aaaaGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACCACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAg                                                                                                               <  1:192168/131‑1 (MQ=255)                                                          aaaaGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACCACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAg                                                                                                               >  2:192168/1‑131 (MQ=255)                                                            aaGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGc                                                                                           >  5:194165/1‑149 (MQ=255)                                                              gCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGAcggc                                                                                                                            >  4:282560/1‑114 (MQ=255)                                                              gCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGAcggc                                                                                                                            <  3:282560/114‑1 (MQ=255)                                                                cTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGa                                                                                       <  6:107970/149‑1 (MQ=255)                                                                   agaACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGAccc                                                                                    <  5:172203/149‑1 (MQ=255)                                                                   agaACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGAcc                                                                                     <  2:348982/148‑1 (MQ=255)                                                                     aaCCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCt                                                                                   >  5:38305/1‑148 (MQ=255)                                                                       ccTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGAt                                                                                >  1:141823/1‑149 (MQ=255)                                                                         tGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGaa                                                                                                                >  7:276508/1‑115 (MQ=255)                                                                         tGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGaa                                                                                                                <  8:276508/115‑1 (MQ=255)                                                                          gACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACCACGACGTTgtcgtcgtc                                                                                                                                     <  2:68926/93‑1 (MQ=255)                                                                          gACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCACCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCG‑CCGCTTctc                                                                                                  <  7:230088/127‑3 (MQ=255)                                                                          gACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACGCGCCCCGCTACGTCCGCCTGCACGACGACGACCACGACGTTgtcgtcgtc                                                                                                                                     >  1:68926/1‑93 (MQ=38)                                                                             aaGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAgg                                                                          >  6:243891/1‑149 (MQ=255)                                                                                       cAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGAccc                                                                                    >  1:155837/1‑129 (MQ=255)                                                                                       cAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGAccc                                                                                    <  2:155837/129‑1 (MQ=255)                                                                                        aaaCATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCaaa                                                               <  2:303139/149‑1 (MQ=255)                                                                                         aaCATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAg                                                              <  6:38305/149‑1 (MQ=255)                                                                                          aCATGCAGTTGATCGAACTTTCCGACGCGCCCCGCTACGTCCGCCTGCAGGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCaaatc                                                             <  5:175669/149‑3 (MQ=255)                                                                                           cATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGcc                                                            <  5:243891/149‑1 (MQ=255)                                                                                           cATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGcc                                                            <  6:194165/149‑1 (MQ=255)                                                                                             tGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCa                                                           >  3:118648/1‑148 (MQ=255)                                                                                               cAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATa                                                         <  4:118648/148‑1 (MQ=255)                                                                                                 gTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGGCCCTGAGCGAAGGCGTGCCGCAAGGccgggagg                                                      >  8:227773/1‑143 (MQ=255)                                                                                                 gTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAgg                                                      >  2:51813/1‑149 (MQ=255)                                                                                                 gTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAgg                                                      >  8:229307/1‑149 (MQ=255)                                                                                                  ttGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGt                                                     <  7:229307/149‑1 (MQ=255)                                                                                                    gATCGAACATTCCGACTCGCCCCGCTACGTGCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTgg                                                   >  1:29671/1‑149 (MQ=255)                                                                                                    gATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTgg                                                   >  4:148416/1‑149 (MQ=255)                                                                                                    gATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTg                                                    <  1:51813/148‑1 (MQ=255)                                                                                                     aTCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGaa                                                                            <  3:31408/123‑1 (MQ=255)                                                                                                     aTCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGaa                                                                            >  4:31408/1‑123 (MQ=255)                                                                                                         aaCATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCaaa                                                               <  7:252249/132‑1 (MQ=255)                                                                                                         aaCATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCaaa                                                               >  8:252249/1‑132 (MQ=255)                                                                                                         aaCATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCAcc                                              >  1:240374/1‑149 (MQ=255)                                                                                                            aTTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCgtg                                           >  8:184298/1‑149 (MQ=255)                                                                                                                  aCTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATc                                     <  7:259338/149‑1 (MQ=255)                                                                                                                   cTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCg                                    <  4:335151/149‑1 (MQ=255)                                                                                                                      gCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCgtgt                                          <  3:2244/140‑1 (MQ=255)                                                                                                                      gCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGcc                                  >  2:284963/1‑148 (MQ=255)                                                                                                                      gCCCCGCTACGTCCGCCTGCACGACGACGAAAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCgtgt                                          >  4:2244/1‑140 (MQ=255)                                                                                                                       ccccGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGcc                                  >  2:162873/1‑147 (MQ=255)                                                                                                                           gCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGc                            <  7:156670/149‑1 (MQ=255)                                                                                                                           gCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGc                            >  2:69038/1‑149 (MQ=255)                                                                                                                                gTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGcc                       >  7:156248/1‑149 (MQ=255)                                                                                                                                    gCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAg                               <  3:231242/137‑1 (MQ=255)                                                                                                                                    gCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAg                               >  4:231242/1‑137 (MQ=255)                                                                                                                                    gCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTg                   >  6:275989/1‑149 (MQ=255)                                                                                                                                          ttgatgacgacCACGACGTTGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTGcgccgc             >  2:88190/6‑149 (MQ=255)                                                                                                                                           cgacgacgaCAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTGCGCCGCt            <  1:162873/149‑1 (MQ=255)                                                                                                                                               gacgacAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTGCGCCGCTACgg        <  5:146140/149‑1 (MQ=255)                                                                                                                                                acgacCACGACGTTGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTGCGCCGCTACgg        <  1:88190/148‑1 (MQ=255)                                                                                                                                                 cgacAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTGCGCCGCTACggg       >  2:60957/1‑148 (MQ=255)                                                                                                                                                  gacAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTGCGCCGCTACGGGc      >  1:29130/1‑148 (MQ=255)                                                                                                                                                    cAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCg                                                                              <  3:91603/74‑1 (MQ=255)                                                                                                                                                    cAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCg                                                                              >  4:91603/1‑74 (MQ=255)                                                                                                                                                     aaCGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTGCGCCGCTACGGGCAAat  <  1:69038/149‑1 (MQ=255)                                                                                                                                                     aaCGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTGCGCCGCTACGGGCAAat  <  7:193634/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTGCGCCGCTACGGGCAAAT  >  NC_002947/4091301‑4091597

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