BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,991,334	T→A	8.4%	I61N (ATT→AAT) ‡	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,991,334	0	T	A	8.4%	608.8 / 33.5	178	I61N (ATT→AAT) ‡	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base T (86/77); new base A (4/11); total (90/88)
Fisher's exact test for biased strand distribution p-value = 6.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.41e-01

Only 100 of 181 total aligned reads displayed.

TGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAG  >  NC_000913/3991187‑3991481

     gCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCt                                                                                                                                               <  4:398212/149‑1 (MQ=255)

      ccATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCtt                                                                                                                                              <  3:154161/149‑1 (MQ=255)

        aTAAATCAATTTCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTa                                                                                                                                            <  3:435652/149‑1 (MQ=255)

        aTAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCttt                                                                                                                                             >  2:622990/1‑148 (MQ=255)

        aTAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCttt                                                                                                                                             <  1:896600/148‑1 (MQ=255)

         tAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTa                                                                                                                                            >  1:47930/1‑148 (MQ=255)

         tAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttgg                                                                                                                                           <  4:552398/149‑3 (MQ=255)

           aaTCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                            <  3:894577/146‑2 (MQ=255)

            aTCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTacac                                                                                                                                         >  2:292565/1‑148 (MQ=255)

               aaTTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCt                                                                                                                                               <  3:188015/139‑1 (MQ=255)

               aaTTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCt                                                                                                                                               >  4:188015/1‑139 (MQ=255)

               aaTTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCt                                                                                                                                     >  4:302105/1‑149 (MQ=255)

                aTTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtg                                                                                                                                    >  3:758115/1‑149 (MQ=255)

                 ttGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtg                                                                                                                                    <  1:622990/148‑1 (MQ=255)

                         ggATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACt                                                                                                                            >  1:538512/1‑148 (MQ=255)

                          gATCGCGCGCTTGCTGCTATGGGGCCTACATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                            >  2:560251/1‑130 (MQ=255)

                          gATCGCGCGCTTGCTGCTATGGGGCCTACATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                            <  1:560251/131‑2 (MQ=255)

                            tCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAAc                                                                                                                        >  3:661929/1‑149 (MQ=255)

                                gcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                            <  4:644517/125‑2 (MQ=255)

                                  gcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACt                                                                                                                  >  3:361960/1‑149 (MQ=255)

                                       ctgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGACGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTTCCTTTACACGCCTGATGAAACTCAACGCCACTACCTg                                                                                                             >  3:587588/1‑149 (MQ=255)

                                        tgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtga                                                                                                                                   >  4:768511/1‑126 (MQ=255)

                                        tgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtga                                                                                                                                   <  3:768511/126‑1 (MQ=255)

                                          ctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGaa                                                                                                           <  1:100950/148‑1 (MQ=255)

                                              ggggCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGc                                                                                                      >  3:173845/1‑149 (MQ=255)

                                              ggggCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGc                                                                                                      <  4:829001/149‑1 (MQ=255)

                                                ggCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCGTTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCt                                                                                                     >  2:126341/1‑148 (MQ=255)

                                                   cTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGaa                                                                                                 <  3:9350/149‑1 (MQ=255)

                                                   cTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGa                                                                                                  <  1:781610/148‑1 (MQ=255)

                                                    tGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAAc                                                                                                >  4:44427/1‑149 (MQ=255)

                                                    tGCATTCCAACAGGTCTACAGTCTACTGCCGACATAGTTGGACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAAGGAGCTTGAAc                                                                                                <  4:587588/149‑1 (MQ=255)

                                                     gCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTAGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAAc                                                                                                          <  4:585626/138‑1 (MQ=255)

                                                     gCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAAc                                                                                                          >  3:585626/1‑138 (MQ=255)

                                                     gCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAAc                                                                                                <  1:63000/148‑1 (MQ=255)

                                                      cATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTAcc                                                                                                               <  3:253210/132‑1 (MQ=255)

                                                      cATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACt                                                                                               >  1:383838/1‑148 (MQ=255)

                                                         tCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTa                                                                                            <  2:498568/148‑1 (MQ=255)

                                                          ccAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTAt                                                                                           >  2:846307/1‑148 (MQ=255)

                                                          ccAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATc                                                                                          <  4:769189/149‑1 (MQ=255)

                                                            aaCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGt                                                                                        >  3:903062/1‑149 (MQ=255)

                                                             aCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGc                                                                                                                                                 <  2:197752/91‑1 (MQ=255)

                                                              cAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTg                                                                                       >  1:262168/1‑148 (MQ=255)

                                                               aGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATCGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGa                                                                                     >  4:522125/1‑149 (MQ=255)

                                                                ggTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGa                                                                                     >  1:96147/1‑148 (MQ=255)

                                                                 gTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGaa                                                                                    <  1:336966/148‑1 (MQ=255)

                                                                    tACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGaa                                                                                    <  4:320658/145‑1 (MQ=255)

                                                                     aCAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTc                                                                                <  1:587549/148‑1 (MQ=255)

                                                                       aGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGt                                                                             >  3:304107/1‑149 (MQ=255)

                                                                          cTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTAc                                                                           >  2:55786/1‑148 (MQ=255)

                                                                            aCTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAg                                                                         >  1:187074/1‑148 (MQ=255)

                                                                             cTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAgg                                                                        <  2:672818/148‑1 (MQ=255)

                                                                                ccGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATcc                                                                    <  4:472597/149‑1 (MQ=255)

                                                                                  gACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATccgc                                                                  <  4:361960/149‑1 (MQ=255)

                                                                                  gACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATccg                                                                   <  2:96147/148‑1 (MQ=255)

                                                                                   aCATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTc                                                                                <  1:648508/134‑1 (MQ=255)

                                                                                      ttgttgCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGaa                                                              >  3:297030/1‑149 (MQ=255)

                                                                                       tgttgCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGaaa                                                             <  4:94830/149‑1 (MQ=255)

                                                                                       tgttgCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGACATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGaa                                                              <  1:861926/148‑1 (MQ=255)

                                                                                            cACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTg                                                         <  1:610121/148‑1 (MQ=255)

                                                                                            cACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGa                                                        <  3:6844/149‑1 (MQ=255)

                                                                                              cTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCt                                                                                                     >  1:227454/1‑102 (MQ=255)

                                                                                              cTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCt                                                                                                     <  2:227454/102‑1 (MQ=255)

                                                                                               tATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCt                                                     >  4:385641/1‑149 (MQ=255)

                                                                                                aTCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCt                                                     >  2:404577/1‑148 (MQ=255)

                                                                                                 tCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCtt                                                    >  1:104784/1‑148 (MQ=255)

                                                                                                  cACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTcc                                                  <  3:862857/149‑1 (MQ=255)

                                                                                                    ccATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACATGCCTGATGCAAGTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTCCAGGATCCGCCGAAAGGTGAGCTTCCa                                                 <  2:652831/148‑1 (MQ=255)

                                                                                                    ccATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCaa                                                <  3:484006/149‑1 (MQ=255)

                                                                                                       tCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAAtt                                              >  2:787696/1‑148 (MQ=255)

                                                                                                        ccGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTa                                             >  1:497689/1‑148 (MQ=255)

                                                                                                          gCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                          <  3:771570/149‑1 (MQ=255)

                                                                                                           cTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTgg                                         <  4:898349/149‑1 (MQ=255)

                                                                                                           cTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                          <  2:580403/148‑1 (MQ=255)

                                                                                                           cTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                          <  2:896398/148‑1 (MQ=255)

                                                                                                            tAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTa                                                                                            <  3:109410/97‑1 (MQ=255)

                                                                                                            tAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTa                                                                                            >  4:109410/1‑97 (MQ=255)

                                                                                                               tGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGta                                     <  4:267873/149‑1 (MQ=255)

                                                                                                                  aacgTTACCATCAAGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAAc                                                                                                >  3:326622/4‑87 (MQ=255)

                                                                                                                   acgTTACCATCAAGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGa                                                                                                        >  4:757355/3‑78 (MQ=255)

                                                                                                                      ttACCATCAAGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCa                              >  3:545149/1‑149 (MQ=255)

                                                                                                                       tACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCAt                             <  3:381758/149‑1 (MQ=255)

                                                                                                                         ccTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATgg                           <  3:363432/149‑1 (MQ=255)

                                                                                                                          cTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATggg                          >  4:372110/1‑149 (MQ=255)

                                                                                                                           ttGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGgt                                        <  1:15079/134‑1 (MQ=255)

                                                                                                                           ttGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGgt                                        >  2:15079/1‑134 (MQ=255)

                                                                                                                           ttGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGgc                         >  3:163470/1‑149 (MQ=255)

                                                                                                                            tGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGgc                         >  1:739319/1‑148 (MQ=255)

                                                                                                                                ggTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCAc                    >  4:673958/1‑149 (MQ=255)

                                                                                                                                 gTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCAcc                   <  4:121205/149‑1 (MQ=255)

                                                                                                                                   aaCGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTa                                             <  1:163048/121‑1 (MQ=255)

                                                                                                                                   aaCGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCt                  >  1:657779/1‑148 (MQ=255)

                                                                                                                                     cGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAgg                                                                        <  2:768897/92‑1 (MQ=255)

                                                                                                                                      gTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGCACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTc              >  3:642768/1‑149 (MQ=255)

                                                                                                                                         cccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCg            <  2:370624/148‑1 (MQ=255)

                                                                                                                                          ccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTa          <  4:881137/149‑1 (MQ=255)

                                                                                                                                           cAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTa          <  1:773645/148‑1 (MQ=255)

                                                                                                                                            aaaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAg         >  1:699028/1‑148 (MQ=255)

                                                                                                                                             aaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAgg        >  2:297611/1‑148 (MQ=255)

                                                                                                                                                 cATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCaa    >  1:814714/1‑148 (MQ=255)

                                                                                                                                                   tttGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAg  <  2:814714/148‑1 (MQ=255)

TGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAG  >  NC_000913/3991187‑3991481

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