BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	130,291	C→G	6.4%	G324R (GGT→CGT) ‡	yacH ←	DUF3300 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	130,291	0	C	G	6.4%	381.4 / 15.1	109	G324R (GGT→CGT) ‡	yacH	DUF3300 family protein
Reads supporting (aligned to +/- strand): ref base C (50/52); new base G (4/3); total (54/55)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-01

Only 100 of 109 total aligned reads displayed.

ATAGGGCACACCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAACGGT  >  NC_000913/130145‑130436

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                                                                                                                                         <  1:840248/148‑1 (MQ=255)

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

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

          ccATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtggt                                                                                                                                        >  1:737767/1‑148 (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)

                 ggTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATa                                                                                                                                 <  2:927874/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)

                                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)

                                        gccatATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgatgat                                                                                                          >  1:433563/1‑148 (MQ=255)

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

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

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

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

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

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

                                                 ttATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtc                                                                                                >  4:697225/1‑149 (MQ=255)

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

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

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

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

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

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

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

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

                                                       gTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtcgtcgtc                                                                                          >  3:836089/1‑149 (MQ=255)

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

                                                        tAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtcgtcgtcg                                                                                         >  4:58090/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                                                                                                                                       >  2:910041/1‑88 (MQ=255)

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

                                                                       tGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGAtgct                                                                           <  2:433563/148‑1 (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)

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

                                                                                tGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGt                                                                                                       <  2:962962/111‑1 (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)

                                                                                            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                                  >  4:820482/1‑149 (MQ=255)

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

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

                                                                                                                   ttGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGaa                              >  4:45506/1‑149 (MQ=255)

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

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

                                                                                                                         tGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCcgcg                        <  3:697225/149‑1 (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)

                                                                                                                                  ttACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCt               <  3:503475/149‑1 (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)

ATAGGGCACACCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAACGGT  >  NC_000913/130145‑130436

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