Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,585,015	0	T	A	20.3%	163.7 / 11.4	59	I86K (ATA→AAA)	yrhB	putative heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base T (39/8);  new base A (0/12);  total (39/20)
Fisher's exact test for biased strand distribution p-value = 1.13e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTATCACGACGCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAG  >  NC_000913/3584767‑3585257                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      ttATCACGACGCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATa                                                                                                                                                                                                                                                      >  2:504545/1‑250 (MQ=255) ttATCACGACGCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATa                                                                                                                                                                                                                                                      >  2:452964/1‑250 (MQ=255)   aTCACGACGCATTCGCGAAAGCGAACCATTGCCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACAAGGACGCGTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTTGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATaaa                                                                                                                                                                                                                                                    <  1:487037/250‑1 (MQ=255)   aTCACGACGCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATaaaa                                                                                                                                                                                                                                                   <  1:41642/251‑1 (MQ=255)   aTCACGACGCAGTCGCGAAAGCGGACCATTACCTTGATGATGCAGATCTCCTGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTTTTTCGAAGCACGAGAATTTCTCTAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTAGTGATAAAGACAGTGGTGAAATTCATTCTCTGGTAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGTATTACGAAATaaaa                                                                                                                                                                                                                                                   <  1:268956/251‑1 (MQ=255)                   gAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCt                                                                                                                                                                                                                                 >  1:525274/1‑251 (MQ=255)                    aaaGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCtt                                                                                                                                                                                                                                <  2:401579/251‑1 (MQ=255)                       gCGAACCATTCCCTTGATGATGCAGATCTCCCGGACGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGCCAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAAGC‑CACCTT‑AGGCTTGcc                                                                                                                                                                                                                             >  2:420378/1‑251 (MQ=255)                          aaCCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTg                                                                                                                                                                                                                          >  1:342104/1‑251 (MQ=255)                                tACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACaaa                                                                                                                                                                                                                     >  1:373357/1‑250 (MQ=255)                                         gatgCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCctct                                                                                                                                                                                                           >  1:323953/1‑251 (MQ=255)                                               gATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCAc                                                                                                                                                                                                     >  2:245354/1‑251 (MQ=255)                                                  ctcCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGa                                                                                                                                                                                                   >  2:89062/1‑250 (MQ=255)                                                    cccGGTCGTCATTACTCTACATGGTCGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGCATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTATCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCGCGaa                                                                                                                                                                                                  <  2:355904/249‑1 (MQ=255)                                                     ccGGTCGTCCTTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTCCGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAACATAGACCTCTTCCCACGAag                                                                                                                                                                                                 >  1:535473/1‑249 (MQ=255)                                                          cgtcATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAGTTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGccc                                                                                                                                                                                          >  2:25519/1‑250 (MQ=255)                                                                  ctctACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCtt                                                                                                                                                                                  >  1:257024/1‑251 (MQ=255)                                                                  ctctACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCtt                                                                                                                                                                                  >  1:57800/1‑251 (MQ=255)                                                                    ctACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCtt                                                                                                                                                                                  >  1:280914/1‑249 (MQ=255)                                                                    ctACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCtt                                                                                                                                                                                  <  1:281724/249‑1 (MQ=255)                                                                    ctACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCtt                                                                                                                                                                                  >  1:194716/1‑249 (MQ=255)                                                                    ctACATGGACGCGTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCTCGAGAATTTCTCGAATCTGGAGATGAGGCCGCGCGCGTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTAGAAGAATATCTACAGGATTACGAAATAAAAAAGTC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCtt                                                                                                                                                                                  <  2:464234/249‑1 (MQ=255)                                                                      aCATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑Gt                                                                                                                                                                              >  2:318842/1‑250 (MQ=255)                                                                       cATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTAc                                                                                                                                                                            <  1:496649/251‑1 (MQ=255)                                                                           gACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCccg                                                                                                                                                                        >  1:269264/1‑251 (MQ=255)                                                                             cGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCccgcc                                                                                                                                                                      >  2:236305/1‑251 (MQ=255)                                                                               cTTTAGACAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGt                                                                                                                                                                    >  2:484910/1‑251 (MQ=255)                                                                                    gCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAat                                                                                                                                                               >  2:91478/1‑251 (MQ=255)                                                                                      cAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAatc                                                                                                                                                              >  1:260871/1‑250 (MQ=255)                                                                                            tGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCAcc                                                                                                                                                       >  1:266066/1‑251 (MQ=255)                                                                                                     tGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCggg                                                                                                                                              >  1:75980/1‑251 (MQ=255)                                                                                                       tttCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTc                                                                                                                                            >  1:148060/1‑251 (MQ=255)                                                                                                            aaccaagagaATTTCGCGAAAGGGGAGGGGGGGCCGCGCGCTTAGCTGGTAACGCACTTTTTATTAATGAAAAAAACGGTGGTAAAATTCATTCTCGGGAAACGGAAAAC‑CC‑GAGGGAAAAATATCTAAAGTATAACAAAAAAAAAAAGGC‑TCCCTT‑CGGCTTGCCCTAAAAAAATCGCCCTCTCCCCAAGAAGGGGGCCGCTATCCTTA‑GCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATg                                                                                                                                       <  2:373357/245‑1 (MQ=255)                                                                                                                  gagaATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                                                 >  1:426786/1‑251 (MQ=255)                                                                                                                         tctcGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTc                                                                                                                          >  1:229932/1‑251 (MQ=255)                                                                                                                                      gaTGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                                              >  1:381619/1‑250 (MQ=255)                                                                                                                                              gccgcgttaaGCTGGTAACGCTCTTTTTATTATTGAAAAAGACAGTGGTAAATTTCTTTCTCTGGAAACGGAAAAC‑CC‑GCGGGAAGAATATCTACAGGATTAAAAAAAAAAAAAGCC‑TCCCTT‑CGGTTTCCCCTGACAAAATCCCCCCCTCCCCACGAAGGGGCCCGCTACCCTTA‑TCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAAt                                                                                                     <  2:229932/242‑1 (MQ=255)                                                                                                                                                  cgcttagcgtgaaacccactttttATAATTGAAAAAGACAGGGGTAAAATTCATTCGC‑GGGAAACGCAAAA‑CCCGCGGAAAGAATATCTAAAGAATAAAAAAAAAAAAAAGGC‑TCCCTT‑GGGCTTGCCCTGAAAAAATACCCCCCTCCCCACAAAGAGGGCCGTTACCCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                                                 <  2:478988/231‑1 (MQ=255)                                                                                                                                                         cTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                                                          >  1:116721/1‑251 (MQ=255)                                                                                                                                                            gTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                                                       >  1:62443/1‑251 (MQ=255)                                                                                                                                                                 gCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATc                                                                                  >  2:423596/1‑251 (MQ=255)                                                                                                                                                                       tttATTATTGATAACGACAGTGATGAAATTCATTCTCTGGGCACGGCAAAC‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGACGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctct                                                                            >  2:478945/1‑251 (MQ=255)                                                                                                                                                                              ttGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcg                                                                     >  1:445500/1‑251 (MQ=255)                                                                                                                                                                               tgaaaaaGACAGTGGTAAAATTCATTCTCGGGGAACGACAAAC‑CG‑GTTGGAAGAATATCTAAAGGATTACAAAAAAAAAAAGCC‑TCCCTT‑CGGCTGCCCCTGAAAAAATAGCCCTCCTCCCACGAAGAGGCCCGCTAACTTAA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcg                                                                     <  1:300079/246‑1 (MQ=255)                                                                                                                                                                                          ggggaaaaaTTCTTTCTCTGGGAACGGAAAAC‑CC‑GCTGAAAAAATACCAACAGGATTAAAAAAAAAAAAAGCC‑TCCCTC‑CGGCTCCCCCAAAAAAAATCCCCCTCCCCCCACGAAGAGGGCCGCTACCCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGAc                                                         <  1:91478/245‑1 (MQ=255)                                                                                                                                                                                            ggtgAAATTCTTTCGC‑GGGGACAGCAAAAGCG‑GCTGAAAGAATATCAACAGGATTACAAAAAAAAAAAGTC‑T‑CCTTGCGGCTTCCCCAGAAAAAATGGCCCCCTTCCCAGGAAGGGGGCCGCAACCCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCt                                                       <  1:504545/251‑1 (MQ=255)                                                                                                                                                                                                 aaTTCATTCTCGGGGAACGGCAAAA‑CC‑GCTGGAAAAATATCTACAGGATTAAAAAAAAAAAAAGGC‑TACCTG‑CGGCTTCCCCTGAAAAAATCCCCCTCTCCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGc                                                  <  1:245354/251‑1 (MQ=255)                                                                                                                                                                                                     cattcTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAACCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                                              >  1:204822/1‑251 (MQ=255)                                                                                                                                                                                                       ttcTCCGGGAACGGCAAAA‑CA‑GCTGGAAGAATCTCTACAGGATTACGAAATAATAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTACCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTACATCATCCACAGAGCGCGGATCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTGGTGTACTGCCCATCCCCCCTTTCAGCCCCACTTTCTGACCTTTTGCTTcccg                                            >  2:384844/1‑249 (MQ=255)                                                                                                                                                                                                        tctcTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTTCTTCcagc                                           >  1:166522/1‑251 (MQ=255)                                                                                                                                                                                                                   ggaaaaa‑CC‑GCGGAAAGAATATCAAGAGGATAAAAAAAAAAAAAAGCC‑TCCCTT‑CGGCTTGCCCTAAAAAAACCCCCCTCTCCCCAAGAGGGGGCCCGCTA‑CCTTAAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGt                                <  2:260871/248‑1 (MQ=255)                                                                                                                                                                                                                      aaaa‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATc                             >  1:313040/1‑251 (MQ=255)                                                                                                                                                                                                                              gttgaaaaaaTATCTACAGGATTAAAAAAAAAAAAGGGC‑CACCTT‑CGGCTCCCCCGAACAAAATACCCCTCTCCCCAGGAAGAGGCCCGCTACCCTAA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCGCTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCt                       <  1:109236/243‑1 (MQ=255)                                                                                                                                                                                                                                 ggAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTATACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTaa                     >  1:264439/1‑250 (MQ=255)                                                                                                                                                                                                                                     gaaTATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCt                >  2:308416/1‑251 (MQ=255)                                                                                                                                                                                                                                     aaaTATCAACAGGATTA‑AAAAAAAAAAAGGCCTCCCTT‑CGGCTCCCCCAGACAAAAAACCCCCCTCCCCCCGAAGGGGGCCCCTACCCTTG‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCt                <  1:89062/250‑1 (MQ=255)                                                                                                                                                                                                                                        tatCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCtttt             >  1:251901/1‑251 (MQ=255)                                                                                                                                                                                                                                          cctaaagatttACAAAAAAAAAAAGGC‑TCCCTT‑CGGCTCGCCCGGAAAAAATCCCCCTCTCCCCACGAGGGGGCCCGCTAACTTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTc           <  1:418679/242‑1 (MQ=255)                                                                                                                                                                                                                                                    taagAAAAAAAAAAGGC‑CCCCTC‑CGGCCTGCCCTAAAAAAAGCGCCCTCCCCCCACGAAGGGGGCCGCAACCCTTT‑GTCCCCCGCCCTTAAATCATCCCCCGAGCGCGGGGCGGATGCGCCGTTCAACTCCCCGTCCGGCCCAACCATTAAGGTTTGTGTTCCGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAg  <  1:198122/247‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      TTATCACGACGCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CC‑GCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTT‑CGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAG  >  NC_000913/3584767‑3585257

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