BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	967,268	A→G	E216E (GAA→GAG)	msbA →	ATP‑binding lipopolysaccharide transport protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	967,268	0	A	G	100.0%	1180.5 / NA	284	E216E (GAA→GAG)	msbA	ATP‑binding lipopolysaccharide transport protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (155/129); total (155/129)

Only 100 of 291 total aligned reads displayed.

TTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAAGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGAAAATGGTTTCAGCCTCTT  >  NC_000913/967171‑967365

 tGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTAt                                                                                               <  1:1636847/101‑1 (MQ=255)

  gTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTAtt                                                                                              <  2:1636836/101‑1 (MQ=255)

   tATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTg                                                                                             >  1:1637191/1‑101 (MQ=255)

     tCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGAt                                                                                           >  1:1637195/1‑101 (MQ=255)

     tCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGAt                                                                                           >  2:1637194/1‑101 (MQ=255)

      cGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGAtt                                                                                          <  1:1636990/101‑1 (MQ=255)

      cGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGAtt                                                                                          <  1:1636824/101‑1 (MQ=255)

       gAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGAttt                                                                                         >  2:1637199/1‑101 (MQ=255)

       gAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGAttt                                                                                         <  1:1636865/101‑1 (MQ=255)

        aaGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGAtttt                                                                                        >  2:1637201/1‑101 (MQ=255)

         aGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTc                                                                                       <  1:1636769/101‑1 (MQ=255)

         aGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTc                                                                                       >  1:1637203/1‑101 (MQ=255)

          gCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGAttt                                                                                         <  1:1636677/98‑1 (MQ=255)

            gTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCggt                                                                                    <  1:1636815/101‑1 (MQ=255)

            gTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCgg                                                                                     >  2:1637208/1‑100 (MQ=255)

             ttttCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGAt                                                                                           <  2:1636743/93‑1 (MQ=255)

             ttttCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCggtg                                                                                   <  1:1636544/101‑1 (MQ=255)

             ttttCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCggtg                                                                                   >  1:1637211/1‑101 (MQ=255)

             ttttCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCggtg                                                                                   >  1:1637214/1‑101 (MQ=255)

             ttttCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCggtg                                                                                   <  1:1636828/101‑1 (MQ=255)

              tttCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCggtgg                                                                                  >  1:1637219/1‑101 (MQ=255)

              tttCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCggtgg                                                                                  >  2:1637218/1‑101 (MQ=255)

                tCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTc                                                                                <  2:1636608/101‑1 (MQ=255)

                  gCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAg                                                                              >  1:1637225/1‑101 (MQ=255)

                  gCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGAATTGATTTTCGGTGGTCAg                                                                              >  2:1637224/1‑101 (MQ=255)

                   cAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAgg                                                                             >  1:1637227/1‑101 (MQ=255)

                    aaCATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGa                                                                            <  2:1636771/101‑1 (MQ=255)

                     aCATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGaa                                                                           <  2:1636944/101‑1 (MQ=255)

                     aCATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGaa                                                                           >  1:1637233/1‑101 (MQ=255)

                     aCATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGaa                                                                           <  1:1636787/101‑1 (MQ=255)

                      cATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAg                                                                          >  2:1637234/1‑101 (MQ=255)

                       aTCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGt                                                                         <  2:1636619/101‑1 (MQ=255)

                        tCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTg                                                                        >  1:1637237/1‑101 (MQ=255)

                         cAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTgg                                                                       >  2:1637241/1‑101 (MQ=255)

                           gTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGgaa                                                                     >  1:1637244/1‑101 (MQ=255)

                            tAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGgaaa                                                                    >  2:1637247/1‑101 (MQ=255)

                             aaaaaCATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGgaaac                                                                   <  1:1637035/101‑1 (MQ=255)

                             aaaaaCATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGgaaac                                                                   <  1:1636892/101‑1 (MQ=255)

                               aaaCATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGgaaacga                                                                 >  1:1637251/1‑101 (MQ=255)

                               aaaCATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGgaaacga                                                                 >  2:1637252/1‑101 (MQ=255)

                                aaCATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGgaaacgaa                                                                >  1:1637255/1‑101 (MQ=255)

                                 aCATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGgaaacgaaa                                                               >  2:1637258/1‑101 (MQ=255)

                                  cATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGgaaacgaaac                                                              >  1:1637261/1‑101 (MQ=255)

                                    tGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGc                                                            >  1:1637265/1‑101 (MQ=255)

                                      cAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCtt                                                          >  2:1637267/1‑101 (MQ=255)

                                      cAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCtt                                                          >  2:1637270/1‑101 (MQ=255)

                                      cAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCtt                                                          <  2:1636838/101‑1 (MQ=255)

                                       aGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCttt                                                         >  1:1637273/1‑101 (MQ=255)

                                        gAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTg                                                        <  1:1636859/101‑1 (MQ=255)

                                         aaCACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGa                                                       >  1:1637276/1‑101 (MQ=255)

                                          acacCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGAt                                                      <  1:1636702/101‑1 (MQ=255)

                                           cacCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATa                                                     >  1:1637283/1‑101 (MQ=255)

                                           cacCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATa                                                     >  2:1637281/1‑101 (MQ=255)

                                           cacCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATa                                                     >  1:1637280/1‑101 (MQ=255)

                                           cacCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATa                                                     <  1:1636798/101‑1 (MQ=255)

                                              cATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAg                                                  >  1:1637288/1‑101 (MQ=255)

                                              cATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAg                                                  <  2:1636827/101‑1 (MQ=255)

                                                  gggCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCa                                               <  2:1636880/100‑1 (MQ=255)

                                                   ggCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGc                                             >  2:1637292/1‑101 (MQ=255)

                                                   ggCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGc                                             <  1:1636763/101‑1 (MQ=255)

                                                    gCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCa                                            >  1:1637293/1‑101 (MQ=255)

                                                      aGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAAc                                          <  1:1636942/101‑1 (MQ=255)

                                                       ggTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAAcc                                         <  2:1636911/101‑1 (MQ=255)

                                                       ggTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAAcc                                         >  1:1637298/1‑101 (MQ=255)

                                                        gTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCg                                        >  2:1637301/1‑101 (MQ=255)

                                                          gACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACTGaa                                      <  1:1636971/101‑1 (MQ=255)

                                                          gACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGaa                                      >  1:1637302/1‑101 (MQ=255)

                                                           accaccAGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAAt                                     <  1:1636733/101‑1 (MQ=255)

                                                             caccaGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGc                                   >  2:1637306/1‑101 (MQ=255)

                                                              accaGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCg                                  >  1:1637308/1‑101 (MQ=255)

                                                               ccaGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGt                                 <  1:1636810/101‑1 (MQ=255)

                                                                caGCGCAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTc                                >  1:1637311/1‑101 (MQ=255)

                                                                   cgcAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTc                             >  1:1637315/1‑101 (MQ=255)

                                                                    gcAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCa                            <  2:1636802/101‑1 (MQ=255)

                                                                    gcAGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCa                            <  2:1636893/101‑1 (MQ=255)

                                                                      aGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAgg                          <  2:1636964/101‑1 (MQ=255)

                                                                      aGAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAgg                          >  2:1637319/1‑101 (MQ=255)

                                                                       gAACAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAggg                         >  1:1637321/1‑101 (MQ=255)

                                                                        aaCAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAgggg                                                                       >  2:1637324/1‑54 (MQ=39)

                                                                          cAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGAt                      <  1:1636910/101‑1 (MQ=255)

                                                                          cAAATGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGAt                      >  1:1637326/1‑101 (MQ=255)

                                                                           aaaTGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATg                     >  2:1637329/1‑101 (MQ=255)

                                                                             aTGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGaa                                      >  2:1637333/1‑82 (MQ=255)

                                                                              tGCTGAAGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGaaa                  <  2:1636928/101‑1 (MQ=255)

                                                                                   aaGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGa                                       >  1:1637337/1‑75 (MQ=255)

                                                                                    aGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGAAAATGGtt            >  2:1637338/1‑101 (MQ=255)

                                                                                    aGGGCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGAAAATGGtt            <  1:1636904/101‑1 (MQ=255)

                                                                                     gggCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGAAAATGGttt           <  2:1636549/101‑1 (MQ=255)

                                                                                      ggCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGAAAAt                <  2:1637064/95‑1 (MQ=255)

                                                                                      ggCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGAAAATGGTTTc          <  1:1636997/101‑1 (MQ=255)

                                                                                       gCCACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGAAAATGGTTTCa         <  1:1636872/101‑1 (MQ=255)

                                                                                        ccACAAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAATGCTTTGATAAAGTCAGCAACCGAAt                                     <  1:1636852/72‑1 (MQ=255)

                                                                                         cacaAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGAAAATGGTTTCAGc       >  1:1637341/1‑101 (MQ=255)

                                                                                          acaAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGAAAATGGTTTCAGcc      >  2:1637343/1‑101 (MQ=255)

                                                                                           caAAGAGGTATTGATTTTCGGTGGTCAGGAAGTGgaaacgaaac                                                              <  2:1636708/44‑1 (MQ=255)

                                                                                             aagaggTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGAAAATGGTTTCAGCctct   <  2:1636982/101‑1 (MQ=255)

                                                                                             aagaggTATTGATTTTCGGTGGTCAGGAAGTGGAAAAGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGAAAATGGTTTCa         >  2:1637348/1‑95 (MQ=255)

                                                                                              agaggTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGAAAATGGTTTCAGCCTCtt  <  1:1637103/101‑1 (MQ=255)

                                                                                              agaggTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGAAAATGGTTTCAGCCTCtt  <  1:1636235/101‑1 (MQ=255)

                                                                                              agaggTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGAAAATGGTTTCAGCCTCtt  >  1:1637350/1‑101 (MQ=255)

TTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTGAAGGGCCACAAAGAAGTATTGATTTTCGGTGGTCAGGAAGTGGAAACGAAACGCTTTGATAAAGTCAGCAACCGAATGCGTCTTCAGGGGATGAAAATGGTTTCAGCCTCTT  >  NC_000913/967171‑967365

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