BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	867,637	G→C	A288P (GCC→CCC)	b0829 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	867,637	0	G	C	96.4%	634.8 / 11.2	194	A288P (GCC→CCC)	b0829	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (6/1); new base C (104/83); total (110/84)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

Only 100 of 211 total aligned reads displayed.

TTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑GCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGT  >  NC_000913/867547‑867726

  aGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑Ccccc                                                                                       <  2:1280471/93‑1 (MQ=255)

  aGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑Ccccc                                                                                       >  1:3497063/1‑93 (MQ=255)

  aGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑Ccccc                                                                                       >  1:3476203/1‑93 (MQ=255)

   gctgctGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑cccccc                                                                                      >  2:2574902/1‑93 (MQ=255)

     tgctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑cccccccc                                                                                    <  1:773616/93‑1 (MQ=255)

      gctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑cccccccca                                                                                   >  1:5493948/1‑92 (MQ=255)

      gctgTTCCCCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑Ccccccccc                                                                                   <  2:2262207/93‑1 (MQ=255)

         gTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACTG‑CCCCCCCCCATc                                                                                <  2:4084324/93‑1 (MQ=255)

         gTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATc                                                                                >  2:3832924/1‑93 (MQ=255)

           tCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGa                                                                              <  1:1522028/93‑1 (MQ=255)

             cGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGc                                                                            >  1:1035390/1‑93 (MQ=255)

             cGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGc                                                                            >  2:4695423/1‑93 (MQ=255)

               cAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCCTCGAGCaa                                                                          >  1:3128112/1‑92 (MQ=255)

                 aCTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGaa                                                                        >  2:6019999/1‑93 (MQ=255)

                 aCTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACA‑‑GCCCCCCCCATCGAGCAGaaa                                                                       >  2:544150/1‑93 (MQ=255)

                  cTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACA‑‑GCCCCCCCCATCGAGCAGaaaa                                                                      >  2:863276/1‑93 (MQ=255)

                   ttGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGaaaa                                                                      >  1:4618706/1‑93 (MQ=255)

                      gTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCCTCGAGCAGAAAACgg                                                                   >  1:82226/1‑93 (MQ=255)

                      gTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACgg                                                                   >  1:6245066/1‑93 (MQ=255)

                      gTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACgg                                                                   >  2:6072753/1‑93 (MQ=255)

                       ggCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACggt                                                                  <  1:1158465/92‑1 (MQ=255)

                        gCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACggtg                                                                 >  2:1810289/1‑93 (MQ=255)

                        gCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACggtg                                                                 <  2:6224546/93‑1 (MQ=255)

                          gATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACggtggt                                                               <  2:64076/93‑1 (MQ=255)

                          gATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACggtggg                                                               >  1:6406264/1‑92 (MQ=255)

                          gATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACggtggg                                                               >  1:13250/1‑92 (MQ=255)

                           aTGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGgtggtg                                                              >  2:5043838/1‑93 (MQ=255)

                             gAAAGGGTTAGATTATCCCCGACTTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGa                                                            <  2:635257/93‑1 (MQ=255)

                             gAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGa                                                            >  1:19057/1‑93 (MQ=255)

                             gAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGa                                                            >  2:5254449/1‑93 (MQ=255)

                              aaaGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGCCCAGAAAACGGtggaggcg                                                           >  1:3871794/1‑91 (MQ=255)

                              aaaGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGag                                                           >  2:1400837/1‑92 (MQ=255)

                              aaaGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGAt                                                           >  2:6030471/1‑93 (MQ=255)

                               aaGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATg                                                          >  2:3876474/1‑93 (MQ=255)

                               aaGGGTAAGATTATCCCCGACCTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGGAGAAAACGGTGGTGGATg                                                          <  1:5971154/93‑1 (MQ=255)

                                   gTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGCGCAGAAAACGGTGGTGGATGGCGa                                                      >  1:1083404/1‑93 (MQ=255)

                                      aGATTATCCCCGCCGTTTCCCGTTGATATCGCGTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                   <  1:4436087/93‑1 (MQ=255)

                                       gATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCt                                                  >  2:356638/1‑93 (MQ=255)

                                       gATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCt                                                  <  1:721448/93‑1 (MQ=255)

                                         ttATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGt                                                <  1:3418558/93‑1 (MQ=255)

                                         ttATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGt                                                <  2:5565445/93‑1 (MQ=255)

                                          tATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGtt                                               <  2:1665118/93‑1 (MQ=255)

                                              cccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTAc                                           >  2:3652624/1‑93 (MQ=255)

                                               ccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAACTGTTTTACg                                          >  1:5568812/1‑93 (MQ=255)

                                                 gACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGCCCCCCCCCCATCGAGCAGAAAACGGGGGGGGGGGGGGAACCTGTTTTACGa                                         >  1:3595406/1‑93 (MQ=255)

                                                 gACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAg                                        >  2:1983645/1‑93 (MQ=255)

                                                  ccGTTTCCCGTTGATATCGCTTGAACATCCAGCCAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGt                                       <  2:2056962/92‑1 (MQ=255)

                                                  aCGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGt                                       <  1:3252086/93‑1 (MQ=255)

                                                  aCGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGt                                       <  2:1352877/93‑1 (MQ=255)

                                                    gTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGc                                     >  1:3305452/1‑93 (MQ=255)

                                                    gTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGc                                     <  2:4854280/93‑1 (MQ=255)

                                                    gTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGc                                     <  2:111341/93‑1 (MQ=255)

                                                     tttCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCg                                    >  1:2307963/1‑93 (MQ=255)

                                                       tCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTa                                  >  1:1014282/1‑93 (MQ=255)

                                                       tCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTa                                  >  1:2363224/1‑93 (MQ=255)

                                                       tCCCGTTGATATCGCTTGAACATCCAGCGAAACA‑‑GCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTaa                                 >  2:4332680/1‑93 (MQ=255)

                                                         ccGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAAt                                <  1:2084029/93‑1 (MQ=255)

                                                         ccGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAAt                                <  1:1800759/93‑1 (MQ=255)

                                                          catcgATATCTCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATc                               <  1:5353129/89‑1 (MQ=255)

                                                          cGTTGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATc                               <  2:2774489/93‑1 (MQ=255)

                                                           gTTGATATCGCTTGACCATCCAGCGAAACAC‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCt                              <  1:3709041/93‑1 (MQ=255)

                                                            ttGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCtt                             <  1:2030435/93‑1 (MQ=255)

                                                            ttGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCtt                             <  1:4259101/93‑1 (MQ=255)

                                                            ttGATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGAGGATGGGGGACCTGTTTTACGAGTGTGTCATCtt                             >  1:4658809/1‑93 (MQ=255)

                                                              gATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGt                           >  1:635195/1‑93 (MQ=255)

                                                              gATATCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGt                           >  2:2952467/1‑93 (MQ=255)

                                                              gATATCGCTTGAACATCCAACGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGt                           >  1:197010/1‑93 (MQ=255)

                                                               atctCGCTGGAACATGCAGGGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCATAATCTTGTc                          <  1:5028649/90‑1 (MQ=255)

                                                               atatCGCTTGACCATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTc                          <  2:5008465/93‑1 (MQ=255)

                                                                tatCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGtta                         >  1:2403563/1‑91 (MQ=255)

                                                                tatCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCa                         <  1:5111255/93‑1 (MQ=255)

                                                                  tCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCAcc                       <  2:4771790/93‑1 (MQ=255)

                                                                  tCGCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGAGGATGGCGAACCTGTTTTACGTGTGCGTCATCTTGTCAcc                       >  1:2573665/1‑93 (MQ=255)

                                                                    gCTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCg                     >  1:3770033/1‑93 (MQ=255)

                                                                     cTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGGGAACCTGTTTTACGAGTGCGTAATCTTTTCACCCGt                    >  2:2617358/1‑93 (MQ=255)

                                                                     cTTGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGt                    <  1:1180166/93‑1 (MQ=255)

                                                                      ttGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGtt                   <  2:5813617/93‑1 (MQ=255)

                                                                      ttGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGtt                   <  1:5245511/93‑1 (MQ=255)

                                                                      ttGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGACGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGtt                   >  2:546366/1‑93 (MQ=255)

                                                                       tGAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGttt                  <  1:4116230/93‑1 (MQ=255)

                                                                        gAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGGACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTc                 >  2:742948/1‑93 (MQ=255)

                                                                        gAACATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTc                 >  1:458645/1‑93 (MQ=255)

                                                                          aCATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTccc               <  1:2229739/93‑1 (MQ=255)

                                                                          aCATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTccc               >  2:1193490/1‑93 (MQ=255)

                                                                           cATCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCt              <  2:2383267/93‑1 (MQ=255)

                                                                            aTCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCtt             <  2:4415778/93‑1 (MQ=255)

                                                                            aTCCAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCtt             >  2:3114145/1‑93 (MQ=255)

                                                                              ccAGCGAAACAG‑CCGCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTg           <  2:3197100/93‑1 (MQ=255)

                                                                              ccAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTg           <  2:45083/93‑1 (MQ=255)

                                                                              ccAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTg           <  2:6432657/93‑1 (MQ=255)

                                                                               cAGCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTgc          <  1:2838565/93‑1 (MQ=255)

                                                                                 gCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTgcgc        <  2:5635776/93‑1 (MQ=255)

                                                                                 gCGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTgcgc        <  1:588947/93‑1 (MQ=255)

                                                                                  cGAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCa       <  2:2346981/93‑1 (MQ=255)

                                                                                   gAAACAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTTCCTTTGCGCAg      >  2:10772/1‑93 (MQ=255)

                                                                                    aaaCAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGc     >  2:3329167/1‑93 (MQ=255)

                                                                                     aaCAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg    >  2:1771317/1‑93 (MQ=255)

                                                                                     aaCAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg    >  1:4911821/1‑93 (MQ=255)

                                                                                     aaCAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg    >  2:3986306/1‑93 (MQ=255)

                                                                                       cAG‑CCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGt  <  1:112554/93‑1 (MQ=255)

TTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAG‑GCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGT  >  NC_000913/867547‑867726

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