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%	1192.9 / ‑5.0	364	intergenic (‑4243/+57)	tldD/cafA	TldD/CafA
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (179/184); total (179/185)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

Only 100 of 373 total aligned reads displayed.

CGCAGGCCGCTAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCTTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTCCTGGTT  >  NC_000913/3393816‑3393995

 gCAGGCCGCTAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGt                                                                                        <  2:7174428/93‑1 (MQ=255)

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

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

  cAGGCCGCTAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTc                                                                                       >  1:2912112/1‑93 (MQ=255)

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

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

          tAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTaaaa                                                                               >  1:5839862/1‑93 (MQ=255)

          tAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTaaaa                                                                               >  2:4718519/1‑93 (MQ=255)

          tAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTaaaa                                                                               >  1:347992/1‑93 (MQ=255)

                gcagGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGa                                                                         <  2:10034536/93‑1 (MQ=255)

                gcagGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGa                                                                         <  2:2813320/93‑1 (MQ=255)

                 cagGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGAc                                                                        >  1:2200804/1‑93 (MQ=255)

                   gggCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTc                                                                      >  1:5162011/1‑93 (MQ=255)

                   gggCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTc                                                                      <  2:1363932/93‑1 (MQ=255)

                   gggCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTc                                                                      >  2:792437/1‑93 (MQ=255)

                   gggCAGAGAACACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTc                                                                      <  1:7064337/93‑1 (MQ=255)

                     gCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACCCCCTGTCCGCTGAAAATCTGACTCaa                                                                    >  1:3312811/1‑93 (MQ=255)

                      cAGCGATCACAACTAGTGCGGCTCTAGTAAGCAGTAAAATCTCCTTCAATCGCCTCACCCGTCACGCCCTGTCTGCTTAAAATGTGACTCaaa                                                                   >  2:8950305/1‑93 (MQ=255)

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

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

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

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

                         cGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCaacaac                                                                >  2:1825497/1‑93 (MQ=255)

                         cGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAAc                                                                >  1:9367632/1‑93 (MQ=255)

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

                          gATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAAcc                                                               <  2:3732572/93‑1 (MQ=255)

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

                           aTCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAAccc                                                              >  1:469730/1‑93 (MQ=255)

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

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

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

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

                              acaACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCttt                                                           <  2:2602707/93‑1 (MQ=255)

                               caACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTg                                                          >  1:3999180/1‑93 (MQ=255)

                                aaCGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGc                                                         >  1:641927/1‑93 (MQ=255)

                                aaCGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGc                                                         <  2:6744760/93‑1 (MQ=255)

                                 aCGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGcc                                                        <  1:9896474/93‑1 (MQ=255)

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

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

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

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

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

                                      cgcgGCGCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATg                                                   <  1:8495669/93‑1 (MQ=255)

                                       gcgGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATgg                                                  <  2:9827206/93‑1 (MQ=255)

                                         ggCTCCAGTAAGCAGTATAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATggcg                                                >  2:8341671/1‑93 (MQ=255)

                                         ggCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATggcg                                                >  2:3163211/1‑93 (MQ=255)

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

                                           cTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATggcggc                                              >  2:803015/1‑93 (MQ=255)

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

                                             ccAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGccc                                            >  2:9552032/1‑93 (MQ=255)

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

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

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

                                                   aGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATc                                      >  1:4252383/1‑93 (MQ=255)

                                                     cAGTAAAAGCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAACCCCTTTGCCGGATGGCGGCCCAGCATCTg                                    <  1:5749216/93‑1 (MQ=255)

                                                      aGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGt                                   >  2:5510904/1‑93 (MQ=255)

                                                      aGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGt                                   >  2:5284497/1‑93 (MQ=255)

                                                      aGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGt                                   <  2:8155546/93‑1 (MQ=255)

                                                       gTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGtt                                  <  1:6892123/93‑1 (MQ=255)

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

                                                       gTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGtt                                  <  1:3891275/93‑1 (MQ=255)

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

                                                       gTAAAAACCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGtt                                  <  1:1411465/93‑1 (MQ=255)

                                                        tAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGttt                                 >  1:5592069/1‑93 (MQ=255)

                                                        tAAAATCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGttt                                 >  1:6224662/1‑93 (MQ=255)

                                                         aaaaTCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTa                                >  2:7504195/1‑93 (MQ=255)

                                                            aTCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTAcat                             <  1:2643757/93‑1 (MQ=255)

                                                             tCCCCGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTAcatc                            >  1:2462481/1‑93 (MQ=255)

                                                              ccccGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTAcatca                           <  1:9522852/93‑1 (MQ=255)

                                                              ccccGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTAcatca                           <  2:3034559/93‑1 (MQ=255)

                                                              ccccGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTAcatca                           <  1:7659653/93‑1 (MQ=255)

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

                                                                 cGGCAATCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTa                        >  2:9764967/1‑93 (MQ=255)

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

                                                                       tCGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGTCCAGCATCTGTTTACATCATTACGAAGt                  >  1:6345123/1‑93 (MQ=255)

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

                                                                        cGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATTGCGGCCCAGCATCTGTTTACATCATTACGACGTc                 <  1:7323061/93‑1 (MQ=255)

                                                                        cGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTc                 <  2:8487895/93‑1 (MQ=255)

                                                                        cGCCTCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTc                 <  1:5289472/93‑1 (MQ=255)

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

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

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

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

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

                                                                            tCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAGCCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAAc             <  1:156711/93‑1 (MQ=255)

                                                                            tCACCCGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAAc             <  2:1616568/93‑1 (MQ=255)

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

                                                                                ccGTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCt         >  2:3257137/1‑93 (MQ=255)

                                                                                  gTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTcc       <  1:3035589/93‑1 (MQ=255)

                                                                                  gTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTcc       <  1:1904636/93‑1 (MQ=255)

                                                                                  gTCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTcc       >  2:2112266/1‑93 (MQ=255)

                                                                                   tCACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTCCt      >  1:432341/1‑93 (MQ=255)

                                                                                    cACTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTCCTg     >  2:448875/1‑93 (MQ=255)

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

                                                                                     aCTCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTCCTgg    <  1:1233927/93‑1 (MQ=255)

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

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

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

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

                                                                                       tCCCTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTCCTGGtt  >  2:6215236/1‑93 (MQ=255)

CGCAGGCCGCTAACCAGCAGGGCAGCGATCACAACGAGCGCGGCTCCAGTAAGCAGTAAAATCCCCGGCAATCGCCTCACCCGTCACTCCTTGTCTGCTAAAAATGTGACTCAAAAACCCTTTGCCGGATGGCGGCCCAGCATCTGTTTACATCATTACGACGTCAAACTGCTCCTGGTT  >  NC_000913/3393816‑3393995

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