BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,857,508	T→C	E79G (GAA→GGA)	yidK ←	putative cotransporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,857,508	0	T	C	99.2%	1167.2 / ‑2.7	373	E79G (GAA→GGA)	yidK	putative cotransporter
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (172/197); total (173/200)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 395 total aligned reads displayed.

AGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAA  >  NC_000913/3857418‑3857597

aGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGGCACCGCTGAAGTCACTccc                                                                                         >  2:7407462/1‑93 (MQ=255)

 gAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTcccc                                                                                        <  1:5651416/93‑1 (MQ=255)

 gAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTcccc                                                                                        <  1:1454021/93‑1 (MQ=255)

  aaaaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCa                                                                                       <  1:1327069/93‑1 (MQ=255)

   aaaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAg                                                                                      >  2:334672/1‑93 (MQ=255)

     aTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGcc                                                                                    >  1:6183795/1‑93 (MQ=255)

      tCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGccc                                                                                   <  1:5456216/93‑1 (MQ=255)

      tCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGccc                                                                                   >  2:2919394/1‑93 (MQ=255)

       cGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCa                                                                                  <  2:3541583/93‑1 (MQ=255)

       cGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCa                                                                                  >  2:370276/1‑93 (MQ=255)

        ggggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCAt                                                                                 <  2:7842926/93‑1 (MQ=255)

         gggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATc                                                                                >  2:1641443/1‑93 (MQ=255)

         gggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATc                                                                                >  2:4057424/1‑93 (MQ=255)

         gggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATc                                                                                <  1:3129258/93‑1 (MQ=255)

         gggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATc                                                                                >  2:4967151/1‑93 (MQ=255)

          ggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCa                                                                               >  1:7743581/1‑93 (MQ=255)

          ggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCa                                                                               <  1:2558938/93‑1 (MQ=255)

           gATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCAc                                                                              <  1:983811/93‑1 (MQ=255)

           gATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCAc                                                                              <  1:5222627/93‑1 (MQ=255)

           gATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCAc                                                                              <  1:7530972/93‑1 (MQ=255)

           gATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCAGAAGTCACTCCCCAGCCCATCAc                                                                              >  2:5392194/1‑93 (MQ=255)

             tggtggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCg                                                                            >  1:572680/1‑93 (MQ=255)

               gtggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGAc                                                                          >  2:4789385/1‑93 (MQ=255)

                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  1:3025464/93‑1 (MQ=255)

                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  2:1808044/93‑1 (MQ=255)

                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  1:3207276/93‑1 (MQ=255)

                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  1:1496257/93‑1 (MQ=255)

                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  1:5857690/93‑1 (MQ=255)

                    aaTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGcc                                                                     >  2:4211161/1‑93 (MQ=255)

                    aaTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGcc                                                                     >  1:4963776/1‑93 (MQ=255)

                     aTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCg                                                                    >  1:620071/1‑93 (MQ=255)

                      tGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGc                                                                   >  2:6941770/1‑93 (MQ=255)

                      tGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGc                                                                   <  1:6947764/93‑1 (MQ=255)

                         cgcgCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCttt                                                                <  1:3761196/93‑1 (MQ=255)

                          gcgcTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCtttt                                                               >  2:3791040/1‑93 (MQ=255)

                              ttCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAg                                                           >  2:2850240/1‑93 (MQ=255)

                              ttCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAg                                                           <  1:7353878/93‑1 (MQ=255)

                               tCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAgg                                                          <  2:5817165/93‑1 (MQ=255)

                                cAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGc                                                         >  1:2018405/1‑93 (MQ=255)

                                 aGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGcc                                                        >  1:1412704/1‑93 (MQ=255)

                                  gATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCt                                                       >  2:2548894/1‑93 (MQ=255)

                                  gATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCt                                                       >  1:1808064/1‑93 (MQ=255)

                                   aTAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTg                                                      >  1:2581288/1‑93 (MQ=255)

                                    tAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGc                                                     >  2:83618/1‑93 (MQ=255)

                                    tAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGc                                                     <  2:2156256/93‑1 (MQ=255)

                                     aGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGcc                                                    >  1:4799098/1‑93 (MQ=255)

                                      gcgcGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGccc                                                   <  1:7405596/93‑1 (MQ=255)

                                        gcgGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCgg                                                 >  1:116050/1‑93 (MQ=255)

                                        gcgGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCgg                                                 <  1:5983128/93‑1 (MQ=255)

                                          ggTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaa                                               >  2:5090050/1‑93 (MQ=255)

                                           gTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaaa                                              <  2:2481524/93‑1 (MQ=255)

                                           gTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaaa                                              <  1:3497454/93‑1 (MQ=255)

                                           gTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaaa                                              <  2:1499190/93‑1 (MQ=255)

                                             aaaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAgg                                            >  2:6041631/1‑93 (MQ=255)

                                               aaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGcc                                          <  1:3408242/93‑1 (MQ=255)

                                               aaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGcc                                          >  2:1580145/1‑93 (MQ=255)

                                               aaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGcc                                          <  2:6338445/93‑1 (MQ=255)

                                                    ttAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCa                                     >  1:3922735/1‑93 (MQ=255)

                                                     tAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAg                                    <  1:3292555/93‑1 (MQ=255)

                                                        cgcgAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAgttg                                 <  1:6404525/93‑1 (MQ=255)

                                                          cgAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAgttgtt                               >  1:279859/1‑93 (MQ=255)

                                                          cgAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAgttgtt                               >  2:1295619/1‑93 (MQ=255)

                                                           gAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTc                              >  2:2734491/1‑93 (MQ=255)

                                                            aGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTcc                             <  2:108956/93‑1 (MQ=255)

                                                            aGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTcc                             >  2:7754358/1‑93 (MQ=255)

                                                            aGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTcc                             <  1:2332159/93‑1 (MQ=255)

                                                              gAAGATTAGCGTCACCGCTGAAGTCACTCCCCAGCCCATTACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCCGTTGTTCCgg                           >  2:2024259/1‑92 (MQ=255)

                                                              gAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGt                           <  2:6707901/93‑1 (MQ=255)

                                                              gAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGt                           >  2:7215674/1‑93 (MQ=255)

                                                              gAAGATCAGCGTCACCGCTGAAGTCACCCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGt                           <  1:2933619/93‑1 (MQ=255)

                                                                aGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTgg                         <  1:167970/93‑1 (MQ=255)

                                                                 gATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGa                        >  1:4683642/1‑93 (MQ=255)

                                                                 gATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGa                        >  2:3165820/1‑93 (MQ=255)

                                                                  aTCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGaa                       >  1:2492729/1‑93 (MQ=255)

                                                                   tCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGaaa                      >  2:451988/1‑93 (MQ=255)

                                                                   tCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGaaa                      <  2:7674560/93‑1 (MQ=255)

                                                                    cAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAg                     <  2:3906870/93‑1 (MQ=255)

                                                                    cAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAg                     <  2:4448401/93‑1 (MQ=255)

                                                                     aGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAgg                    <  2:1940676/93‑1 (MQ=255)

                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 <  1:4324291/93‑1 (MQ=255)

                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 <  1:4416453/93‑1 (MQ=255)

                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 <  1:5033471/93‑1 (MQ=255)

                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 >  1:7813217/1‑93 (MQ=255)

                                                                           aCCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggctggtt              >  2:4557554/1‑93 (MQ=255)

                                                                            ccGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTa             <  1:6800480/93‑1 (MQ=255)

                                                                            ccGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTa             <  2:232072/93‑1 (MQ=255)

                                                                            ccGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGGTGTTCCGTGGAAAGGTTGGTTa             >  2:5891660/1‑93 (MQ=255)

                                                                               cTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCa          <  1:4131884/93‑1 (MQ=255)

                                                                                tGAAGTCACTCCCCAGCCCATCAGCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAt         <  1:837973/93‑1 (MQ=255)

                                                                                tGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAt         >  2:5216515/1‑93 (MQ=255)

                                                                                tGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAt         >  1:614588/1‑93 (MQ=255)

                                                                                 gAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAtt        <  1:6500130/93‑1 (MQ=255)

                                                                                  aaGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTa       <  1:2226039/93‑1 (MQ=255)

                                                                                   aGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTaa      <  2:5562852/93‑1 (MQ=255)

                                                                                   aGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTaa      >  1:7122681/1‑93 (MQ=255)

                                                                                   aGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTaa      <  1:102477/93‑1 (MQ=255)

                                                                                    gTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAc     >  2:5686564/1‑93 (MQ=255)

                                                                                    gTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAc     <  1:186206/93‑1 (MQ=255)

                                                                                      cACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGa   <  2:7530405/93‑1 (MQ=255)

                                                                                       aCTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGaa  <  2:4328171/93‑1 (MQ=255)

AGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAA  >  NC_000913/3857418‑3857597

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: `‑`