BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,253,830	C→G	7.1%	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	7.1%	481.0 / 16.7	138	S708T (AGT→ACT)	ycgV	putative adhesin
Reads supporting (aligned to +/- strand): ref base C (77/51); new base G (4/6); total (81/57)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

GCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAG  >  NC_000913/1253683‑1253977

gCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACt                                                                                                                                                    <  2:438454/149‑1 (MQ=255)

  ggTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTgg                                                                                                                                                  >  2:855864/1‑149 (MQ=255)

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

        gTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAgg                                                                                                                                            <  4:295204/149‑1 (MQ=255)

         tCGCCTCCGCTATAGTCCGGCGATGCATGTGTTTAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGa                                                                                                                                           <  3:63885/149‑1 (MQ=255)

         tCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGa                                                                                                                                           >  4:1005/1‑149 (MQ=255)

         tCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGa                                                                                                                                           <  1:484933/149‑1 (MQ=255)

          cGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCAGTGGCaaactga                                                                                                                                          <  2:735740/149‑5 (MQ=255)

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

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

           gCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCAGTGGCaaactgag                                                                                                                                         <  1:395193/149‑6 (MQ=255)

               ccGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACATACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACGCATGTCAAAGCCGCTCAGTTTTCCACTGGCAAAGGAGTCCAg                                                                                                                                     <  4:621193/149‑1 (MQ=255)

                        tCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTccccc                                                                                                                            >  2:77227/1‑149 (MQ=255)

                             cGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAAc                                                                                                                       >  4:444104/1‑149 (MQ=255)

                              gATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACt                                                                                                                      >  2:907248/1‑149 (MQ=255)

                                tgtaTGTGGTGAGCTAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTACCCCCACCAAACTGGATACCGCTGTAACCGATGTCAAAGCCCCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                    <  4:666034/146‑1 (MQ=255)

                                tGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgt                                                                                                                    <  4:863935/149‑2 (MQ=255)

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

                                tGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                    >  2:512485/1‑149 (MQ=255)

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

                                tGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                    >  4:696449/1‑149 (MQ=255)

                                        ttGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAg                                                                                                            >  3:326102/1‑149 (MQ=255)

                                        ttGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAg                                                                                                            >  3:128882/1‑149 (MQ=255)

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

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

                                            gCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCAGTGGCaaactga                                                                                                                                          <  2:846732/115‑5 (MQ=255)

                                            gCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCAGTGGCaaactga                                                                                                                                          >  1:846732/1‑111 (MQ=255)

                                            gCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                    >  2:987860/1‑137 (MQ=255)

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

                                             ccAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATAtt                                                                                                       >  2:98917/1‑149 (MQ=255)

                                             ccAATATAAAGACCTACATACAACGGCATTAGATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCCGCATCCCATGTCAAATCCGCTGAGTTTGCCACTGCCAATGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATAtt                                                                                                       <  3:668257/149‑1 (MQ=255)

                                               aaTATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTAc                                                                                                     <  2:1046525/149‑1 (MQ=255)

                                               aaTATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTAc                                                                                                     <  1:855864/149‑1 (MQ=255)

                                                    aCAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCa                                                                                                             <  3:272676/136‑1 (MQ=255)

                                                    aCAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCt                                                                                                >  1:779930/1‑149 (MQ=255)

                                                    aCAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCt                                                                                                >  1:28114/1‑149 (MQ=255)

                                                      aGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCttt                                                                                              >  3:727866/1‑149 (MQ=255)

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

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

                                                          cGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTActc                                                                                          <  1:18562/149‑1 (MQ=255)

                                                          cGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTActc                                                                                          >  3:467352/1‑149 (MQ=255)

                                                                  aaCGGCATTACATCGGAGAGACGTTTATCCCCACCAAACTGGATACCGCCGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAACGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATGTTACCGTCTTTACTCTGATTTCg                                                                                  >  4:613253/1‑149 (MQ=255)

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

                                                                       cATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGt                                                                             <  3:283677/149‑1 (MQ=255)

                                                                           aCATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCAcc                                                                         >  2:526633/1‑149 (MQ=255)

                                                                            cATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCGATGACAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGACTTCGCGGGTCAccc                                                                        >  2:217920/1‑149 (MQ=255)

                                                                              tCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAgg                                                                              <  4:770509/141‑1 (MQ=255)

                                                                              tCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAgg                                                                              >  3:770509/1‑141 (MQ=255)

                                                                              tCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCAt                                                                      >  1:522875/1‑149 (MQ=255)

                                                                                atagagaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACCGGCAAAGGAGTCCAGGCTTCCCCCGTAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATcc                                                                    >  2:93478/3‑149 (MQ=255)

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

                                                                                    agaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTg                                                                >  2:190119/1‑149 (MQ=255)

                                                                                      aCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCa                                                              >  2:132679/1‑149 (MQ=255)

                                                                                      aCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCa                                                              >  3:387507/1‑149 (MQ=255)

                                                                                        gTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATc                                                            >  2:916351/1‑149 (MQ=255)

                                                                                         tttATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCa                                                                                >  2:877206/1‑128 (MQ=255)

                                                                                         tttATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCa                                                           >  1:500596/1‑149 (MQ=255)

                                                                                         tttATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCa                                                           <  1:403995/149‑1 (MQ=255)

                                                                                            aTCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCg                                                        >  1:183983/1‑149 (MQ=255)

                                                                                              ccccACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTa                                                      >  1:839154/1‑149 (MQ=255)

                                                                                              ccccACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTa                                                      >  4:258265/1‑149 (MQ=255)

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

                                                                                                 caccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTa                                                                                                      >  3:294292/1‑98 (MQ=255)

                                                                                                 caccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTa                                                                                                      <  4:294292/98‑1 (MQ=255)

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

                                                                                                  accaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCa                                                                                <  3:30257/119‑1 (MQ=255)

                                                                                                  accaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTg                                                                >  1:52401/1‑135 (MQ=255)

                                                                                                  accaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTg                                                                <  2:52401/135‑1 (MQ=255)

                                                                                                  accaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGt                                                  >  1:102298/1‑149 (MQ=255)

                                                                                                     aaaCTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGtttt                                               >  4:698421/1‑149 (MQ=255)

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

                                                                                                        cTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCaa                                            <  1:49803/149‑1 (MQ=255)

                                                                                                        cTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCaa                                            >  4:257954/1‑149 (MQ=255)

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

                                                                                                         tGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAAc                                           <  4:326102/149‑1 (MQ=255)

                                                                                                            aTACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATa                                        <  3:755439/149‑1 (MQ=255)

                                                                                                               ccGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGtt                                     <  4:128882/149‑1 (MQ=255)

                                                                                                                 gCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCa                                   >  2:45444/1‑149 (MQ=255)

                                                                                                                   tGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGtt                                                                 >  4:333495/1‑117 (MQ=255)

                                                                                                                   tGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAAt                                 <  2:827419/149‑1 (MQ=255)

                                                                                                                       aCCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATaa                               >  3:965937/1‑147 (MQ=255)

                                                                                                                       aCCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATaa                               <  4:965937/147‑1 (MQ=255)

                                                                                                                         ccATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTc                                                                            <  3:400943/100‑1 (MQ=255)

                                                                                                                         ccATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTc                                                                            >  4:400943/1‑100 (MQ=255)

                                                                                                                         ccATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTaa                           >  2:69656/1‑149 (MQ=255)

                                                                                                                         ccATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTaa                           >  1:815672/1‑149 (MQ=255)

                                                                                                                              tCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGAc                      >  4:496976/1‑149 (MQ=255)

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

                                                                                                                                  aGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACAt                    <  3:886001/147‑1 (MQ=255)

                                                                                                                                  aGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACAt                    >  4:886001/1‑147 (MQ=255)

                                                                                                                                    ccGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGCCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTgag                <  1:98917/149‑1 (MQ=255)

                                                                                                                                      gCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGAt              >  2:291576/1‑149 (MQ=255)

                                                                                                                                          cctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGc                                                         >  1:76455/3‑102 (MQ=255)

                                                                                                                                           ctttGCCTGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCaa                                  <  4:502822/123‑1 (MQ=255)

                                                                                                                                           ctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCaa                                  >  3:502822/2‑124 (MQ=255)

                                                                                                                                                 cACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCa   <  1:69656/149‑1 (MQ=255)

                                                                                                                                                 cACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCa   >  3:17994/1‑149 (MQ=255)

                                                                                                                                                  aCTGGCAAAGGAGTCCGGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAg  >  2:663889/1‑149 (MQ=255)

                                                                                                                                                  aCTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAg  >  2:703772/1‑149 (MQ=255)

                                                                                                                                                   cTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTg         <  4:288608/141‑1 (MQ=255)

GCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAG  >  NC_000913/1253683‑1253977

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