BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,253,830	C→G	5.3%	S708T (AGT→ACT)	ycgV ←	putative adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,253,830	0	C	G	5.3%	474.4 / 13.6	132	S708T (AGT→ACT)	ycgV	putative adhesin
Reads supporting (aligned to +/- strand): ref base C (65/60); new base G (4/3); total (69/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.13e-01

Only 100 of 133 total aligned reads displayed.

GCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAG  >  NC_000913/1253683‑1253977

  ggTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTg                                                                                                                                                   <  3:641176/148‑1 (MQ=255)

      ccGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCaaa                                                                                                                                              <  1:142772/149‑1 (MQ=255)

      ccGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCaa                                                                                                                                               <  3:176866/148‑1 (MQ=255)

       cGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAg                                                                                                                                             <  1:710540/149‑1 (MQ=255)

         tCGCCTCCGCTATAGTCCGGCGATGTGTGGGTTGAGGCAATATACAGACCGACATACAAGGGCCTTACATCAGATAGACGCTTATCCCCACCAAACTGGATACCGCTGTAACCAATGTTAAAGCAGCTCAGTTTGCCACTGGCAAAgg                                                                                                                                            <  4:656027/148‑1 (MQ=255)

         tCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCAGTGGCaaac                                                                                                                                             >  1:1015764/1‑146 (MQ=255)

         tCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCAGTGGCaaac                                                                                                                                             <  2:1015764/147‑2 (MQ=255)

          cGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAg                                                                                                                                          <  2:1073409/149‑1 (MQ=255)

             cTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTg                                                                                                                                                   <  1:326523/137‑1 (MQ=255)

             cTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTg                                                                                                                                                   >  2:326523/1‑137 (MQ=255)

             cTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACAAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTcc                                                                                                                                       >  1:591164/1‑149 (MQ=255)

              tCCGCTATGGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATGCAACGGCATTACATCAGAGAGACGTTTATCCACACCAAACTGGATACCGCTGTAACCCACGTCAAAGCCGCACAGTTTGCCACTGGCAAAGGAGTcc                                                                                                                                       <  3:657378/148‑1 (MQ=255)

              accGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTcc                                                                                                                                       <  4:376309/147‑1 (MQ=255)

                 gCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACACCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGc                                                                                                                                   <  2:1001560/149‑1 (MQ=255)

                        tCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTcccc                                                                                                                             <  4:971329/148‑1 (MQ=255)

                               aTGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACt                                                                                                                      >  4:936944/1‑148 (MQ=255)

                                tGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                    <  1:161380/149‑1 (MQ=255)

                                   aTGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAg                                                                                                                                             >  2:764638/1‑121 (MQ=255)

                                   aTGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAg                                                                                                                                             <  1:764638/121‑1 (MQ=255)

                                       gtTGAGCCCATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCa                                                                                                             >  2:568787/1‑149 (MQ=255)

                                       gtTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCa                                                                                                             >  1:450924/1‑149 (MQ=255)

                                        ttGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGTGACGTTTATCCCCACCAAACAGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTACCCCCATTACTGCGCAACCa                                                                                                             <  4:371361/148‑1 (MQ=255)

                                        ttGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAg                                                                                                            >  2:1042248/1‑149 (MQ=255)

                                          gAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGat                                                                                                          >  1:143703/1‑149 (MQ=255)

                                           aGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGata                                                                                                         <  1:1042248/149‑1 (MQ=255)

                                            gCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGatat                                                                                                        >  1:596795/1‑149 (MQ=255)

                                             ccAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGatat                                                                                                        >  4:1048036/1‑148 (MQ=255)

                                              cAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATAtt                                                                                                       >  3:145260/1‑148 (MQ=255)

                                               aaTATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATa                                                                                                                         >  4:583874/1‑129 (MQ=255)

                                               aaTATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTAc                                                                                                     >  2:37143/1‑149 (MQ=255)

                                                atatACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTAc                                                                                                     <  4:492935/148‑1 (MQ=255)

                                                       gACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTa                                                                                             >  1:437391/1‑149 (MQ=255)

                                                        aCCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTAc                                                                                            >  1:114428/1‑149 (MQ=255)

                                                            aCATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGGTACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACAGCGCAACCAGATATTACCGTCTTTACTCTg                                                                                        >  1:696573/1‑149 (MQ=255)

                                                             cATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTg                                                                                        >  3:892826/1‑148 (MQ=255)

                                                             cATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGa                                                                                       <  2:359537/149‑1 (MQ=255)

                                                             cATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGa                                                                                       <  1:284350/149‑1 (MQ=255)

                                                             cATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGa                                                                                       <  2:209295/149‑1 (MQ=255)

                                                                aCAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGAttt                                                                                    >  2:1009342/1‑149 (MQ=255)

                                                                  aaCGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCtt                                                                                               <  1:262404/136‑1 (MQ=255)

                                                                  aaCGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCtt                                                                                               >  2:262404/1‑136 (MQ=255)

                                                                  aaCGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCg                                                                                  >  2:701017/1‑149 (MQ=255)

                                                                   aCGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGc                                                                                 >  1:183545/1‑149 (MQ=255)

                                                                    cGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCa                                                                                <  2:6775/149‑1 (MQ=255)

                                                                     ggCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCa                                                                                >  3:382849/1‑148 (MQ=255)

                                                                             aTCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCa                                                                       >  2:462644/1‑149 (MQ=255)

                                                                               cAGAGAGACGTTTATCCCCACCAATCTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACCGGCAAAgg                                                                                                                                            <  3:739340/78‑1 (MQ=255)

                                                                               cAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAgg                                                                                                                                            >  4:739340/1‑78 (MQ=255)

                                                                               cAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGt                                                                             <  3:380414/141‑1 (MQ=255)

                                                                                agagagACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATcc                                                                    >  1:1033791/1‑149 (MQ=255)

                                                                                   gagaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGtt                                                                 <  1:178192/149‑1 (MQ=255)

                                                                                   gagaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGtt                                                                 <  2:540962/149‑1 (MQ=255)

                                                                                   gagaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGt                                                                  <  4:407196/148‑1 (MQ=255)

                                                                                          ttATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCa                                                           >  4:221510/1‑148 (MQ=255)

                                                                                          ttATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCa                                                           >  4:137657/1‑148 (MQ=255)

                                                                                           tATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAg                                                          <  4:378795/148‑1 (MQ=255)

                                                                                                ccaccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTAc                                                     >  3:27696/1‑148 (MQ=255)

                                                                                                ccaccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACg                                                    >  2:864363/1‑149 (MQ=255)

                                                                                                ccaccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACg                                                    >  2:6903/1‑149 (MQ=255)

                                                                                                 caccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACgg                                                   >  1:458931/1‑149 (MQ=255)

                                                                                                 caccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACg                                                    >  4:1023430/1‑148 (MQ=255)

                                                                                                 caccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACg                                                    <  3:612861/148‑1 (MQ=255)

                                                                                                   ccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCg                                                                                                   >  3:240581/1‑99 (MQ=255)

                                                                                                   ccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCAcc                                                                         >  3:479767/1‑125 (MQ=255)

                                                                                                   ccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCAcc                                                                         <  4:479767/125‑1 (MQ=255)

                                                                                                    caAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGtt                                                 <  4:27696/148‑1 (MQ=255)

                                                                                                      aaCTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTc                                              >  2:419654/1‑149 (MQ=255)

                                                                                                       aCTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTc                                              >  3:215357/1‑148 (MQ=255)

                                                                                                        cTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGtttt                                               <  4:748222/146‑1 (MQ=255)

                                                                                                        cTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCa                                             <  3:338084/148‑1 (MQ=255)

                                                                                                         tGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAAc                                           <  2:986513/149‑1 (MQ=255)

                                                                                                            aTACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACAt                                         <  3:740884/148‑1 (MQ=255)

                                                                                                              aCCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGt                                      >  2:807568/1‑149 (MQ=255)

                                                                                                               ccGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGt                                      >  4:340462/1‑148 (MQ=255)

                                                                                                                 gCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTc                                    <  4:1023221/148‑1 (MQ=255)

                                                                                                                  cTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCa                                   <  3:547282/148‑1 (MQ=255)

                                                                                                                   tGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAAt                                 >  1:540613/1‑149 (MQ=255)

                                                                                                                     tAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATaa                               <  2:187027/149‑1 (MQ=255)

                                                                                                                       aCCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGtttt                                               >  3:225764/1‑131 (MQ=255)

                                                                                                                       aCCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGtttt                                               <  4:225764/131‑1 (MQ=255)

                                                                                                                          cATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTGCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTaa                           >  4:584220/1‑148 (MQ=255)

                                                                                                                          cATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTaa                           >  3:14575/1‑148 (MQ=255)

                                                                                                                           aTGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAAcc                         >  1:337123/1‑149 (MQ=255)

                                                                                                                               cAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGAc                      >  3:620909/1‑148 (MQ=255)

                                                                                                                               cAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGAc                      >  3:535848/1‑148 (MQ=255)

                                                                                                                                 aaGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGTTTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACAt                    >  3:419793/1‑148 (MQ=255)

                                                                                                                                 aaGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACAtt                   <  2:540613/149‑1 (MQ=255)

                                                                                                                                 aaGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACAt                    >  3:138004/1‑148 (MQ=255)

                                                                                                                                   gCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTga                 <  2:122526/149‑1 (MQ=255)

                                                                                                                                   gCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTg                  >  3:806634/1‑148 (MQ=255)

                                                                                                                                        tCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATa             <  3:1023430/148‑1 (MQ=255)

                                                                                                                                         tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATaa            >  4:861075/4‑148 (MQ=255)

                                                                                                                                          tctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACAt                                         <  3:861502/116‑1 (MQ=255)

                                                                                                                                          aGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAAtt          <  2:458931/149‑1 (MQ=255)

                                                                                                                                           gTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAg                                                          >  4:1056371/1‑100 (MQ=255)

                                                                                                                                           gTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAg                                                          <  3:1056371/100‑1 (MQ=255)

                                                                                                                                            tttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGt                                                                  >  3:984206/1‑91 (MQ=255)

                                                                                                                                            tttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGt                                                                  <  4:984206/91‑1 (MQ=255)

                                                                                                                                                 cACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTgccgcc    <  4:357534/148‑1 (MQ=255)

                                                                                                                                                 cACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTgccgcc    <  4:145260/148‑1 (MQ=255)

GCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAG  >  NC_000913/1253683‑1253977

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