BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,325,494	C→G	A261G (GCG→GGG)	yjdA →	putative vimentin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,325,494	0	C	G	99.6%	805.7 / NA	231	A261G (GCG→GGG)	yjdA	putative vimentin
Reads supporting (aligned to +/- strand): ref base C (1/0); new base G (117/113); total (118/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 238 total aligned reads displayed.

AGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGC  >  NC_000913/4325404‑4325583

 gTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGt                                                                                        >  1:1642427/1‑93 (MQ=255)

 gTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGt                                                                                        >  1:5499908/1‑93 (MQ=255)

 gTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGt                                                                                        >  1:2861119/1‑93 (MQ=255)

  ttGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTa                                                                                       >  2:5254723/1‑93 (MQ=255)

   tGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTAc                                                                                      <  2:1904771/93‑1 (MQ=255)

    gACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACt                                                                                     <  2:6200698/93‑1 (MQ=255)

    gACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACt                                                                                     >  2:1142413/1‑93 (MQ=255)

    gACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACt                                                                                     <  1:4734783/93‑1 (MQ=255)

    gACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACt                                                                                     <  1:4301262/93‑1 (MQ=255)

     aCGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTg                                                                                    >  2:2675581/1‑93 (MQ=255)

      cGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTgg                                                                                   <  1:5946903/93‑1 (MQ=255)

       gTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGc                                                                                  >  1:1885632/1‑93 (MQ=255)

       gTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGc                                                                                  <  1:2395653/93‑1 (MQ=255)

        ttACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCg                                                                                 >  2:3880317/1‑93 (MQ=255)

          aCTGGATCCCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGt                                                                               <  1:630681/93‑1 (MQ=255)

          aCTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGt                                                                               <  1:2720612/93‑1 (MQ=255)

           cTGGATACCCCCGGGCGAAATGAAGCCGGGCAACCGCATCTGCAAAAAATTCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTg                                                                              <  1:4617549/93‑1 (MQ=255)

            tGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGc                                                                             <  2:673756/93‑1 (MQ=255)

            tGGATACCCCCGCGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGc                                                                             <  2:1484840/93‑1 (MQ=255)

             ggATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCt                                                                            <  2:1805429/93‑1 (MQ=255)

                tccccccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGa                                                                         <  1:4440500/91‑1 (MQ=255)

                tACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGa                                                                         >  1:729540/1‑93 (MQ=255)

                 aCCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGAt                                                                        >  1:3125375/1‑93 (MQ=255)

                 aCCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGCGCTgggt                                                                        >  2:1690130/1‑91 (MQ=255)

                  cccccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGAtt                                                                       >  2:3273179/1‑93 (MQ=255)

                   ccccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATta                                                                      <  2:4349679/93‑1 (MQ=255)

                    cccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATtat                                                                     <  2:3014539/93‑1 (MQ=255)

                     ccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATtata                                                                    >  2:3779413/1‑93 (MQ=255)

                     ccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATtata                                                                    >  2:1616478/1‑93 (MQ=255)

                     ccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATtata                                                                    <  1:2036530/93‑1 (MQ=255)

                     ccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGGGCGGGATtata                                                                    >  2:4780205/1‑93 (MQ=255)

                      cGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATAc                                                                   <  2:618304/93‑1 (MQ=255)

                       gggCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACg                                                                  >  2:5403658/1‑93 (MQ=255)

                        ggCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGc                                                                 >  2:6024706/1‑93 (MQ=255)

                          ccAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCaa                                                               >  2:488974/1‑93 (MQ=255)

                           cAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAAc                                                              <  1:4864585/93‑1 (MQ=255)

                           cAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAAc                                                              <  1:5601344/93‑1 (MQ=255)

                             aaTGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTg                                                            >  1:4971260/1‑93 (MQ=255)

                             aaTGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTg                                                            <  2:1212892/93‑1 (MQ=255)

                             aaTGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTg                                                            >  1:1511172/1‑93 (MQ=255)

                              aTGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGa                                                           <  1:507415/93‑1 (MQ=255)

                              aTGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGa                                                           <  1:3521957/93‑1 (MQ=255)

                              aTGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGAGCTGGATTATACGCAACTGa                                                           >  2:1836148/1‑93 (MQ=255)

                                gAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGaaa                                                         <  1:1115745/93‑1 (MQ=255)

                                 aaGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAAt                                                        <  2:394753/93‑1 (MQ=255)

                                 aaGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAAt                                                        >  2:2900306/1‑93 (MQ=255)

                                   gCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCg                                                      <  1:3225722/93‑1 (MQ=255)

                                   gCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCg                                                      >  1:4150983/1‑93 (MQ=255)

                                   gCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCg                                                      >  1:305498/1‑93 (MQ=255)

                                    ccGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGa                                                     <  1:5781800/93‑1 (MQ=255)

                                     cGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAt                                                    >  1:5267109/1‑93 (MQ=255)

                                      gggCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtc                                                   <  2:237840/93‑1 (MQ=255)

                                       gacAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtct                                                  <  1:3297686/91‑1 (MQ=255)

                                          aaCCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCg                                               <  1:6383753/93‑1 (MQ=255)

                                          aaCCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCg                                               >  1:3766846/1‑93 (MQ=255)

                                          aaCCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGGGCTGGATTATACGCAACTGAAATCGACATCCg                                               >  1:2835622/1‑93 (MQ=255)

                                           aCCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGa                                              >  2:2401093/1‑93 (MQ=255)

                                           aCCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGa                                              <  2:4995908/93‑1 (MQ=255)

                                             cGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATg                                            <  1:6118954/93‑1 (MQ=255)

                                             cGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATg                                            >  1:853980/1‑93 (MQ=255)

                                              gCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGa                                           >  2:5737254/1‑93 (MQ=255)

                                              gCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGa                                           >  2:5383377/1‑93 (MQ=255)

                                                 tCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAaga                                        >  1:1957837/1‑93 (MQ=255)

                                                 tCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAaga                                        >  1:2542153/1‑93 (MQ=255)

                                                 tCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAaga                                        <  2:6377852/93‑1 (MQ=255)

                                                 tCTGCAAAAAATGCTTAACCAGCAGCTGGAACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAaga                                        >  1:4293073/1‑93 (MQ=255)

                                                   tGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAgg                                      >  2:3711330/1‑93 (MQ=255)

                                                    gCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAggc                                     >  2:6077435/1‑92 (MQ=255)

                                                    gCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGt                                     >  1:3396576/1‑93 (MQ=255)

                                                    gCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGt                                     <  1:2100078/93‑1 (MQ=255)

                                                     cAAAAAATGCTTAGCCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTc                                    <  2:5606885/93‑1 (MQ=255)

                                                     cAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTc                                    >  1:2290723/1‑93 (MQ=255)

                                                      aaaaaaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGGAGAGGTcc                                   >  1:4917150/1‑93 (MQ=255)

                                                      aaaaaaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTcc                                   >  1:3730281/1‑93 (MQ=255)

                                                      aaaaaaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTcc                                   <  2:2578446/93‑1 (MQ=255)

                                                      aaaaaaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTcc                                   <  1:1002849/93‑1 (MQ=255)

                                                      aaaaaaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTcc                                   >  1:895217/1‑93 (MQ=255)

                                                      aaaaaaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTcc                                   <  1:4818818/93‑1 (MQ=255)

                                                       aaaaaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCg                                  >  2:354548/1‑93 (MQ=255)

                                                        aaaaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGt                                 <  2:2236866/93‑1 (MQ=255)

                                                        aaaaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGt                                 >  1:1663144/1‑93 (MQ=255)

                                                         aaaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTg                                <  1:6087508/93‑1 (MQ=255)

                                                          aaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTgg                               >  2:4646940/1‑92 (MQ=255)

                                                          aaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGa                               <  1:5737817/93‑1 (MQ=255)

                                                          aaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGa                               <  1:3915534/93‑1 (MQ=255)

                                                           aTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAg                              >  1:1180852/1‑93 (MQ=255)

                                                                tAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGCTGAAGAGGTCCGTGAGGCGAt                         >  2:5138521/1‑93 (MQ=255)

                                                                tAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGAt                         <  2:5312430/93‑1 (MQ=255)

                                                                   ccAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGAtttt                      <  2:4385078/93‑1 (MQ=255)

                                                                    cagcagCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTg                     >  2:150181/1‑93 (MQ=255)

                                                                        agcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGcag                 >  2:4907289/1‑91 (MQ=255)

                                                                         gcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGt                >  2:3312323/1‑93 (MQ=255)

                                                                          cTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTg               >  1:4019724/1‑93 (MQ=255)

                                                                          cTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTg               <  2:1537912/93‑1 (MQ=255)

                                                                           tGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTgg              <  2:2743462/93‑1 (MQ=255)

                                                                           tGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTgg              <  2:4122549/93‑1 (MQ=255)

                                                                                  cgcCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATc       <  2:4513875/93‑1 (MQ=255)

                                                                                   gcCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCg      >  2:1423805/1‑93 (MQ=255)

                                                                                     cTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGt    <  2:5807887/93‑1 (MQ=255)

                                                                                       cgggggTACTAGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGc  <  2:820875/93‑1 (MQ=255)

AGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGC  >  NC_000913/4325404‑4325583

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