Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,580,095	0	A	C	24.1%	126.2 / 10.3	58	G456G (GGT→GGG)	hsdS	type I restriction enzyme EcoKI specificity protein
Reads supporting (aligned to +/- strand):  ref base A (31/13);  new base C (0/14);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 1.98e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAGGATGGA  >  NC_000913/4579957‑4580229                                                                                                                                           |                                                                                                                                        aaaaTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAc                                                                                                                                        >  1:453776/1‑140 (MQ=255)  aaaTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAcc                                                                                                                                       <  2:1376108/140‑1 (MQ=255)  aaaTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAcc                                                                                                                                       <  2:453776/140‑1 (MQ=255)   aaTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAccc                                                                                                                                      >  2:1262610/1‑140 (MQ=255)   aaTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAccc                                                                                                                                      >  2:1507742/1‑140 (MQ=255)               tGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTgcg                                                                                                                          >  1:1137754/1‑140 (MQ=255)                  gCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGt                                                                                                                       >  1:672258/1‑140 (MQ=255)                                           atTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AAcg                                                                                             >  2:620047/1‑140 (MQ=255)                                             ttATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AAcgcg                                                                                           >  2:1510919/1‑140 (MQ=255)                                                tgtgGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACgcggcg                                                                                        <  2:1332026/140‑1 (MQ=255)                                                tgtgGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACgcggcg                                                                                        <  2:1311393/140‑1 (MQ=255)                                                  tgGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACgcggcggc                                                                                      <  2:300228/140‑1 (MQ=255)                                                    ggTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGgcgc                                                                                    >  2:950093/1‑140 (MQ=255)                                                    ggTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGgcgc                                                                                    >  2:575043/1‑140 (MQ=255)                                                    ggTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGGGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGgcgc                                                                                    >  1:1127484/1‑140 (MQ=255)                                                      ttATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTg                                                                                  >  1:714680/1‑140 (MQ=255)                                                      ttATGATTGGCGTGAAATAAT‑AAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCGCCTAGCTGGGCGTTCACCTTTG‑TTTTTTCCCGCAAACGCGGCGGCGCTGt                                                                                 >  2:217322/1‑140 (MQ=255)                                                          gATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGtttt                                                                              >  1:1217108/1‑140 (MQ=255)                                                              ggCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTcc                                                                          >  2:346789/1‑140 (MQ=255)                                                               gCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCg                                                                         <  2:628245/140‑1 (MQ=255)                                                               gCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCg                                                                         >  2:91793/1‑140 (MQ=255)                                                               gCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCg                                                                         <  2:352755/140‑1 (MQ=255)                                                                cGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGc                                                                        >  1:1512161/1‑140 (MQ=255)                                                                          taaAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATc                                                              >  1:132426/1‑140 (MQ=255)                                                                                 cgcACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGttt                                                       >  2:873954/1‑140 (MQ=255)                                                                                  ggcccGGAAATGTGCCCCAAGAAATAATTTTCCGGATTTTTTACGTGGGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGtttt                                                      <  1:217322/137‑1 (MQ=255)                                                                                  gcACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGtttt                                                      <  1:553215/140‑1 (MQ=255)                                                                                  gcACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGtttt                                                      <  2:903447/140‑1 (MQ=255)                                                                                   cACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTc                                                     >  1:511898/1‑140 (MQ=255)                                                                                       ggAAGGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGGTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGcc                                                 <  1:991915/140‑1 (MQ=255)                                                                                       ggAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGTTTTTCGGcc                                                 <  2:672258/140‑1 (MQ=255)                                                                                              ggcgcCAAAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACt                                          <  2:1087993/139‑1 (MQ=255)                                                                                                cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGTTTTTTTACCCCCGCTAGCTGCGCTTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCCCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg                                        <  1:1052596/140‑1 (MQ=255)                                                                                                cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGTTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGTTTTTCGGCCCGCCACTgg                                        <  2:1512161/140‑1 (MQ=255)                                                                                                cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg                                        >  2:954839/1‑140 (MQ=255)                                                                                                cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg                                        >  1:1031328/1‑140 (MQ=255)                                                                                                cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg                                        >  2:1273671/1‑140 (MQ=255)                                                                                                cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg                                        >  1:1259034/1‑140 (MQ=255)                                                                                                  ccAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGTTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTTTTTTCTCCGCTGATCAAATCCGGTTTTTCGGCCCGCCACTGGGc                                      <  2:1137754/140‑1 (MQ=255)                                                                                                  ccAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGtttt                                                      >  2:697011/1‑124 (MQ=255)                                                                                                  ccAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGtttt                                                      <  1:697011/124‑1 (MQ=255)                                                                                                        aataatGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAg                                >  1:1096777/1‑140 (MQ=255)                                                                                                            atGTTCAGGATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTca                            <  1:135791/140‑1 (MQ=255)                                                                                                                  aGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACgg                      >  1:1015179/1‑140 (MQ=255)                                                                                                                   ggATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGTTTTCCGGCCCGCCACTGGGCGGTAAGTTCACCACGGa                     <  2:261785/140‑1 (MQ=255)                                                                                                                    gATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGaa                    >  1:1526612/1‑140 (MQ=255)                                                                                                                      ttttttACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACg                  <  1:1510919/140‑1 (MQ=255)                                                                                                                        ttttAGGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑ACCGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCGGGTTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCt                <  2:959267/140‑1 (MQ=255)                                                                                                                        ttttACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCt                <  2:1031328/140‑1 (MQ=255)                                                                                                                         tttACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCtt               >  2:1046766/1‑140 (MQ=255)                                                                                                                            aCGTGAGGTTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTCCCAGC‑ACCGCGGCGGCGCTTTTTTCTCCGCTGATCAAATCGGGTTTTTGGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTg            <  2:175368/140‑1 (MQ=255)                                                                                                                              gTGAGGTTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTCCCAGC‑AACGCGGCGGCGCTTTTTTCTCCGCTGATCAAACCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGcc          <  2:1015179/140‑1 (MQ=255)                                                                                                                              gTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGTTTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGTTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGcc          <  2:256785/140‑1 (MQ=255)                                                                                                                              gTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGcc          <  1:131161/140‑1 (MQ=255)                                                                                                                              gTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTCCCAGC‑ACCGCGGCGGCGCTGTTTTCCCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGcc          <  1:156811/140‑1 (MQ=255)                                                                                                                              gTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGcc          >  1:1132983/1‑140 (MQ=255)                                                                                                                              gTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGcc          >  2:131161/1‑140 (MQ=255)                                                                                                                                 gggCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAgg       <  1:324176/139‑1 (MQ=255)                                                                                                                                   gCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAggat     >  1:784070/1‑140 (MQ=255)                                                                                                                                      tttttCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAggatgga  <  1:557764/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                        AAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGC‑AACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAGGATGGA  >  NC_000913/4579957‑4580229

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