Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,671,752	A→G	intergenic (+68/‑25)	asr → / → ydgU	acid shock‑inducible periplasmic protein/stationary phase‑induced protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,671,752	0	A	G	98.5%	252.9 / ‑4.1	68	intergenic (+68/‑25)	asr/ydgU	acid shock‑inducible periplasmic protein/stationary phase‑induced protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (32/35);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-01

GCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTAGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCC  >  NC_000913/1671607‑1671899                                                                                                                                                  |                                                                                                                                                    gcgcAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCAcc                                                                                                                                                >  1:37427/1‑151 (MQ=255)     aggcagCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCt                                                                                                                                            <  2:703551/151‑1 (MQ=255)      ggcagCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCtt                                                                                                                                           >  2:576872/1‑151 (MQ=255)            gaagaaACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCgg                                                                                                                                     >  1:593540/1‑151 (MQ=255)             aagaaACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCgg                                                                                                                                     <  1:658518/150‑1 (MQ=255)             aagaaACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGa                                                                                                                                    <  1:734387/151‑1 (MQ=255)              agaaACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCgg                                                                                                                                     >  1:842974/1‑149 (MQ=255)              agaaACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAg                                                                                                                                   <  1:460528/151‑1 (MQ=255)                aaaCACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCGCCGCCTTTAGCGGa                                                                                                                                    <  1:551948/148‑1 (MQ=255)                  acacGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGgt                                                                                                                                >  2:15173/1‑150 (MQ=255)                         aaGAAACACAGCCATCAGCAACCGGCAAAGCCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTAtggt                                                                                                                        <  2:1015894/151‑1 (MQ=255)                            aaaCACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTggg                                                                                                                     <  2:1181451/151‑1 (MQ=255)                             aaCACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGc                                                                                                                    <  1:973025/151‑1 (MQ=255)                              acacAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGACGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGAAACGGTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTggg                                                                                                                     >  1:1027616/1‑149 (MQ=255)                                   gCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCGGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTa                                                                                                               <  1:1184892/150‑1 (MQ=255)                                    ccATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATc                                                                                                             >  2:67808/1‑151 (MQ=255)                                     cATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCg                                                                                                            >  2:983690/1‑151 (MQ=255)                                        cAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCtt                                                                                                         >  2:502256/1‑151 (MQ=255)                                         aGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCttt                                                                                                        <  1:359419/151‑1 (MQ=255)                                          gCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTg                                                                                                       >  2:16407/1‑151 (MQ=255)                                          gCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTg                                                                                                       >  2:1163138/1‑151 (MQ=255)                                             aCCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGt                                                                                                    >  1:706525/1‑151 (MQ=255)                                              ccGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGtt                                                                                                   >  2:382906/1‑151 (MQ=255)                                                   aaaaCCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTc                                                                                              >  1:1092167/1‑151 (MQ=255)                                                         tgctgcACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGatc                                                                                         >  1:1065290/1‑150 (MQ=255)                                                          gctgcACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAGGCGGGTAACGTTTGTCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGatc                                                                                         <  1:4859/149‑1 (MQ=255)                                                            tgcACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATcc                                                                                     >  1:201642/1‑151 (MQ=255)                                                             gcACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGatcatc                                                                                      <  2:1075870/149‑1 (MQ=255)                                                              cacaACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCtt                                                                                   >  1:782145/1‑151 (MQ=255)                                                              cacaACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCt                                                                                    <  1:659027/150‑1 (MQ=255)                                                                caACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTAt                                                                                 >  1:574560/1‑151 (MQ=255)                                                                 aaCCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTAtt                                                                                <  2:37431/151‑1 (MQ=255)                                                                  aCCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTAttt                                                                               <  1:1087688/151‑1 (MQ=255)                                                                   cccGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTAtttt                                                                              <  2:935408/151‑1 (MQ=255)                                                                   cccGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTAtttt                                                                              <  1:402575/151‑1 (MQ=255)                                                                     cGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTAt                                                                            >  2:1018647/1‑151 (MQ=255)                                                                       cAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATgc                                                                          >  2:1045814/1‑151 (MQ=255)                                                                       cAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATgc                                                                          >  2:403670/1‑151 (MQ=255)                                                                        atcGTAAGTTTTCAACAGTAATGCTGCCGCGCCCCCTCGCGCCTCAAAATTACGCTGCTAAGCGTGTAACGTTTGGCACCTCCTTTAGCGGAGGGTATTATGATTGTCCGTTATCGCTTTGATTTCATTCTGATCATCCTTATTTTAtcc                                                                          <  1:348371/148‑3 (MQ=255)                                                                         gCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATgcgc                                                                        >  1:511702/1‑151 (MQ=255)                                                                          cGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATgcgc                                                                        >  1:1046003/1‑150 (MQ=255)                                                                           gTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCAc                                                                      >  1:496546/1‑151 (MQ=255)                                                                                ttttCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGAtt                                                                 >  1:249942/1‑151 (MQ=255)                                                                                   tCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTAc                                                               <  2:1012030/150‑1 (MQ=255)                                                                                    cAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCg                                                             <  2:683508/151‑1 (MQ=255)                                                                                      aCAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGcc                                                           <  1:349956/151‑1 (MQ=255)                                                                                       cAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGccc                                                          <  2:749173/151‑1 (MQ=255)                                                                                           aaTGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGtt                                                       <  2:1176350/150‑1 (MQ=255)                                                                                            aTGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGtttt                                                     <  2:235112/151‑1 (MQ=255)                                                                                            aTGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGtttt                                                     >  2:1173032/1‑151 (MQ=255)                                                                                               cTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTAt                                                  >  1:396927/1‑151 (MQ=255)                                                                                                 ggCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGAGTACCGCCCGTTTTTAt                                                  <  2:1099508/149‑1 (MQ=255)                                                                                                   cgcgcgCCAGCGCGCCTGAAAATAACGGTCCTGAGCTGGTAACGTTTGGCTCCCCCTTTAGCGGCGGGTCGTCTGGTGGGCCGTTAGCGCTTTGAGTTTGTTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTTTCttt                                              <  1:1168892/151‑1 (MQ=255)                                                                                                    gcgcCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTc                                             >  1:796379/1‑151 (MQ=255)                                                                                                              gcgcCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAg                                   <  2:292407/151‑1 (MQ=255)                                                                                                                gcCTGAAAATTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGtt                                 >  2:1143398/1‑151 (MQ=255)                                                                                                                     aaaaTTACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTACATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTg                            >  2:341157/1‑151 (MQ=255)                                                                                                                          tACGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGAttt                        <  2:758543/150‑1 (MQ=255)                                                                                                                           aCGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGAtttt                       <  2:1084860/150‑1 (MQ=255)                                                                                                                           aCGGTGCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTa                      <  2:757007/151‑1 (MQ=255)                                                                                                                                gCTAAGCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTccc                 <  1:490045/151‑1 (MQ=255)                                                                                                                                     gCGGGTAACGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACttt            <  1:200897/151‑1 (MQ=255)                                                                                                                                          ttaCGTTTGGCAGCGCCTTTAGCGTAGGTTAGTATGGTCCGCCGTTAACGCTTTGAGTTCATTCTGATCATGCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATATTTCCAGATTGTAGTTATCTGATTTTACTCCCACTTTCAtc        <  1:94773/148‑1 (MQ=255)                                                                                                                                           aaCGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCAtccc      <  1:1006604/151‑1 (MQ=255)                                                                                                                                           aaCGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCAtcc       <  2:186507/150‑1 (MQ=255)                                                                                                                                           aaCGTTTAGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTGTCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCAtccc      <  2:963386/151‑1 (MQ=255)                                                                                                                                             cGTTTGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCAtcccgt    >  1:393864/1‑151 (MQ=255)                                                                                                                                               tttGGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCAtcccgtcc  >  1:727238/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACAGTAATGCTGGCGCGCCCCCTCGCGCCTGAAAATTACGGTGCTAAGCGGGTAACGTTTAGCACCGCCTTTAGCGGAGGGTAGTATGGTGGGCCGTTATCGCTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCC  >  NC_000913/1671607‑1671899

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