Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,021,663	(A)5→4	coding (1737/2055 nt)	helD →	DNA helicase IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,021,659	0	A	.	98.0%	201.9 / ‑1.5	50	coding (1733/2055 nt)	helD	DNA helicase IV
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (33/16);  total (34/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

TGGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATC  >  NZ_CP009273/1021575‑1021744                                                                                     |                                                                                      tGGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGaaaagc                                                                                  <  2:152013/90‑3 (MQ=255)       aGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAac                                                                           >  2:1566458/1‑90 (MQ=255)       aGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAac                                                                           >  2:895404/1‑90 (MQ=255)       aGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAac                                                                           >  2:586738/1‑90 (MQ=255)            tctGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCt                                                                      >  2:1477782/1‑90 (MQ=255)            tctGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCt                                                                      >  2:1362501/1‑90 (MQ=255)            tctGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCt                                                                      >  2:1008477/1‑90 (MQ=255)                  tATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGa                                                                >  2:478247/1‑90 (MQ=255)                   aTGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGaa                                                               >  2:961018/1‑90 (MQ=255)                    tGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAg                                                              <  1:586738/90‑1 (MQ=255)                    tGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAg                                                              <  1:1566458/90‑1 (MQ=255)                    tGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAg                                                              >  1:688844/1‑90 (MQ=255)                         aaCCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGa                                                                >  2:1193956/1‑83 (MQ=255)                         aaCCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGa                                                                <  1:1193956/83‑1 (MQ=255)                         aaCCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCa                                                         >  1:292237/1‑90 (MQ=255)                         aaCCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCa                                                         >  2:1130960/1‑90 (MQ=255)                             ggAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATc                                                     >  1:816957/1‑90 (MQ=255)                              gAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCg                                                    <  1:143564/90‑1 (MQ=255)                              gAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCg                                                    <  1:895404/90‑1 (MQ=255)                                agagCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGAc                                                  <  2:79009/90‑1 (MQ=255)                                   gcgcATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACttt                                               >  1:1135165/1‑90 (MQ=255)                                   gcgcATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACttt                                               >  2:65947/1‑90 (MQ=255)                                   gcgcATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACttt                                               >  2:1225484/1‑90 (MQ=255)                                            gATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCAtt                                      >  1:627117/1‑90 (MQ=255)                                            gATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCAtt                                      >  1:464631/1‑90 (MQ=255)                                             aTCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTc                                     >  2:1279646/1‑90 (MQ=255)                                                cTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATg                                  >  2:621231/1‑90 (MQ=255)                                                 tGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGc                                 >  1:726731/1‑90 (MQ=255)                                                  ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGc                                 >  1:636536/1‑90 (MQ=255)                                                  ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCg                                                    >  1:1270303/1‑70 (MQ=255)                                                  ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCg                                                    <  2:1270303/70‑1 (MQ=255)                                                  ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGcc                                <  1:287591/90‑1 (MQ=255)                                                  ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGcc                                >  1:1485347/1‑90 (MQ=255)                                                  ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGcc                                >  2:608340/1‑90 (MQ=255)                                                       gTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCaa                           >  2:543938/1‑90 (MQ=255)                                                       gTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCaa                           <  2:1095405/90‑1 (MQ=255)                                                         tcccATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAg                         <  1:278728/88‑1 (MQ=255)                                                              tCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCaa                    <  2:726731/90‑1 (MQ=255)                                                                   tGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGc               <  1:1052310/90‑1 (MQ=255)                                                                   tGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGc               <  1:759476/90‑1 (MQ=255)                                                                     aGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCgg             >  1:516166/1‑90 (MQ=255)                                                                      ggCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGa            >  2:855309/1‑90 (MQ=255)                                                                        ccTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGAtt          >  2:1321692/1‑90 (MQ=255)                                                                        ccTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGAtt          >  1:354090/1‑90 (MQ=255)                                                                        ccTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGAtt          <  2:938314/90‑1 (MQ=255)                                                                         cTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTa         >  2:577526/1‑90 (MQ=255)                                                                         cTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTa         >  1:936338/1‑90 (MQ=255)                                                                           gCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACg       >  1:1134515/1‑90 (MQ=255)                                                                               gCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAAc                   <  2:1303825/74‑1 (MQ=255)                                                                               gCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAAc                   >  1:1303825/1‑74 (MQ=255)                                                                                ccTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGtcatc  <  2:816957/90‑1 (MQ=255)                                                                                     |                                                                                      TGGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATC  >  NZ_CP009273/1021575‑1021744

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