BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	660,180	A→G	6.4%	Y13H (TAC→CAC)	lipA ←	lipoyl synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	660,180	0	A	G	6.4%	429.1 / 13.7	141	Y13H (TAC→CAC)	lipA	lipoyl synthase
Reads supporting (aligned to +/- strand): ref base A (58/74); new base G (7/2); total (65/76)
Fisher's exact test for biased strand distribution p-value = 8.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 146 total aligned reads displayed.

TTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTATATCATTGAAACGGACC  >  NC_000913/660032‑660329

     gcgcATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAac                                                                                                                                                 <  1:670101/150‑1 (MQ=255)

        cATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTAGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACAcc                                                                                                                                              >  2:1908689/1‑150 (MQ=255)

         aTTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACcg                                                                                                                                             >  1:2317815/1‑150 (MQ=255)

         aTTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACcg                                                                                                                                             <  1:1520251/150‑1 (MQ=255)

          ttGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACcgc                                                                                                                                            <  1:1691910/150‑1 (MQ=255)

            gCGGCTTTGATGCCCTGGATACGTGTAGATTACGCTGGAAGCATTATATTAATAAAATAATGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGt                                                                                                                                          >  2:919308/1‑150 (MQ=255)

               gCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTcc                                                                                                                                       >  2:1856258/1‑150 (MQ=255)

               gCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACAGCGTTcc                                                                                                                                       >  1:2362809/1‑150 (MQ=255)

                cTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTGTTAAATACcgcgatgcc                                                                                                                                      <  1:420924/150‑6 (MQ=255)

                  ttGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATc                                                                                                                                    <  1:424578/150‑1 (MQ=255)

                  ttGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATc                                                                                                                                    <  1:49715/150‑1 (MQ=255)

                   tGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATca                                                                                                                                   >  1:806318/1‑150 (MQ=255)

                   tGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATca                                                                                                                                   >  2:1121427/1‑150 (MQ=255)

                     aTGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATcaca                                                                                                                                 <  1:1920024/150‑1 (MQ=255)

                        cccTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATg                                                                                                                              >  1:229870/1‑150 (MQ=255)

                          cTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATggg                                                                                                                            <  2:2322982/150‑1 (MQ=255)

                           tGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGt                                                                                                                           <  1:1585525/150‑1 (MQ=255)

                            ggATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGtt                                                                                                                          <  2:1233042/150‑1 (MQ=255)

                               tACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTAc                                                                                                                       <  2:966315/150‑1 (MQ=255)

                                         gTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAgcgtg                                                                                                             >  2:1206758/1‑150 (MQ=255)

                                         gTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAgcgtg                                                                                                             >  1:247127/1‑150 (MQ=255)

                                         gTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGACTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAgcgtg                                                                                                             >  2:1638549/1‑150 (MQ=255)

                                                    gCTTGATTTTCATCCATTCCGGCTAGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGt                                                                                                  >  1:1435147/1‑150 (MQ=255)

                                                        gATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCg                                                                                              <  2:1806944/150‑1 (MQ=255)

                                                         aTTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGa                                                                                             <  2:2458633/150‑1 (MQ=255)

                                                           tttCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAAt                                                                                           <  2:419324/150‑1 (MQ=255)

                                                            ttCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAAtt                                                                                          >  2:1818234/1‑150 (MQ=255)

                                                               aTCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACg                                                                                       >  1:2282437/1‑150 (MQ=255)

                                                                tCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTGCAGTTGCGAATTACGa                                                                                      <  2:245996/150‑1 (MQ=255)

                                                                tCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGa                                                                                      <  1:90669/150‑1 (MQ=255)

                                                                tCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGa                                                                                      >  1:1044796/1‑150 (MQ=255)

                                                                 ccATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGaa                                                                                     >  1:2231988/1‑150 (MQ=255)

                                                                  cATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAg                                                                                    <  1:286865/150‑1 (MQ=255)

                                                                  cATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAg                                                                                    <  1:216340/150‑1 (MQ=255)

                                                                   aTTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAgg                                                                                   >  2:2136508/1‑150 (MQ=255)

                                                                     tCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGaa                                                                                 <  1:802239/150‑1 (MQ=255)

                                                                      ccGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGaaa                                                                                >  1:1005844/1‑150 (MQ=255)

                                                                      ccGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGaaa                                                                                <  2:944109/150‑1 (MQ=255)

                                                                       cGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAg                                                                               <  1:385775/150‑1 (MQ=255)

                                                                        ggCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGc                                                                              >  2:1549243/1‑150 (MQ=255)

                                                                             gcgcAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTttc                                                                         <  2:1931618/150‑1 (MQ=255)

                                                                             gcgcAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTttc                                                                         <  1:1319004/150‑1 (MQ=255)

                                                                             gcgcAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTttc                                                                         >  1:134192/1‑150 (MQ=255)

                                                                             gcgcAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTttc                                                                         <  1:1947321/150‑1 (MQ=255)

                                                                              cgcAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATCGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTttca                                                                        >  1:2078286/1‑150 (MQ=255)

                                                                              cgcAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTttca                                                                        >  2:45915/1‑150 (MQ=255)

                                                                                cagcagGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTttcaat                                                                      <  2:1634631/150‑1 (MQ=255)

                                                                                cagcagGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTttcaat                                                                      >  1:25117/1‑150 (MQ=255)

                                                                                 agcagGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTttcaatt                                                                     <  2:81754/150‑1 (MQ=255)

                                                                                   cagGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTttcaattca                                                                   >  1:1957353/1‑150 (MQ=255)

                                                                                     gggCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTtcaattcaat                                                                 >  2:1554178/1‑150 (MQ=255)

                                                                                     gggCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTtcaattcaat                                                                 <  2:922168/150‑1 (MQ=255)

                                                                                     gggCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTtcaattcaat                                                                 <  2:1033742/150‑1 (MQ=255)

                                                                                      ggCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATa                                                                <  2:2282437/150‑1 (MQ=255)

                                                                                            gcgcTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGAttgt                                                          >  2:1957435/1‑150 (MQ=255)

                                                                                            gcgcTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGAttgt                                                          >  1:884247/1‑150 (MQ=255)

                                                                                             cgcTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTTCGAATTACGAAGGAAAGCTTTTCAATTCAATAGAttgtt                                                         >  1:557145/1‑150 (MQ=255)

                                                                                             cgcTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGGGAATTACGAAGGAAAGCGTTTCAATTCAATAGAttgtt                                                         >  1:283928/1‑150 (MQ=255)

                                                                                             cgcTCTGTTGCCACGGTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGGTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGAttgtt                                                         <  2:856861/150‑1 (MQ=255)

                                                                                              gcTCTGTTGCCACGTTTTTAACCTGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGAttgttg                                                        >  2:2113065/1‑150 (MQ=255)

                                                                                                 ctGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTaa                                                     <  2:457154/150‑1 (MQ=255)

                                                                                                  tGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAAt                                                    >  2:649982/1‑150 (MQ=255)

                                                                                                   gTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTttcaattc                                                                    >  1:723612/1‑133 (MQ=255)

                                                                                                   gTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAAtt                                                   <  1:1383299/150‑1 (MQ=255)

                                                                                                   gTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAAtt                                                   <  1:98989/150‑1 (MQ=255)

                                                                                                     tGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTAt                                                 >  1:1198835/1‑150 (MQ=255)

                                                                                                       ccACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCa                                               <  1:1977699/150‑1 (MQ=255)

                                                                                                                 aaCCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTg                                     <  1:1121427/150‑1 (MQ=255)

                                                                                                                 aaCCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTg                                     <  1:1674613/150‑1 (MQ=255)

                                                                                                                  aCCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGCAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGa                                    >  2:1625336/1‑150 (MQ=255)

                                                                                                                  aCCGGGATAAGGGCCATCTTATCGGCAGCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGa                                    <  2:2182494/150‑1 (MQ=255)

                                                                                                                    cGGGATAAGGGCCATCTTATAGGCATAGCGGTATTTAACACCGCGTTCCATCACAAAAGGTTTAATCATAGAGTGCGTGTTCCACTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATAAACTATTTTTGAAt                                  >  2:1402612/1‑150 (MQ=255)

                                                                                                                        aTAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAAc                              <  1:2424178/150‑1 (MQ=255)

                                                                                                                         tAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACg                             >  2:2042435/1‑150 (MQ=255)

                                                                                                                           aGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGAc                           <  2:1985226/150‑1 (MQ=255)

                                                                                                                            gggCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACt                          <  2:884247/150‑1 (MQ=255)

                                                                                                                             ggCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTg                         >  1:548024/1‑150 (MQ=255)

                                                                                                                              gCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTgg                        <  2:1824237/150‑1 (MQ=255)

                                                                                                                               ccATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGc                       <  1:1638549/150‑1 (MQ=255)

                                                                                                                               ccATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGc                       <  2:229870/150‑1 (MQ=255)

                                                                                                                               ccATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGc                       <  1:1647007/150‑1 (MQ=255)

                                                                                                                               ccATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGc                       >  1:1277295/1‑150 (MQ=255)

                                                                                                                                 aTCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAg                     <  1:436647/150‑1 (MQ=255)

                                                                                                                                  tCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGt                    <  1:1282145/150‑1 (MQ=255)

                                                                                                                                   cTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGta                   <  2:2231988/150‑1 (MQ=255)

                                                                                                                                    ttATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGtat                  >  2:704327/1‑150 (MQ=255)

                                                                                                                                      aTCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTatat                <  2:1001750/150‑1 (MQ=255)

                                                                                                                                        atggaacgcgGTGTTAAATACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTATATCa              >  2:227602/7‑150 (MQ=255)

                                                                                                                                         tggaacgcgGTGTTAAATACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTATATCAt             >  1:361714/6‑150 (MQ=255)

                                                                                                                                          ggaacgcgGTGTTAAATACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTATATCAtt            >  2:2076796/5‑150 (MQ=255)

                                                                                                                                          ggaacgcgGTGTTAAATACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTATATCAtt            >  1:926839/5‑150 (MQ=255)

                                                                                                                                           gaacgcgGTGTTAAATACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTATATCATTg           >  2:386909/4‑150 (MQ=255)

                                                                                                                                            aTCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTATATCATTGa          <  1:1576950/150‑1 (MQ=255)

                                                                                                                                            aTCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTATATCATTGa          <  1:1819498/150‑1 (MQ=255)

                                                                                                                                              cgcgGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTATATCATTGaaa        >  1:639817/1‑150 (MQ=255)

                                                                                                                                               gcgGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTATATCATTGAAAc       >  2:358317/1‑150 (MQ=255)

                                                                                                                                                cgGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTATATCATTGAAACg      <  1:1206758/150‑1 (MQ=255)

                                                                                                                                                 ggTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTATATCATTGAAACgg     <  2:896656/150‑1 (MQ=255)

                                                                                                                                                    aTTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTATATCATTGAAACGGAcc  <  2:134192/150‑1 (MQ=255)

                                                                                                                                                    aTTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTATATCATTGAAACGGAcc  <  2:1957353/150‑1 (MQ=255)

TTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTGATTTTCATCCATTCCGGCTTGCGCAGCAGGGCTTCGCGCTCTGTTGCCACGTTTTTAACCGGGATAAGGGCCATCTTATCGGCATCGCGGTATTTAACACCGCGTTCCATCACAATGGGTTTACTCATAGCGTGCGTGTTCCAGTTGCGAATTACGAAGGAAAGCGTTTCAATTCAATAGATTGTTGTAATTATCAACTATTTTTGAATTAACGACTGGCAGTATATCATTGAAACGGACC  >  NC_000913/660032‑660329

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: `‑`