BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	130,289	A→G	6.5%	G324G (GGT→GGC) ‡	yacH ←	DUF3300 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	130,289	0	A	G	6.5%	337.0 / 10.9	108	G324G (GGT→GGC) ‡	yacH	DUF3300 family protein
Reads supporting (aligned to +/- strand): ref base A (50/51); new base G (4/3); total (54/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-01

Only 100 of 110 total aligned reads displayed.

CATGATAGGGCACACCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAACGGT  >  NC_000913/130141‑130436

cATGATAGGGCACACCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGa                                                                                                                                                     <  4:120444/149‑1 (MQ=255)

    aTAGGGCACACCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCg                                                                                                                                                  >  2:627088/1‑148 (MQ=255)

    aTAGGGCACACCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGc                                                                                                                                                 >  4:579548/1‑149 (MQ=255)

        ggCACACCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCAt                                                                                                                                              <  1:417584/148‑1 (MQ=255)

         gCACACCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATc                                                                                                                                             >  1:318477/1‑148 (MQ=255)

           acacCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtg                                                                                                                                          <  3:102002/149‑1 (MQ=255)

            cacCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtgg                                                                                                                                         <  4:584224/149‑1 (MQ=255)

             acCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtgg                                                                                                                                         <  2:158653/148‑1 (MQ=255)

             acCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtgg                                                                                                                                         <  1:635870/148‑1 (MQ=255)

             acCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtgg                                                                                                                                         <  1:778975/148‑1 (MQ=255)

             acCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtgg                                                                                                                                         <  1:840248/148‑1 (MQ=255)

              ccATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtggt                                                                                                                                        >  1:555575/1‑148 (MQ=255)

               gaTTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGtga                                                                                                                                      >  4:823941/2‑149 (MQ=255)

               cATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtggtg                                                                                                                                       <  1:343753/148‑1 (MQ=255)

                aTTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGtgat                                                                                                                                     <  4:287669/149‑1 (MQ=255)

                aTTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGtgat                                                                                                                                     >  4:80812/1‑149 (MQ=255)

                aTTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGtga                                                                                                                                      <  2:694587/148‑1 (MQ=255)

                 ttACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGtgat                                                                                                                                     >  2:350043/1‑148 (MQ=255)

                  tACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGtgatga                                                                                                                                   <  4:846042/149‑1 (MQ=255)

                  tACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGtgatg                                                                                                                                    <  1:792431/148‑1 (MQ=255)

                      gTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAat                                                                                                                               >  4:328291/1‑149 (MQ=255)

                      gTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAat                                                                                                                               >  4:503475/1‑149 (MQ=255)

                          tttGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAatcatc                                                                                                                           >  3:516278/1‑149 (MQ=255)

                                  gtgCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCa                                                                                                                    <  2:652503/148‑1 (MQ=255)

                                   tgCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCat                                                                                                                   <  1:205038/148‑1 (MQ=255)

                                    gccgccATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatga                                                                                                                 <  4:6763/149‑1 (MQ=255)

                                    gccgccATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatg                                                                                                                  <  1:604134/148‑1 (MQ=255)

                                         catgccatATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATGACAatgatgatt                                                                                                            <  3:150006/149‑2 (MQ=255)

                                           tgccatATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATGACAatgatgatta                                                                                                           <  1:306463/148‑3 (MQ=255)

                                            gccatATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgatgat                                                                                                          <  1:153922/148‑1 (MQ=255)

                                             ccatATTCTTATCAGTAAGATGCTCCCCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATgg                                                                                                        >  3:602721/1‑149 (MQ=255)

                                              catATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATgg                                                                                                        >  1:631529/1‑148 (MQ=255)

                                               atatTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGGCGGTCATCgtga                                                                                                                                         <  4:906449/114‑2 (MQ=255)

                                               atatTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGGCGGTCATCgtga                                                                                                                                         >  3:906449/1‑113 (MQ=255)

                                               atatTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGt                                                                                                       >  2:83987/1‑148 (MQ=255)

                                                  ttCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATg                                                                                                   >  3:350625/1‑149 (MQ=255)

                                                   tCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATg                                                                                                   <  1:497127/148‑1 (MQ=255)

                                                     ttATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtc                                                                                                <  4:557754/149‑1 (MQ=255)

                                                      tATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtc                                                                                                >  2:353464/1‑148 (MQ=255)

                                                      tATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtc                                                                                                >  1:518295/1‑148 (MQ=255)

                                                       aTCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGGCGGTCATCgtga                                                                                                                                         >  2:496345/1‑105 (MQ=255)

                                                       aTCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGGCGGTCATCgtga                                                                                                                                         <  1:496345/106‑2 (MQ=255)

                                                        tCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtcgtc                                                                                             >  4:794418/1‑149 (MQ=255)

                                                        tCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtcgt                                                                                              <  1:570116/148‑1 (MQ=255)

                                                        tCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGACGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtcgt                                                                                              <  1:802163/148‑1 (MQ=255)

                                                            tAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtcgtcgtcg                                                                                         >  4:58090/1‑149 (MQ=255)

                                                            tAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtcgtcgtcg                                                                                         >  3:496249/1‑149 (MQ=255)

                                                              aGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATcgtcgtcgt                                                                                        <  2:631529/148‑1 (MQ=255)

                                                               gATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTcc                                                                                      <  4:215409/149‑1 (MQ=255)

                                                                aTGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTcc                                                                                      >  2:238377/1‑148 (MQ=255)

                                                                     cACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTc                                                                                 >  1:69379/1‑148 (MQ=255)

                                                                      aCCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCg                                                                                >  2:10308/1‑148 (MQ=255)

                                                                       ccGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatg                                                                                                               >  2:180546/1‑116 (MQ=255)

                                                                       ccGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatg                                                                                                               <  1:180546/116‑1 (MQ=255)

                                                                       ccGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGAt                                                                              <  4:236726/149‑1 (MQ=255)

                                                                         ggTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGAtg                                                                             <  1:238377/148‑1 (MQ=255)

                                                                           tGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtggtg                                                                                                                                       <  1:910041/88‑1 (MQ=255)

                                                                           tGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtggtg                                                                                                                                       >  2:910041/1‑88 (MQ=255)

                                                                             aTACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATgctgct                                                                        <  4:497493/149‑1 (MQ=255)

                                                                                cGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGa                                                                      <  1:764174/148‑1 (MQ=255)

                                                                                  gTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATa                                                                   >  3:571485/1‑149 (MQ=255)

                                                                                   ttGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATa                                                                   >  2:155716/1‑148 (MQ=255)

                                                                                   ttGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAg                                                                  >  3:646173/1‑149 (MQ=255)

                                                                                    tGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGt                                                                                                       <  2:962962/111‑1 (MQ=255)

                                                                                     gAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTg                                                                >  4:387121/1‑149 (MQ=255)

                                                                                      aaaTTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGc                                                               <  4:818538/149‑1 (MQ=255)

                                                                                        aTTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCg                                                              >  2:78813/1‑148 (MQ=255)

                                                                                            ttgACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACg                                                          >  2:146187/1‑148 (MQ=255)

                                                                                              gACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTgg                                                       <  4:372947/149‑1 (MQ=255)

                                                                                              gACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTgg                                                       <  3:328291/149‑1 (MQ=255)

                                                                                                cGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGt                                                      <  1:934283/148‑1 (MQ=255)

                                                                                                   cGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCa                                                  <  4:569380/149‑1 (MQ=255)

                                                                                                    gATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCaa                                                 >  3:7033/1‑149 (MQ=255)

                                                                                                     atatTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCaa                                                 >  1:338322/1‑148 (MQ=255)

                                                                                                     atatTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAAc                                                <  3:585283/149‑1 (MQ=255)

                                                                                                      tatTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgt                                                                                                                                           >  4:295982/1‑57 (MQ=255)

                                                                                                      tatTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgt                                                                                                                                           <  3:295982/57‑1 (MQ=255)

                                                                                                       atTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGc                                              >  4:409239/1‑149 (MQ=255)

                                                                                                          gATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCa                                           >  4:822761/1‑149 (MQ=255)

                                                                                                           aTGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCa                                           <  1:83987/148‑1 (MQ=255)

                                                                                                                   gCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGc                                  <  3:395998/149‑1 (MQ=255)

                                                                                                                   gCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGc                                  >  4:820482/1‑149 (MQ=255)

                                                                                                                      gTTTTGTTGCCAGCCATTACCGTCACGATGACCGCCAGCGTGGTGCTGATAATCATCATTGTCACGATGATGAGGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGGACGTGGTAGATACGCCCATGCTATAGCCg                                >  2:256942/1‑148 (MQ=255)

                                                                                                                          tgtTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATc                            <  1:155716/148‑1 (MQ=255)

                                                                                                                          tgtTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATc                            >  2:658181/1‑148 (MQ=255)

                                                                                                                            ttGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCcgc                         <  4:602721/149‑1 (MQ=255)

                                                                                                                             tGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCcgcg                        <  3:697225/149‑1 (MQ=255)

                                                                                                                             tGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCcgc                         >  1:784726/1‑148 (MQ=255)

                                                                                                                                 aaccaTTACCGTCACGATGACCGCCATCGTGGGGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGGCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACa                    >  4:922282/3‑149 (MQ=255)

                                                                                                                                  gccaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTc                                                                                 >  3:835929/1‑87 (MQ=255)

                                                                                                                                  gccaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTc                                                                                 <  4:835929/87‑1 (MQ=255)

                                                                                                                                     aTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGc                <  3:395111/149‑1 (MQ=255)

                                                                                                                                      ttACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCt               >  4:53263/1‑149 (MQ=255)

                                                                                                                                           accACGATGGCGGTCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAAc                                                >  2:798069/3‑111 (MQ=255)

                                                                                                                                           accACGATGGCGGTCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAAc                                                <  1:798069/109‑1 (MQ=255)

                                                                                                                                           accACGATGGCGGTCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTcaa          >  3:206768/3‑149 (MQ=255)

                                                                                                                                            tCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTcaa          >  2:337377/1‑148 (MQ=255)

                                                                                                                                                gATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAAc     <  3:51927/149‑1 (MQ=255)

                                                                                                                                                  tGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAACgg   <  3:579548/149‑1 (MQ=255)

                                                                                                                                                    aCCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAACGGt  <  1:627088/148‑1 (MQ=255)

CATGATAGGGCACACCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAACGGT  >  NC_000913/130141‑130436

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