BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,382,003	A→G	6.3%	A19A (GCA→GCG)	ompG →	outer membrane porin G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,382,003	0	A	G	6.3%	357.4 / 9.0	109	A19A (GCA→GCG)	ompG	outer membrane porin G
Reads supporting (aligned to +/- strand): ref base A (54/48); new base G (4/3); total (58/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 111 total aligned reads displayed.

ATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCTGG  >  NC_000913/1381860‑1382151

aTGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAgg                                                                                                                                                  <  3:341113/148‑1 (MQ=255)

   tCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCg                                                                                                                                               >  3:380184/1‑148 (MQ=255)

     aTTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGagg                                                                                                                                            >  2:495857/1‑149 (MQ=255)

      tttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGagg                                                                                                                                            >  3:431853/1‑148 (MQ=255)

      tttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGagg                                                                                                                                            >  4:1031210/1‑148 (MQ=255)

      tttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGagg                                                                                                                                            >  4:11039/1‑148 (MQ=255)

       ttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGccattccc                                                                                                                                          <  1:773379/149‑7 (MQ=255)

        tCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                         <  1:1004373/149‑1 (MQ=255)

          tttGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaagg                                                                                                                                       >  2:1014373/1‑149 (MQ=255)

            tGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaaggaa                                                                                                                                     <  1:942570/149‑1 (MQ=255)

            tGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaagga                                                                                                                                      <  4:425890/148‑1 (MQ=255)

             gATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAAc                                                                                                                                    <  2:1052578/149‑1 (MQ=255)

                gCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGAc                                                                                                                                 >  1:676478/1‑149 (MQ=255)

                 cAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGAc                                                                                                                                 >  3:753983/1‑148 (MQ=255)

                   gAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTg                                                                                                                               >  4:618703/1‑148 (MQ=255)

                    aaaCGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGc                                                                                                                             >  1:399462/1‑149 (MQ=255)

                      aCGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCAc                                                                                                                           <  2:1067491/149‑1 (MQ=255)

                       cGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaagga                                                                                                                                      >  4:478692/1‑137 (MQ=255)

                       cGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaagga                                                                                                                                      <  3:478692/137‑1 (MQ=255)

                       cGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCAc                                                                                                                           >  3:784951/1‑148 (MQ=255)

                        ggAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGGCTGGCACtt                                                                                                                         >  1:891867/1‑149 (MQ=255)

                          aaaTAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACttt                                                                                                                        <  3:396408/148‑1 (MQ=255)

                           aaTAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTaa                                                                                                                      <  1:472828/149‑1 (MQ=255)

                               gCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAatat                                                                                                                   <  3:142899/148‑1 (MQ=255)

                                cAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATc                                                                                                                  <  4:798790/148‑1 (MQ=255)

                                    tCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTTTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGgcgc                                                                                                             <  2:264213/149‑1 (MQ=255)

                                    tCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGAGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGgcg                                                                                                              >  3:407063/1‑148 (MQ=255)

                                         aaaTACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATg                                                                                                         >  4:374384/1‑148 (MQ=255)

                                             cctgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACg                                                                                                     >  4:732273/4‑148 (MQ=255)

                                                 ggggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGgcg                                                                                                              <  3:392624/135‑1 (MQ=255)

                                                 ggggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGgcg                                                                                                              >  4:392624/1‑135 (MQ=255)

                                                 ggggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCAGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAAt                                                                                                 <  3:749379/148‑1 (MQ=255)

                                                  gggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATCACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATa                                                                                                >  4:805599/1‑148 (MQ=255)

                                                   ggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGa                                                                                              >  2:452640/1‑149 (MQ=255)

                                                   ggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGa                                                                                              >  2:323712/1‑149 (MQ=255)

                                                    gAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGaa                                                                                             >  2:941494/1‑149 (MQ=255)

                                                     aaGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGaa                                                                                             >  3:343221/1‑148 (MQ=255)

                                                       ggCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAAc                                                                                          <  1:519572/149‑1 (MQ=255)

                                                       ggCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAAc                                                                                          >  2:966821/1‑149 (MQ=255)

                                                        gCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAAc                                                                                          >  3:944205/1‑148 (MQ=255)

                                                          aTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTc                                                                                       >  2:993974/1‑149 (MQ=255)

                                                              ccccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGAGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAg                                                                                    <  4:407063/148‑1 (MQ=255)

                                                               cccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGa                                                                                                                                  >  3:817860/1‑101 (MQ=255)

                                                               cccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGa                                                                                                                                  <  4:817860/101‑1 (MQ=255)

                                                                ccAGGATAATACAAGGATCAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTCAATATCGGCGCGATGTACGAAATAGAAAACGTCGAggg                                                                                  >  4:82398/1‑148 (MQ=255)

                                                                ccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAggg                                                                                  >  3:569764/1‑148 (MQ=255)

                                                                 cAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGt                                                                                 >  3:384388/1‑148 (MQ=255)

                                                                   ggATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTa                                                                               <  3:432232/148‑1 (MQ=255)

                                                                    gATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATg                                                                             >  1:582473/1‑149 (MQ=255)

                                                                     ataataCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATg                                                                             <  4:343221/148‑1 (MQ=255)

                                                                         taCAAGGAACAGTAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGa                                                                         <  3:686273/148‑1 (MQ=255)

                                                                         taCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGa                                                                         <  3:11039/148‑1 (MQ=255)

                                                                          aCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGaa                                                                        >  4:134878/1‑148 (MQ=255)

                                                                           cAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGa                                                                      <  1:658075/149‑1 (MQ=255)

                                                                            aaGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGat                                                                     <  1:169653/149‑1 (MQ=255)

                                                                              ggAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGatat                                                                   >  1:1067144/1‑149 (MQ=255)

                                                                                aaCAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGatat                                                                   <  3:448137/147‑1 (MQ=255)

                                                                                aaCAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGatat                                                                   >  4:448137/1‑147 (MQ=255)

                                                                                 aCAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatgga                                                                >  2:171173/1‑149 (MQ=255)

                                                                                 aCAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatgg                                                                 <  4:784951/148‑1 (MQ=255)

                                                                                      aatGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGTGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATggg                                                            <  3:161098/148‑1 (MQ=255)

                                                                                         gAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTgg                                                        <  1:891633/149‑1 (MQ=255)

                                                                                          aaaaaGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTgg                                                        >  4:424260/1‑148 (MQ=255)

                                                                                          aaaaaGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTgg                                                        >  3:336835/1‑148 (MQ=255)

                                                                                           aaaaGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCg                                                      >  1:58973/1‑149 (MQ=255)

                                                                                             aaGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCgg                                                     <  4:569764/148‑1 (MQ=255)

                                                                                              aGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGa                                                    <  3:857450/148‑1 (MQ=255)

                                                                                                  attaCCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCtt                                               >  1:812449/1‑149 (MQ=255)

                                                                                                    taCCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGAAGCCTTc                                              <  3:974985/148‑1 (MQ=255)

                                                                                                         tGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTa                                        >  1:322061/1‑149 (MQ=255)

                                                                                                          gTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTa                                        >  3:185509/1‑148 (MQ=255)

                                                                                                          gTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTAt                                       <  1:323712/149‑1 (MQ=255)

                                                                                                              cGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTa                                   >  1:888605/1‑149 (MQ=255)

                                                                                                               gCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTaa                                  <  2:58973/149‑1 (MQ=255)

                                                                                                               gCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTa                                   <  4:494544/148‑1 (MQ=255)

                                                                                                                    ggTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATgcc                              >  3:1046032/1‑148 (MQ=255)

                                                                                                                      tGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATgccgcc                           >  2:314799/1‑149 (MQ=255)

                                                                                                                      tGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATgccgc                            >  4:90071/1‑148 (MQ=255)

                                                                                                                         tgtgtgCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAAt                                                                                                 >  3:742951/1‑76 (MQ=255)

                                                                                                                         tgtgtgCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAAt                                                                                                 <  4:742951/76‑1 (MQ=255)

                                                                                                                              gCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTgg                                                        <  3:529861/112‑1 (MQ=255)

                                                                                                                              gCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTgg                                                        >  4:529861/1‑112 (MQ=255)

                                                                                                                               tttcctcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGCTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCg                  >  1:588097/8‑149 (MQ=255)

                                                                                                                                gggAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGt                 <  1:141284/149‑1 (MQ=255)

                                                                                                                                  gAATGGCCTGCGCACAGGCCGAGGAAAGGTACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCg                                                                          >  1:854460/1‑90 (MQ=255)

                                                                                                                                  gAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCg                                                                          <  2:854460/90‑1 (MQ=255)

                                                                                                                                   ctcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatggatg                                                              >  1:25853/4‑101 (MQ=255)

                                                                                                                                   ctcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatggatg                                                              <  2:25853/98‑1 (MQ=255)

                                                                                                                                   ctcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGga              >  1:674039/4‑149 (MQ=255)

                                                                                                                                   aaTGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGga              <  2:716771/149‑1 (MQ=255)

                                                                                                                                    tcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatggatgg                                                             >  1:930229/3‑101 (MQ=255)

                                                                                                                                    tcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatggatgg                                                             <  2:930229/99‑1 (MQ=255)

                                                                                                                                     tGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCt                                                          <  1:147814/103‑1 (MQ=255)

                                                                                                                                     tGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCt                                                          >  2:147814/1‑103 (MQ=255)

                                                                                                                                      ggCCTGCGCACAGGCCGAGGAAAGGAACGGCTGGCACTTTAATATCGGCGCGATGTACGAAATAGACAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCCGAGCCTTCAGTCTATTTTAATGCCGCAAACGGGCCGTGGAGaa           <  2:891867/149‑1 (MQ=255)

                                                                                                                                      ggCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGgaga            <  3:90071/148‑1 (MQ=255)

                                                                                                                                       gCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAAt          <  1:887886/149‑1 (MQ=255)

                                                                                                                                          tGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTg        >  4:946530/1‑148 (MQ=255)

                                                                                                                                            cgcACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGctc     <  1:452640/149‑1 (MQ=255)

                                                                                                                                               acaGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCTgg  <  1:371969/149‑1 (MQ=255)

ATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCTGG  >  NC_000913/1381860‑1382151

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