BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,393,906	T→C	intergenic (‑4243/+57)	tldD ← / ← cafA	TldD/CafA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,393,906	0	T	C	99.7%	1009.0 / ‑5.3	320	intergenic (‑4243/+57)	tldD/cafA	TldD/CafA
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (145/173); total (146/174)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 324 total aligned reads displayed.

CGCAGGCCGCTAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCTTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTCCTGGTT  >  NC_000913/3393816‑3393995

cGCAGGCCGCTAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTg                                                                                         <  1:691471/93‑1 (MQ=255)

  cAGGCCGCTAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTc                                                                                       <  1:1167750/93‑1 (MQ=255)

  cAGGCCGCTAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTc                                                                                       <  1:7061855/93‑1 (MQ=255)

   aGGCCGCTAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCt                                                                                      <  1:3971467/93‑1 (MQ=255)

    ggCCGCTAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTg                                                                                     >  2:2726753/1‑93 (MQ=255)

    ggCCGCTAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCACCGGCAATCGCCTCACCCGTCACTCCCCGTCtc                                                                                     >  1:2877492/1‑92 (MQ=255)

       cGCTAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTa                                                                                  <  2:5502203/93‑1 (MQ=255)

         cTAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTaaa                                                                                >  2:6783113/1‑93 (MQ=255)

         cTAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTaaa                                                                                >  2:3726886/1‑93 (MQ=255)

           aaCCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTaaaaa                                                                              <  2:1690162/93‑1 (MQ=255)

           aaCCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTaaaaa                                                                              <  1:1861554/93‑1 (MQ=255)

            aCCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAAt                                                                             >  1:5240422/1‑93 (MQ=255)

             ccAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAAtg                                                                            <  2:1848975/93‑1 (MQ=255)

             ccAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAAtg                                                                            >  1:7720843/1‑93 (MQ=255)

               agcagGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAAtgtg                                                                          >  1:2657965/1‑93 (MQ=255)

               agcagGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAAtgtg                                                                          <  2:4238881/93‑1 (MQ=255)

               agcagGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAAtgtg                                                                          >  1:6117940/1‑93 (MQ=255)

               agcagGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAAtgtg                                                                          >  1:4348247/1‑93 (MQ=255)

                gcagGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGa                                                                         >  1:7159841/1‑93 (MQ=255)

                  agGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACt                                                                       <  2:2221400/93‑1 (MQ=255)

                  agGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACt                                                                       >  2:7809174/1‑93 (MQ=255)

                  agGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATACCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACt                                                                       >  1:136666/1‑93 (MQ=255)

                    ggCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCa                                                                     <  1:3495944/93‑1 (MQ=255)

                       aGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAACCCCCGGCAATCGCCTCACCCGGCACTCCCTGTCTGCTAAAAATGTGACTCaaaa                                                                  <  2:1067103/93‑1 (MQ=255)

                        gCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCaaaaa                                                                 <  1:2580509/93‑1 (MQ=255)

                        gCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCaaaaa                                                                 <  1:6277952/93‑1 (MQ=255)

                         cGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAAc                                                                <  2:3438378/93‑1 (MQ=255)

                          gATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAAcc                                                               >  2:4761536/1‑93 (MQ=255)

                           aTCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAAccc                                                              <  1:1433953/93‑1 (MQ=255)

                            tCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCt                                                             <  1:4400200/93‑1 (MQ=255)

                            tCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCt                                                             >  1:2005939/1‑93 (MQ=255)

                             cacaACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCtt                                                            <  2:1966665/93‑1 (MQ=255)

                               caACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTg                                                          <  2:7433524/93‑1 (MQ=255)

                                  cGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCg                                                       <  1:7554382/93‑1 (MQ=255)

                                    aGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGa                                                     >  1:6387111/1‑93 (MQ=255)

                                    aGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGa                                                     <  2:1087397/93‑1 (MQ=255)

                                    aGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGa                                                     <  2:3872924/93‑1 (MQ=255)

                                     gcgcGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGAt                                                    <  1:2422406/93‑1 (MQ=255)

                                      cgcgGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATg                                                   >  2:5882554/1‑93 (MQ=255)

                                        cgGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATggc                                                 <  2:2253051/93‑1 (MQ=255)

                                         ggCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATggcg                                                <  1:3651101/93‑1 (MQ=255)

                                          gCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATggcgg                                               <  1:6861184/93‑1 (MQ=255)

                                          gCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATggcgg                                               >  2:7068179/1‑93 (MQ=255)

                                            tCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGcc                                             <  2:7918824/93‑1 (MQ=255)

                                            tCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGcc                                             >  1:671563/1‑93 (MQ=255)

                                              cAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCa                                           >  1:7366106/1‑93 (MQ=255)

                                              cAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCa                                           <  1:4197793/93‑1 (MQ=255)

                                              cAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCa                                           <  1:6678791/93‑1 (MQ=255)

                                               aGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAg                                          >  1:2325605/1‑93 (MQ=255)

                                                gTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGc                                         >  2:237238/1‑93 (MQ=255)

                                                 tAAGCAGTAAAATCCCCGGCAATCGCCTCACACGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCa                                        <  2:3769363/93‑1 (MQ=255)

                                                   aGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATc                                      <  2:7732338/93‑1 (MQ=255)

                                                   aGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATc                                      <  2:3080238/93‑1 (MQ=255)

                                                     cAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTg                                    <  1:2164994/93‑1 (MQ=255)

                                                       gTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGtt                                  >  1:1747310/1‑93 (MQ=255)

                                                       gTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGtt                                  >  2:5226035/1‑93 (MQ=255)

                                                          aaaTCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTAc                               <  2:7596291/93‑1 (MQ=255)

                                                           aaTCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTAca                              >  1:1022500/1‑93 (MQ=255)

                                                           aaTCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTAca                              <  1:6143481/93‑1 (MQ=255)

                                                            acccccGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAACCCCTTTGCCGGATGGCGGCCCAGCATCTGTTTAcat                             <  2:6385506/91‑1 (MQ=255)

                                                            aTCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTAcat                             >  1:581955/1‑93 (MQ=255)

                                                             tCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTAcatc                            <  2:7396010/93‑1 (MQ=255)

                                                              ccccGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTAcatca                           >  2:885539/1‑93 (MQ=255)

                                                              ccccGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTAcatca                           >  2:3512035/1‑93 (MQ=255)

                                                               cccGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTAcatcat                          <  2:7246191/93‑1 (MQ=255)

                                                               cccGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTAcatcat                          >  1:4987703/1‑93 (MQ=255)

                                                                ccGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCAtt                         <  1:6549891/93‑1 (MQ=255)

                                                                ccGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCAtt                         >  2:6009227/1‑93 (MQ=255)

                                                                 cGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTa                        <  1:7464881/93‑1 (MQ=255)

                                                                 cGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTATTTACATCATTa                        <  2:5127715/93‑1 (MQ=255)

                                                                  ggCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTac                       <  2:1320313/93‑1 (MQ=255)

                                                                  ggCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTac                       >  2:4560442/1‑93 (MQ=255)

                                                                  ggCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTac                       <  2:5628423/93‑1 (MQ=255)

                                                                   gCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTacg                      <  1:820316/93‑1 (MQ=255)

                                                                     aaTCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTacgac                    <  1:4383716/93‑1 (MQ=255)

                                                                     aaTCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTacgac                    <  2:3858680/93‑1 (MQ=255)

                                                                     aaTCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTacgac                    <  2:5717447/93‑1 (MQ=255)

                                                                      aTCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAACCCCTTTGCCGGATGGCGGCCCAGCACCTGTTTACATCATTacgacg                   <  2:434894/93‑1 (MQ=255)

                                                                      aTCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTacgacg                   >  2:3193087/1‑93 (MQ=255)

                                                                      aTCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTacgacg                   >  2:6640820/1‑93 (MQ=255)

                                                                       tCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGt                  <  2:2864034/93‑1 (MQ=255)

                                                                         gCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCa                >  1:5587762/1‑93 (MQ=255)

                                                                          ccTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCaa               <  2:6340629/93‑1 (MQ=255)

                                                                          ccTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCaa               >  2:1273865/1‑93 (MQ=255)

                                                                          ccTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCaa               <  1:1311653/93‑1 (MQ=255)

                                                                           cTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCaaa              >  2:4844602/1‑93 (MQ=255)

                                                                           cTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCaaa              >  2:2451237/1‑93 (MQ=255)

                                                                           cTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCaaa              <  1:3347330/93‑1 (MQ=255)

                                                                           cTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCaaa              <  1:2290328/93‑1 (MQ=255)

                                                                            tCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGTCGTCAAAc             <  2:6235719/93‑1 (MQ=255)

                                                                             cACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACt            <  2:6488395/93‑1 (MQ=255)

                                                                             cACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACt            <  1:1685614/93‑1 (MQ=255)

                                                                               cccGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGc          <  1:5105649/93‑1 (MQ=255)

                                                                               cccGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGc          >  2:4149029/1‑93 (MQ=255)

                                                                                    cACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTCCTg     <  2:6965694/93‑1 (MQ=255)

                                                                                    cACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTCCTg     <  1:591786/93‑1 (MQ=255)

                                                                                    cACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTCCTg     <  1:7687589/93‑1 (MQ=255)

                                                                                     aCTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTCCTgg    >  1:5842610/1‑93 (MQ=255)

                                                                                      cTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTCCTGGt   <  1:5331705/93‑1 (MQ=255)

                                                                                      cTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTCCTGGt   >  1:3165256/1‑93 (MQ=255)

CGCAGGCCGCTAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCTTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTCCTGGTT  >  NC_000913/3393816‑3393995

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