BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,209,171	T→C	Y142H (TAT→CAT)	dnaG →	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,209,171	0	T	C	83.9%	585.1 / 105.1	261	Y142H (TAT→CAT)	dnaG	DNA primase
Reads supporting (aligned to +/- strand): ref base T (21/21); new base C (119/100); total (140/121)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

Only 100 of 264 total aligned reads displayed.

TCAGAGGCAAACGCTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACGTCCTGAAGCGGTTT  >  NC_000913/3209080‑3209263

  agagGCAAACGCTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATc                                                                                           <  2:6331975/93‑1 (MQ=255)

  agagGCAAACGCTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATc                                                                                           <  1:2666299/93‑1 (MQ=255)

    agGCAAACGCTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTg                                                                                         <  1:3033939/93‑1 (MQ=255)

    agGCAAACGCTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTg                                                                                         <  1:5849332/93‑1 (MQ=255)

    agGCAAACGCTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTg                                                                                         <  2:6632444/93‑1 (MQ=255)

      gCAAACGCTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGa                                                                                       <  2:3057311/93‑1 (MQ=255)

      gCAAACGCTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGa                                                                                       <  2:4866176/93‑1 (MQ=255)

         aaCGCTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGaaaa                                                                                    >  1:63144/1‑93 (MQ=255)

          aCGCTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGaaaaa                                                                                   >  1:5459690/1‑93 (MQ=255)

           cGCTTTATCAGTTGATGGATGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAAc                                                                                  <  2:40320/93‑1 (MQ=255)

           cGCTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAAc                                                                                  >  2:1650028/1‑93 (MQ=255)

            gCTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAAcg                                                                                 <  1:6563682/93‑1 (MQ=255)

             cTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAAcgc                                                                                <  2:3244538/93‑1 (MQ=255)

             cTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAAcgc                                                                                >  1:3081294/1‑93 (MQ=255)

                tATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGa                                                                             <  1:2619905/93‑1 (MQ=255)

                tATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGa                                                                             >  2:4215659/1‑93 (MQ=255)

                 aTCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGAt                                                                            <  1:932777/93‑1 (MQ=255)

                 aTCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGAt                                                                            >  1:3523832/1‑93 (MQ=255)

                 aTCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGAt                                                                            <  1:5545062/93‑1 (MQ=255)

                  tCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGAtt                                                                           <  1:4128796/93‑1 (MQ=255)

                  tCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGAtt                                                                           >  1:4461001/1‑93 (MQ=255)

                    aGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTaa                                                                         <  2:5424617/93‑1 (MQ=255)

                    aGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTaa                                                                         <  2:3616096/93‑1 (MQ=255)

                      ttGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGc                                                                       <  1:7384041/93‑1 (MQ=255)

                      ttGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGc                                                                       <  2:1716831/93‑1 (MQ=255)

                        gATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCa                                                                     <  2:2553562/93‑1 (MQ=255)

                          tGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACg                                                                   <  2:2045287/93‑1 (MQ=255)

                          tGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACg                                                                   >  1:1698334/1‑93 (MQ=255)

                           ggACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGa                                                                  <  1:6635573/93‑1 (MQ=255)

                            gACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAg                                                                 <  1:7232406/93‑1 (MQ=255)

                            gACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAg                                                                 >  2:3773178/1‑93 (MQ=255)

                            gACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAg                                                                 >  1:4877184/1‑93 (MQ=255)

                             aCGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCTACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAgg                                                                <  2:6819971/93‑1 (MQ=255)

                             aCGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAgg                                                                >  2:2443334/1‑93 (MQ=255)

                                gTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTa                                                             >  2:3411988/1‑93 (MQ=255)

                                gTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTa                                                             >  1:5992177/1‑93 (MQ=255)

                                   tGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTAtcg                                                          >  2:1951931/1‑93 (MQ=255)

                                   tGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTAtcg                                                          >  2:2992155/1‑93 (MQ=255)

                                   tGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTAtcg                                                          >  1:3787669/1‑93 (MQ=255)

                                    gAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTAtcgc                                                         <  1:1371781/93‑1 (MQ=255)

                                     aaTACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTAtcgct                                                        <  1:5142893/93‑1 (MQ=255)

                                     aaTACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTAtcgct                                                        <  2:513454/93‑1 (MQ=255)

                                      aTACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTAtcgctc                                                       >  1:3198155/1‑93 (MQ=255)

                                        aCGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTAtcgctcgc                                                     >  1:2868590/1‑93 (MQ=255)

                                          gTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCtt                                                   >  2:4613270/1‑93 (MQ=255)

                                          gTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCtt                                                   <  1:2155635/93‑1 (MQ=255)

                                           ttttACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCttt                                                  <  1:4049929/93‑1 (MQ=255)

                                           ttttACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCttt                                                  >  2:1032424/1‑93 (MQ=255)

                                             ttACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGc                                                <  1:6113741/93‑1 (MQ=255)

                                             ttACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGc                                                <  2:1251967/93‑1 (MQ=255)

                                               aCCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGa                                              >  1:4335517/1‑93 (MQ=255)

                                                ccAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGAt                                             >  1:2143236/1‑93 (MQ=255)

                                                 caacaaTCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGAtt                                            <  1:1954185/93‑1 (MQ=255)

                                                 caacaaTCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGAtt                                            >  2:1158622/1‑93 (MQ=255)

                                                  aacaaTCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTg                                           >  1:5569084/1‑93 (MQ=255)

                                                   acaaTCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTgg                                          >  1:5735758/1‑93 (MQ=255)

                                                   acaaTCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTgg                                          <  1:3736728/93‑1 (MQ=255)

                                                    caaTCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGt                                         <  2:3904816/93‑1 (MQ=255)

                                                      aTCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGttt                                       >  2:5730608/1‑93 (MQ=255)

                                                        cTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTg                                     >  1:1158954/1‑93 (MQ=255)

                                                        cTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTg                                     >  1:7004024/1‑93 (MQ=255)

                                                         tttACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTgc                                    <  2:3835314/93‑1 (MQ=255)

                                                         tttACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTgc                                    >  2:2654141/1‑93 (MQ=255)

                                                          ttACAACAACCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTgcg                                   >  2:2726092/1‑93 (MQ=255)

                                                              aacaacCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGcccc                               >  1:2598186/1‑93 (MQ=255)

                                                               acaacCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGccccc                              >  2:6092339/1‑93 (MQ=255)

                                                               acaacCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGccccc                              >  2:896147/1‑93 (MQ=255)

                                                               acaacCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGccccc                              >  1:143427/1‑93 (MQ=255)

                                                               acaacCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGccccc                              >  1:6836741/1‑93 (MQ=255)

                                                                cacccTGTTGCCCCGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGcccccc                             <  1:6436182/90‑1 (MQ=255)

                                                                caacCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGcccccc                             <  2:2617282/93‑1 (MQ=255)

                                                                caacCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGcccccc                             <  2:1469107/93‑1 (MQ=255)

                                                                 aacCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCg                            >  2:2988706/1‑93 (MQ=255)

                                                                 aacCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCg                            >  1:2498411/1‑93 (MQ=255)

                                                                  acCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGcccccccg                           >  2:4353969/1‑91 (MQ=255)

                                                                  acCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCgg                           >  2:291595/1‑93 (MQ=255)

                                                                  acCTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCgg                           >  2:2260332/1‑93 (MQ=255)

                                                                   ccTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGc                          <  2:7047031/93‑1 (MQ=255)

                                                                    cTGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCt                         >  1:1087874/1‑93 (MQ=255)

                                                                     tGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTg                        <  2:774862/93‑1 (MQ=255)

                                                                     tGTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTg                        >  2:4800965/1‑93 (MQ=255)

                                                                      gTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTgg                       >  2:2850994/1‑93 (MQ=255)

                                                                      gTTGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTgg                       <  2:5055038/93‑1 (MQ=255)

                                                                       ttGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTggg                      <  1:2567304/93‑1 (MQ=255)

                                                                       ttGCCACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTggg                      >  2:319521/1‑93 (MQ=255)

                                                                         gCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGAc                    >  2:5542997/1‑93 (MQ=255)

                                                                          ccACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACa                   >  2:6870327/1‑93 (MQ=255)

                                                                           cACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACaa                  <  2:2119584/93‑1 (MQ=255)

                                                                           cACGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACaa                  <  2:5997336/93‑1 (MQ=255)

                                                                            aCGTCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAAc                 >  2:2379659/1‑93 (MQ=255)

                                                                               tCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACGTc              >  2:4748900/1‑93 (MQ=255)

                                                                               tCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACGTc              >  1:2955206/1‑93 (MQ=255)

                                                                               tCTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACGTc              >  1:5344188/1‑93 (MQ=255)

                                                                                cTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACGTcc             <  1:1847557/93‑1 (MQ=255)

                                                                                cTGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACGTcc             <  1:2309722/93‑1 (MQ=255)

                                                                                 tGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACGTCCt            <  2:6533868/93‑1 (MQ=255)

                                                                                 tGCGCGCCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACGTCCt            >  2:7473355/1‑93 (MQ=255)

                                                                                  gcgcgcCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACGTCCTg           >  1:174285/1‑93 (MQ=255)

                                                                                    gcgcCAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACGTCCTGaa         <  2:58305/93‑1 (MQ=255)

                                                                                        cAGCATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACATCCTGAAGCgg     >  1:389523/1‑93 (MQ=255)

TCAGAGGCAAACGCTTTATCAGTTGATGGACGGTCTGAATACGTTTTACCAACAATCTTTACAACAACCTGTTGCCACGTCTGCGCGCCAGTATCTGGAAAAACGCGGATTAAGCCACGAGGTTATCGCTCGCTTTGCGATTGGTTTTGCGCCCCCCGGCTGGGACAACGTCCTGAAGCGGTTT  >  NC_000913/3209080‑3209263

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