Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,565,203	G→A	6.0%	S177S (AGC→AGT)	dosP ←	oxygen‑sensing c‑di‑GMP phosphodiesterase DosP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,565,203	0	G	A	6.0%	203.7 / 2.4	66	S177S (AGC→AGT)	dosP	oxygen‑sensing c‑di‑GMP phosphodiesterase DosP
Reads supporting (aligned to +/- strand):  ref base G (24/38);  new base A (2/2);  total (26/40)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

GCTTTAATCCAGATTTTTTCACCGGTGCGCGTCAACAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCC  >  NC_000913/1565057‑1565347                                                                                                                                                   |                                                                                                                                                 gCTTTAATCCAGATTTTTTCACCGGTGCGCGTCAACAGCAGAAATTCGGCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTg                                                                                                                                               <  2:1295052/150‑1 (MQ=255)  cTTTAATCCAGATTTTTTCACCGGTGCGCGTCAACAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTgg                                                                                                                                              <  2:991959/150‑1 (MQ=255)   tttAATCCAGATTTTTTCACCGGTGCGCGTCAACAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGc                                                                                                                                             <  2:1332035/150‑1 (MQ=255)   tttAATCCAGATTTTTTCACCGGTGCGCGTCAACAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGc                                                                                                                                             >  1:612383/1‑150 (MQ=255)     tAATCCAGATTTTTTCACCGGTGCGCGTCAACAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCACTGGCtt                                                                                                                                           >  2:369123/1‑150 (MQ=255)     tAATCCAGATTTTTTCACCGGTGCGCGTCAACAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCACTGGCtt                                                                                                                                           >  2:362086/1‑150 (MQ=255)      aaTCCAGATTTTTTCACCGGTGCGCGTCAACAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTc                                                                                                                                          >  1:1580888/1‑150 (MQ=255)            gATTTTTTCACCGGTGCGCGTCAACAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCAGGCTTCGCTAAt                                                                                                                                    >  2:122094/1‑150 (MQ=255)              ttttttCACCGGTGCGCGTCAACAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGc                                                                                                                                  <  1:1493567/150‑1 (MQ=255)                   tCACCGGTGCGCGTCAACAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATAGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTaa                                                                                                                             >  2:408853/1‑150 (MQ=255)                      ccGGTGCGCGTCAACAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCa                                                                                                                          <  1:1284811/150‑1 (MQ=255)                            cgcgTCAACAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTa                                                                                                                              <  1:739269/140‑1 (MQ=255)                            cgcgTCAACAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTa                                                                                                                              >  2:739269/1‑140 (MQ=255)                                tCAACAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCgg                                                                                                                <  1:395746/150‑1 (MQ=255)                                tCAACAGCAGAAATTCGTCCTCATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCgg                                                                                                                <  2:400421/150‑1 (MQ=255)                                  aaCAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTa                                                                                                              <  1:1668844/150‑1 (MQ=255)                                   aCAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTaa                                                                                                             <  1:563088/150‑1 (MQ=255)                                    cagcagAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACAtt                                                                                                                    <  1:247368/142‑1 (MQ=255)                                    cagcagAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACAtt                                                                                                                    >  2:247368/1‑142 (MQ=255)                                     agcagAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAAt                                                                                                           >  1:13149/1‑150 (MQ=255)                                      gcagAAATTCGACCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAAAGCCGTTATCGGCAGGGATTTCAGGAATGTTCAGGAGTGTATCGCGCTGCAAACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATg                                                                                                          >  2:879669/1‑150 (MQ=255)                                        agAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATgcg                                                                                                        >  2:1465943/1‑150 (MQ=255)                                         gAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACAt                                                                                                                     >  1:1483742/1‑136 (MQ=255)                                         gAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACAt                                                                                                                     <  2:1483742/136‑1 (MQ=255)                                                   cTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATtgcactg                                                                                             <  2:1580888/150‑1 (MQ=255)                                                                gTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGt                                                                                >  2:1496213/1‑150 (MQ=255)                                                                  tttCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTc                                                                              <  1:1373998/150‑1 (MQ=255)                                                                   ttCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTcc                                                                             <  2:1839395/150‑1 (MQ=255)                                                                   ttCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTcc                                                                             <  2:1084314/150‑1 (MQ=255)                                                                    tCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCg                                                                            <  1:1496213/150‑1 (MQ=255)                                                                         aGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCg                                                                       <  2:1716358/150‑1 (MQ=255)                                                                           cAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAg                                                                     <  2:520937/150‑1 (MQ=255)                                                                            aaCTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAgg                                                                    <  2:423878/150‑1 (MQ=255)                                                                              cTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGAc                                                                  <  2:857213/150‑1 (MQ=255)                                                                              cTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGAc                                                                  >  2:1641438/1‑150 (MQ=255)                                                                              cTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGCGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGAc                                                                  <  1:969162/150‑1 (MQ=255)                                                                                  tgtAAACGAATGCGGTTATCGGCAGCGAAATCACGAATGTTAAGGAGTGTAACTGGCTGCAAACCGCTGGCATCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAAAATGGAGTTCCGGATCGAGGACAAta                                                              >  2:132309/1‑149 (MQ=255)                                                                                   gtAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCa                                                             >  2:1852900/1‑150 (MQ=255)                                                                                     aaaCGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCAcc                                                           <  2:612383/150‑1 (MQ=255)                                                                                      aaCGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcg                                                          <  1:799791/150‑1 (MQ=255)                                                                                      aaCGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcg                                                          <  1:408853/150‑1 (MQ=255)                                                                                             gCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCggtcggtc                                                   <  2:432054/150‑1 (MQ=255)                                                                                                   aTCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATg                                                                                    <  1:1620423/111‑1 (MQ=255)                                                                                                   aTCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATg                                                                                    >  2:1620423/1‑111 (MQ=255)                                                                                                    tCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATgg                                            <  1:1580434/150‑1 (MQ=255)                                                                                                      ggCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTc                                          >  1:389109/1‑150 (MQ=255)                                                                                                      ggCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTc                                          >  1:1480890/1‑150 (MQ=255)                                                                                                       gCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCa                                         >  2:629671/1‑150 (MQ=255)                                                                                                         aGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAAc                                       <  1:1022294/150‑1 (MQ=255)                                                                                                            gAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGc                                    >  2:1405203/1‑150 (MQ=255)                                                                                                            gAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGc                                    <  1:338903/150‑1 (MQ=255)                                                                                                            gAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGc                                    <  1:331715/150‑1 (MQ=255)                                                                                                            gAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGc                                    >  1:1224833/1‑150 (MQ=255)                                                                                                               ttCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCaat                                 <  1:879669/150‑1 (MQ=255)                                                                                                               ttCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCCGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCaat                                 <  1:122094/150‑1 (MQ=255)                                                                                                                  aGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCAGTGCACAATATGGCGTTCCGGATCGAgg                                                                    <  2:1983119/112‑1 (MQ=255)                                                                                                                  aGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAgg                                                                    >  1:1983119/1‑112 (MQ=255)                                                                                                                  aGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCaataat                              <  2:13149/150‑1 (MQ=255)                                                                                                                      aTGTGCAGGCGTGTATCGGTCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAAt                          <  2:1223625/150‑1 (MQ=255)                                                                                                                        gTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTg                        >  1:1658894/1‑150 (MQ=255)                                                                                                                           cAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCg                     >  2:244943/1‑150 (MQ=255)                                                                                                                              gAGTGTATCGGGCTGCATACCACTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGt                  <  1:369123/150‑1 (MQ=255)                                                                                                                              gAGTGTATCGGGCTGCATACCACTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGt                  <  1:362086/150‑1 (MQ=255)                                                                                                                                       gggCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCttt         >  1:268850/1‑150 (MQ=255)                                                                                                                                       gggCTGCATACCGCAGGCTTGGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTCCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCttt         <  2:1681264/150‑1 (MQ=255)                                                                                                                                              aTACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGcc  <  2:1224833/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GCTTTAATCCAGATTTTTTCACCGGTGCGCGTCAACAGCAGAAATTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCC  >  NC_000913/1565057‑1565347

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