BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	130,293	C→T	6.5%	G323D (GGC→GAC)	yacH ←	DUF3300 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	130,293	0	C	T	6.5%	342.7 / 10.5	108	G323D (GGC→GAC)	yacH	DUF3300 family protein
Reads supporting (aligned to +/- strand): ref base C (49/52); new base T (4/3); total (53/55)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

Only 100 of 111 total aligned reads displayed.

ATAGGGCACACCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAACGGTTCTC  >  NC_000913/130145‑130440

aTAGGGCACACCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGc                                                                                                                                                     >  4:579548/1‑149 (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:778975/148‑1 (MQ=255)

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

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

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

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

                      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)

                                        gccatATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgatgat                                                                                                              >  1:433563/1‑148 (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:697225/1‑149 (MQ=255)

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

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

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

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

                                                   aTCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGGCGGTCATCgtga                                                                                                                                             <  1:496345/106‑2 (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)

                                                                 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)

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

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

                                                                                tGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGt                                                                                                           >  1:962962/1‑111 (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                                                           <  3:328291/149‑1 (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                                                                                                                                               <  3:295982/57‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                                                                  ttACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCt                   <  3:503475/149‑1 (MQ=255)

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

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

                                                                                                                                       accACGATGGCGGTCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAAc                                                    >  2:798069/3‑111 (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)

                                                                                                                                                   gCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAACGGTtct   >  1:12469/1‑148 (MQ=255)

                                                                                                                                                    ccATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAACGGTtctc  >  1:769080/1‑148 (MQ=255)

                                                                                                                                                    ccATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAACGGTtctc  <  2:318477/148‑1 (MQ=255)

ATAGGGCACACCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAACGGTTCTC  >  NC_000913/130145‑130440

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